LOAN  DOCUMENT 


PHOTOGRAPH  THIS  SHEET 


LEVEL 

DOCUMENT  IDENTIFICATION 

M*Y  % 


distribution  statement  a 

Api PTPY®d  for  Public  Release" 
Distribution  Unlimited 


DISTRIBUTION  STATEMENT 


iTTPHIWflM  - 
NTI*  6101 

DTIC  TRAC 

UNANNOUNCED  Q 

JUSTIFICATION 


BY 


DISTRIBUTION/ 


AVAILABILITY  CODES 


DISTRIBUTION 


AVAILABILITY  AND/OR  SPECIAL 


DATE  ACCESSIONED 


A-/  I  I 

DISTRIBUTION  STAMP 


DATE  RETURNED 


1  99905  1  7047 

DATE  RECEIVED  IN  DTIC 


REGISTERED  OR  CERTIFIED  NUMBER 


PHOTOGRAPH  THIS  SHEET  AND  RETURN  TO  DTIC-FDAC 


DTIC 


FORM 

JUNtO 


70A 


DOCUMENT  PROCESSING  SHEET 

LOAN  DOCUMENT 


STOCK  IS  EXHAUSTED. 


w » > n  w r  ©  2 > a 


AD 


Army 
Materiel 
Systems 
Analysis 
Activity 


PROCEEDINGS 

FIRST  MEETING  OF  SPECIAL  INTEREST  GROUP 
ON  CONTROL  THEORY 

22  May  1980 

ABERDEEN  PROVING  GROUND,  MD  21005 


ill! 


Approved  for  public  release 
distribution  unlimited. 


U  S  ARMY  MATERIEL  SYSTEMS  ANALYSIS  ACTIVITY 

ABERDEEN  MOVING  GROUND,  MARYLAND  21005 


DISPOSITION 


Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to 
the  originator. 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an 
official  Department  of  the  Army  position  unless  so  specified 
by  other  official  documentation. 


WARNING 

Information  and  data  contained  in  this  document  are  based  on 
the  input  available  at  the  time  of  preparation.  The  results  may  be 
subject  to  change  and  should  not  be  construed  as  representing  the 
DARCOM  position  unless  so  specified. 


TRADE  NAMES 

The  use  of  trade  names  in  this  report  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial 
hardware  or  software.  The  report  may  not  be  cited  for  purposes 
of  advertisement. 


u 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

Proceeding  of  the  First  Meeting  of  the  Special 
Interest  Group  on  Control  Theory  (22  May  80) 
Aberdeen  Proving  Ground,  MD 

5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

Conference 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR(s) 

H.  Cohen  (Chairman) 

8.  CONTRACT  OR  GRANT  NUMBER(s) 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Director 

US  Army  Materiel  Systems  Analysis  Activity 

Aberdeen  Proving  Ground,  MD  21005 

to.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  &  WORK  UNIT  NUMBERS 

DA  Project  No. 

1 R665706M541 

H.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Director 

US  Army  Materiel  Systems  Analysis  Activity 

ATTN:  DRXsy-MP  „ 

12.  REPORT  DATE 

May  1980 

13.  NUMBER  OF  PAGES 

122 _ 

1 V  mJ$98R| D RESSfi f'dinUlit  1  from  Controlline  Office) 

Cdr,  US  Army  Materiel  Development  &  Readiness 

Comd,  5001  Eisenhower  Avenue,  Alexandria,  VA 

22333 

15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

15a.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  this  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 

18.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (’Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

Control  Theory,  Kalman  Filtering,  Man-Model,  Maneuvering  Target,  Fire  Control, 
Missile  Guidance  and  Control 

20.  ABSTRACT  ("Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

Report  documents  papers  presented  at  first  meeting  of  special  interest  group 
on  control  theory  with  emphasis  on  military  systems  and  recommendations  of 
participants. 

DD 


FORM 
1  JAN  73 


1473 


EDITION  OF  I  NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


TABLE  OF  CONTENTS 


PAGE  NO. 

AGENDA  .  i 

DISCUSSION  AND  SUMMARY  .  iii 

RESEARCH  FOR  US  ARMY  ADVANCED  G§C  SYSTEM .  1 

SUB-OPTIMAL  STATE  ESTIMATION  AS  RELATED  TO  PREDICTIVE 

FIRE  CONTROL  SYSTEM  DESIGN  .  9 

DETECTION  AND  INITIATION  OF  FIRING  COMMAND  FOR  AN 
ACCELERATION  PREDICTOR  FIRE  CONTROL  SYSTEM 

ENGAGING  MANEUVERING  TARGETS . 21 

HYBRID  COMPUTER  SIMULATION  OF  COMBAT  TANK  DRIVEN  RETICLE 

FIRE  CONTROL . 47 

ROBUST  AUTOREGRESSIVE  MODELS  FOR  PREDICTING  AIRCRAFT 

MOTION  FROM  NOISY  DATA  . .  59 

A  DESIGN  METHODOLOGY  FOR  ESTIMATORS  AND  PREDICTORS  IN 

FIRE  CONTROL  SYSTEMS .  75 

AN  ADAPTIVE  LEAD  PREDICTION  ALGORITHM  FOR  MANEUVERING.  .  82 

TARGET  ENGAGEMENT 

APPLICATION  OF  DELAY  FEEDBACK  IN  CONTROL  SYSTEM  DESIGN  .  94 


DISTRIBUTION  LIST 


121 


FIRST  MEETING 


SPECIAL  INTEREST  GROUP  ON  CONTROL  THEORY 

22  May  1980 

MTD  Bldg  400 
APG,  MD  2100S 

AGENDA 


WELCOMING  STATEMENT  -  Keith  A.  Myers 

Assistant  Director,  US  Army  Materiel  Systems  Analysis 
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Aberdeen  Proving  Ground,  MD 

Robust  Autoregressive  Models  for  Predicting  Aircraft 
Motor  From  Noisy  Data 
by  Walter  Dziwak 

US  Army  Armament  Research  $  Development  Command 
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A  Design  Methodology  for  Estimators  and  Predictors  in 
Fire  Control  System 
by  James  F.  Leathrum 
Consultant  to  AMSAA 


Design  §  Implementation  of  Modem  Control  Concept 
to  Improve  the  XM97  Helicopter  Turret  System 
Doo  J.  Lee 

US  Army  Armament  Research  5  Development  Command 
Dover,  NJ  (Paper  not  available  at  this  time) 

Future  Plans  -  Discussion 
Herbert  E.  Cohen 

DARCOM  Mathematics  Program  Office 

US  Army  Materiel  Systems  Analysis  Activity 

Aberdeen  Proving  Ground,  MD 


DISCUSSION  AND  SUMMARY 


The  chairman  of  the  special  interest  group  on  control  theory  solicited 
from  those  attending  their  views  on  the  following  issues,  namely, 

a.  How  can  we  trigger  a  greater  DARCOM-wide  participation  to  this 
important  subject  to  the  Army. 

b.  Was  this  first  meeting  constructive  and  useful.  Was  there  a 
need  for  another  meeting*  If  so,  when,  where  and  who  should  attend. 

All  participants  considered  the  first  meeting  very  useful  addressing 
significant  problem  areas  in  fire  control  and  missile  guidance.  There 
was  unanimous  agreement  to  hold  a  second  meeting  in  early  FY81  (Nov/Dec  80) 
at  APG  and  general  agreement  to  invite  the  Navy  and  Air  Force  to  participate. 
Recommendation  was  made  and  favorably  indorsed  by  all  to  change  the  name  of 
the  group  to  "Coordinating  Group  on  Modern  Control  Theory"  to  emphasize 
the  application  of  "modem  control  theory  and  techniques"  to  military 
problems  as  compared  to  classical  control  theory  techniques.  Suggestion  was 
made  that  the  chairman  explore  the  feasibility  of  using  the  JTCG/ME  by  which 
the  Navy  and  Air  Force  could  participate. 

It  was  strongly  felt  that  travel  funds  severely  limited  number  of 
participants  and  that  DARCOM  HQ  support  should  be  solicited  to  encourage 
subordinate  commands  to  actively  participate  in  this  coordinating  group. 
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RESEARCH  FOR  US  ARMY  ADVANCED  G&C  SYSTEM 


Dr.  Harold  L.  Pastrick 
Guidance  &  Control  Directorate 
US  Army  Missile  Laboratory 
Redstone  Arsenal,  Alabama  35809 


Abstract 

In  1973  the  US  Army  Missile  Command  embarked 
upon  a  task  to  develop  an  advanced  guidance  and  con¬ 
trol  system  for  future  missiles.  It  was  intended 
to  "leapfrog"  systems  currently  under  development  in 
order  to  meet  the  stringent  demands  and  constraints 
imposed  by  targets  with  predicted  characteristics  of 
the  1990s  and  beyond  and  by  the  predicated  battle¬ 
field  environment  of  that  time.  * 

In  this  paper  the  problem  is  redefined  and  the 
latest  development  program  presented.  Results  that 
have  been  achieved  to  date  are  described,  particu¬ 
larly  in  the  areas  of  mathematical  models  of  the 
missiles  and  their  guidance  and  control  (G&C)  sys¬ 
tems  being  used  for  analyses  and  simulations;  aero¬ 
dynamics;  propulsion;  guidance  laws  being  developed 
and  analyzed;  status  of  development  of  Disturbance 
Accommodating  Control;  signal  processing  to  locate 
and  track  the  target(s);  and  digital  design  tool 
development.  The  paper  is  concluded  with  a  section 
of  future  plans,  to  include  contractor  support. 

I.  Introduction 

The  US  Army  Missile  Command  (MICOM)  recently 
began  a  task  to  develop  an  advanced  G&C  system  for 
Future  Army  Modular  Missiles.  The  intent  is  to 
"leapfrog"  systems  currently  under  development. 

The  purpose  of  this  paper  is  to  describe  the  work 
that  has  been  completed  within  this  new  task  and  to 
provide  an  indication  of  future  efforts  that  are 
now  planned.- 

The  first- step  In  implementing  this  task  was 
to  conduct  a  literature  survey  to  establish  a  tech¬ 
nology  base  starting  point.  Following  this  survey, 
guidance  laws  were^placed  In  five  categories  and 
defined  mathematically.  The  implementation  and 
predicted  performance  of  each  category  was  then  in¬ 
vestigated  and  compared  in  light  of  current  and 
predicted  hardware  and  software  capabilities.  ‘ 

This  work  was  subsequently  updated  in  1979.  2 

The  program  objectives  are  three-fold; 

•  To  develop  and  prove  a  G&C  system  that  is 
capable  of  guiding  and  controlling  future  US  Army 
missiles  (generally  defined  as  air  defense  and  sur- 
face-to-surface  general  support)  to  destroy  pre-? 
scribed  lines  of  future  targets ,,  This  must  be 
accomplished  under  the  predicted  severe  battlefield 
environment  of  the  future. 

•  To  "leapfrog"  systems  under  current  develop¬ 
ment. 


r- his  paper  is  declared  a  work  of  the  US  Govern¬ 
ment  and  therefore  is  in  the  public  domain. 
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•  To  broaden  and  deepen  the  existing  G&C  sys¬ 
tem  technology  and  design  base  within  the  G&C 
Directorate  of  MICCM. 

In  the  sequel  the  expected  threat  and  its 
characteristics  will  be  summarized.  This  is  follow¬ 
ed  by  a  description  of  the  development  plan  for  the 
advanced  G&C  system.  The  progress  to  date  is  pre¬ 
sented,  and  the  paper  is  concluded  with  a  section 
on  future  plans. 

II .  Expected  Threat 

Theater  defense  typically  is  provided  by  a 
mixture  of  ground-based  and  airborne  defense  systems 
supported  by  radars,  command  and  control  systems, 
electronic  warfare  equipment,  and  passive  measures 
such  as.  camouflage,  decoys,  and  equipment  disper¬ 
sion.  The  air  defense  objective  of  ground  based 
systems  is  to  limit  the  opponent's  effectiveness  by 
attacking  his  critical  assets  so  that  land  forces 
may  maneuver  with  a  minimum  of  interference  from 
the  enemy  air  weaponry. 

For  many  years  now,  the  enemy  doctrine  has 
emphasized  large  mass  and  brute  force,  and  his  air 
attacks  will  provide  no  exception.  It  is  entirely 
feasible  to  assume  that  an  attack  in  the  Central 
Europe  area  will  oe  accompanied  by  several  thousand 
combat  aircraft.  2  In  addition,  his  doctrine  calls 
for  the  massing  of  large  quantities  of  artillery 
fire  on  a  section  selected  for  a  tank-led  break¬ 
through.  It  is  unlikely  that  NATO  forces  either 
now  or  in  the  near  future  will  match  the  Warsaw 
Pact  forces  in  terms  of  numbers  of  weapons,  nor  is 
it  the  intent  to  aim  toward  that  end.  Rather,  it  is 
important  to  optimize  the  effectiveness  of  our* 
smaller  force  to  meet  the  anticipated  threat. 

The  Army  is  attempting  to  maximize  the  effec¬ 
tiveness  of  its  current  family  of  air  defense 
weapons  while  concurrently  developing  a  new  family 
to  meet  the  threat  of  the  I990's.  In  the  near 
term,  there  will  be  continued  modification  of  cur¬ 
rent  systems  as  necessary,  and  while  still  feasible, 
to  overcome  qualitative  and  quantitative  deficien¬ 
cies.  Longer-term  replacements  continue  in  devel¬ 
opment  or  procurement  for  all  the  major  field  army 
air  defense  systems.  Examples  of  this  strategy 
include  the  following;  high  to  medium  altitude 
missile  systems  -  PATRIOT  for  NIKE  HERCULES  and 
HAWK;  short  range  miss il e  systms  -  U.  S.  ROLAND  for 
CHAPARRAL;  manportable  missiles’-  STINGER  for 
REDEYE;  mobile  gun  systems  -  DIYAO  given  to  VULCAN. 
The  systems  will  provide  the  effective  aerial 
umbrella  needed  by  our  forces  to  not  only  survive, 
but  to  fight  effectively. 

For  security  reasons,  it  is  impossible  in  this 
forum  to  describe  specifically  the  air  threat  that 
will  be  encountered  in  the  scenario  described  above. 
However,  in  order  to  Quantify  the  problem  somewnat, 
we  shall  attempt  to  attribute  vehicle  characteris¬ 
tics  to  the  enemy  based  on  our  current  technology 
in  the  field  of  air  defense  targets.  The  Targets 
Management  Office,  US  Army  Missile  Command  since 


1 


1964  has  published  an  extensive  library  of  target 
program  reports.  In  particular,  they  have  classi¬ 
fied  a  variety  of  aerial  targets  as  test  and  evalua¬ 
tion  (T&E)  targets  used  for  air  defense  weapon 
systems.  Needless  to  say,  targets  used  for  this 
purpose  must  exercise  an  air  defense  weapon  system 
to  the  limits  of  its  capabilities.  5 

A  particularly  Interesting  target  Is  known  as 
HAHST,  an  acronym  for  High  Altitude  High  Speed 
Target.  It  is  designed  to  achieve  speeds  up  to 
Mach  4  at  altitudes  up  to  100,000  feet.  Additional 
performance  characteristics  for  HAHST  and  other 
existing  targets  are  given  in  References  2  and  5. 

The  effectiveness  of  any  missile  system  is  con¬ 
ditioned  on  its  ability  to  function  in  an  Electronic 
Countermeasures  environment.  Stand  off  jammers, 
barrage  jammers  and  even  dispensed  chaff  will  be 
used  to  deny  the  air  defense  tracking  radars  the 
capabilities  needed  to  be  effective  in  their  sector. 
The  jammers  are  intended  to  reduce  the  acquisition 
range  of  the  radars  and.  If  perfected,  will  elimi¬ 
nate  accurate  tracking  entirely.  Additionally,  the 
threat,  aircraft  and  missiles  will  be  fabricated  to 
present  the  smallest  possible  radar  cross  section. 
The  state-of-the-art  In  this  field  is  beyond  the 
scope  of  presentation  in  this  paper. 


utilize  different  -nodules  such  as  propulsion,  guid¬ 
ance,  warhead,  etc. 

A  program  plan  initiated  by  the  G&C  Director¬ 
ate,  MICCM,  was  undertaken  approximately  two  years 
ago  with  the  program  objectives  enumerated  above. 
The  program  plan  contains  six  intermediate  objec¬ 
tives,  the  end  of  which  each  constitutes  a  program 
milestone.  They  are: 

§  Define  overall  program; 

•  Collect  elements  that  will  form  the  candi¬ 
date  GiC  systems; 

•  Define  candidate  GiC  systems; 

•  Evaluate  the  candidate  G&C  systems: 
select  best  system; 

•  Design  and  fabricate  the  selected  system; 

•  Demonstrate  "Prcof-of-Concept"  of  the 
selected  system. 

These  separate  program  elements  (or  intermediate) 
objectives  are  described  oelow  in  more  detail.  An 
accompanying  milestone  chart  is  provided  as  Fig.  1. 


From  a  defensive  viewpoint,  the  effect  of 
enemy  jamming  of  the  air  defense  radars  and  the 
minimization  of  enemy  attack  aircraft  radar  cross 
sections  have  a  profound  impact  on  the  air  defense 
missile  system.  The  rationale  Is  reasonable  and 
straightforward.  If  the  enemy  does  indeed  have 
aircraft  and  attack  missiles  and  remote  pilotless 
vehicles  (RPVs)  either  with  or  better  than  the 
characteristics  attributed  to  them  via  the  method 
above,  and  If  the  enemy  minimizes  his  radar  cross 
section  to  the  current  state-of-the-art,  his  attack 
vehicles  will  be  extremely  difficult  to  acquire  at 
long  range.  The  effect  of  the  combined  high  speed 
and  high  agility  (i.e.,  high  g-maneuver  capability) 
with  low  radar  cross  section  yields  precious  little 
reaction  time  to  the  air  defense  system.  The  close- 
in  acquisition  will  seriously  degrade  the  existing 
air  defense  missile's  G&C  system's  performance, 
since  most  are  based  on  a  proportional  navigation 
and  guidance  (PNG)  law.  An  environment  such  as 
presented  above,  however,  can  be  better  addressed 
in  terms  of  guidance  laws  explicitly  tailored  to 
this  type  of  threat.  Thus,  the  optimally  guided 
and  controlled,  highly  maneuvering,  defensive  mis¬ 
sile  using  terminal  guidance  sensors  chosen  from 
across  a  wide  range  of  the  frequency  spectrum  must 
be  initiated  Into  the  development  cycle.  As  a 
necessary  first  step,  the  research  described  in  the 
sequel  addresses  that  problem. 

III.  Development  Program  Plan 

From  an  overall  systems  viewpoint,  this  pro¬ 
gram  shall  address  the  issue  of  creating  new  theory 
in  the  G&C  area  to  meet  the  high  performance  threat 
of  the  future  as  one  of  the  leading  technology 
items.  Closely  associated  with  it  and  In  parallel 
with  the  G&C  effort,  weapon  system  work  shall  be 
undertaken  to  modify  airframe  and  propulsion  to  be 
capable  of  engaging  the  threat  of  the  1990's.  Gen¬ 
eral  support  weapons  shall  be  viewed  initially  as  a 
subset  of  the  air  defense  system(s),  whereas  pre¬ 
viously,  these  two  classes  of  weapons  were  develop¬ 
ed  independently.  This  research  shall  attempt  to 
view  them  as  potentially  similar  systems  that 


NO.  TASK 

1  DEFINE  PROGRAM  -a 

2  COLLECT  G&C  ELEMENTS 

3  DEFINE  CANDIDATE  SYSTEMS 

4  EVALUATE  CANDIDATES  I  ..  j 

5  DESIGN  &  FABRICATE  C 

6  PROOF  OF  CONCEPT 


—  -  1  fc  L  ■  *  * . — 

FY79  FY30  FY81  FY82  FY83  FY84 
Figure  1.  Program  Milestone  Chart. 

Task  1.  Define  Overall  Program 

In  this  first  task  it  is  necessary  to  define 
the  program  objectives,  goals,  and  constraints. 

This  Is  followed  by  a  survey  of  representative  mis¬ 
sile  plants  and  the  characteristics  of  sensors  and 
effectors  that  might  be  used  (either  available  or 
under  development).  Target  dynamics  and  Initial 
conditions  must  be  defined,  as  must  missions  for 
the  missile  system(s)  using  the  advanced  G&C  sys¬ 
tems  that  emanates  from  this  program.  Also  includ¬ 
ed  in  the  first  program  element  is  the  beginning  of 
coordination  with  other  missile-developing  agencies 
and  services.  Finally,  the  time  frame  predicted 
usage  of  the  G&C  system  must  be  defined.  As  des¬ 
cribed  in  this  paper,  most  of  Task  1  has  been  com¬ 
pleted. 

Task  2.  Collect  Elements  for  Candidate  G&C  Systems 

In  this  task  detailed  definitions  of  the 
characteristics  of  candidate  guidance  laws,  candi¬ 
date  autopilots  {control  laws),  state  estimation 
techniques,  state  truncation  techniques,  and  Dis¬ 
turbance  Accommodation  Control  (DAC)  shall  be 
accomplished.  5  Further,  characteristics  of  all 
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expected  disturbances,  both  external  (to  the  mis¬ 
sile)  and  Internal,  must  be  collected  and  evaluated. 
The  need  for  DAC  and  for  state  estimation  will  be 
evaluated.  Characteristics  of  all  expected  signif¬ 
icant  system  nonl Inearitles  will  be  collected  and 
their  dynamic  Importance  evaluated. 

Task  3.  Define  Candidate  G&C  Systems 

The  components  and  subsystems  that  have  been 
identified  in  the  earlier  tasks;  will  now  be  combin¬ 
ed  Into  candidate  G&C  systems.  Missile  structural 
dynamics  will  be  determined  (if  not  already  com¬ 
pleted  within  the  plant  definition  of  the  earlier 
tasks)  or  refined,  as  necessary.  The  simulation 
program  objectives  will  be  defined;  the  simulation 
program  will  then  be  defined;  and  finally,  develop¬ 
ment  of  the  final  system  simulation  will  begin. 
Investigations  of  the  dynamics  of  the  candidate  G&C 
systems  will  begin  using  both  mathematical  analysis 
and  computer  simulation.  A  figure  of  merit  (cost 
functional)  will  be  developed  during  this  task,  as 
will  model  error  effects  and  criteria. 

Task  4.  Evaluate  Candidate  G&C  Systems 

The  development  of  the  system  simulation  and 
the  investigation  of  the  candidate  G&C  systems, 
begun  In  Task  3,  will  be  completed.  The  candidate 
G&C  systems  then  will  be  evaluated  with  respect  to 
the  Figure  of  Merit  developed  In  Task  3,  using 
mathematical  analysis  and  computer  simulation.  The 
strengths  and  weaknesses  of  distributed  versus  cen¬ 
tralized  controllers  will  be  assessed,  probably  as 
enhanced  by  the  use  of  microcomputers.  Finally, 
the  best  G&C  system  will  be  selected  and  Integrated 
into  a  missile  airframe. 

Task  5.  Design  and  Fabrication 

A  detailed  design  of  the  selected  G&C  system 
will  be  performed.  The  testbed(s)  selected  for  use 
in  the  ?roof-of-Concept  phase  will  be  fabricated 
and/or  assembled. 

Task  6.  Prop f-of -Concept 

The  Proof-of-Concept  will  be  demonstrated  with 
hardware  firings  during  this  task.  These  firings 
will  be  augnented  by  simulations  (computer  and 
hardware-in-the-loop  as  deemed  necessary)  and  anal¬ 
yses  as  necessary. 

IV.  Progress  to  Date 

The  accomplishments  to  date  have  been  achieved 
by  elements  of  MICOM  and  through  the  use  of  re¬ 
search  contracts.  It  has  been  supplemented  through 
coordination  with  the  US  Air  Force  Armaments  Labo¬ 
ratory  at  Eglin  Air  Force  Base,  the  US  Army  Bal¬ 
listic  Missile  Defense  Systems  Command  at  Hunts¬ 
ville,  Alabama,  and  the  Office  of  the  Under  Secre¬ 
tary  of  Defense.  This  coordination  has  been  and 
will  continue  to  be  carried  out  to  eliminate  dupli¬ 
cated  development  effort  on  similar  projects  within 
the  Department  of  Defense.  To  augment  the  in-house 
research  and  engineering  capability,  several  re¬ 
search  contracts  have  been  initiated  in  specialized 
areas  currently  including:  the  Computer  Sciences 
Corporation,  Huntsville,  Alabama;  the  (Dynamic 
Systems  Research  and  Training  Corporation,  Hunts¬ 
ville,  Alabama;  the  University  of  Florida;  Western 
Kentucky  University;  and  the  Control  Dynamics  Com¬ 
pany,  Huntsville,  Alabama.  Additional  contracts 


with  other  organizations  are  anticipated  as  the 
scope  of  the  program  grows. 

The  technical  work  comprising  Tasks  1  and  2 
has  been  apportioned  to  various  members  of  the 
MICOM  and  contractor  team.  As  indicated,  most  of 
Task  1  has  been  completed,  and  effort  on  Task  2 
work  Is  underway. 

A.  Target  Definition 

A  comprehensive  investigation  of  predicted 
future  targets  and  their  dynamics  has  been  complet¬ 
ed.  This  investigation  included  reviewing  and  dis¬ 
cussing  material  available  within  the  sources  of 
the  US  Army  and  US  Air  Force  and  included  inputs 
from  several  industrial  organizations .  Finally, 
the  collected  information  was  discussed  with  ele¬ 
ments  of  the  Office  of  the  Under  Secretary  of 
Defense,  Research  and  Engineering.  While  most  of 
the  collected  information  is  classified  for  secu¬ 
rity  reasons,  the  results  of  this  most  important 
phase  of  the  study  lead  to  emphasizing  targets  in 
three  categories:  highly  maneuverable  aircraft, 
cruise  missiles,  and  tactical  ballistic  missiles. 
Not  addressed  in  this  study  are  RPV's  and  surface 
targets.  It  is  felt  that  use  of  a  sophisticated 
missile  system  to  engage  numerous  RPV's  would  not 
be  cost  effective  and  that  missiles  currently 
under  development  will  be  able  to  combat  surface 
targets.  Except  for  periodic  updating,  this  sub¬ 
task  Is  now  complete. 

B.  Missile  Plant  Characteristics 

In  order  to  get  the  other  subtasks  underway, 
standard  equations  of  motion  for  missile  bodies 
with  some  structural  flexibility  have  been  used. 

The  missiles  are  assumed  to  be  acted  upon  by  the 
usual  aerodynamic  forces  and  torques.  Currently, 
variations  of  the  SPRINT  missile  family  are  among 
the  principal  contenders  for  the  airframe  and  pro¬ 
pulsion  system. 

C.  Aerodynamics 

A  detailed  plan  is  being  generated  for  future 
actions  in  gathering  and  generating  aerodynamic 
data  to  be  used  in  this  program.  Novel  aerodynamic 
shapes  are  under  consideration  and  evaluation  for 
use  in  developing  control  authority  for  missiles 
using  advanced  G&C  systems.  These  include  missiles 
using  configurations  shown  in  Table  1  (along  with 
their  characteristics) . 

0.  Error  Sources 

Expected  error  sources  have  been  categorized 
into  five  detailed  groups,  to  include  predicted 
mean  and  standard  deviation  values.  Several  are 
presented  as  representative  of  the  groups.  These 
will,  of  course,  change  as  the  choices  are  narrowed 
and  the  program  develops. 

1 .  Missil e .  Included  in  this  group  are  para¬ 
meter  uncertainties  such  as  thrust  magnitude  and 
misalignment,  C.G.  location  and  offset,  mass,  and 
transverse  and  axial  moments  of  inertia. 

2.  Aerodynamic.  Normal  and  axial  forces; 
pitch,  yaw,  and  roll  moments;  pitch  and  yaw  damping 
derivatives;  roll  damping  coefficient  due  to  fin 
deflection;  and  axial  drag  are  factors. 
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Table  1.  Aerodynamic  Configurations  and  Characteristics. 


Configuration 

Advantages 

Disadvantages 

Mach 

No. 

Altitude 

(ft.) 

Maneuverability 

(9's) 

1.  Deployable  wings 
with  all -moveable 
tails 

a.  Simple  con¬ 
trols 

b.  Good  roll 
control 

a.  Control  force 

In  opposite  direc¬ 
tion  from  maneuver 

b.  High  angle  of 
attack 

H 

0-100  K 

10-15 

2.  Low  aspect 
ratio,  long-chord 
delta  wing  with 
tab  controls 

a.  Simple  con¬ 
trols 

b.  Good  roll 
control 

a.  Control  force 

In  opposite  direc¬ 
tion  from  maneuver 

b.  High  angle  of 
attack 

4-7 

0-140  K 

20-30 

3.  Flared  skirt- 
stabilized  missile 
with  all  moveable 
wings 

a.  Control  force 
in  direction  of 
maneuver 

b.  Lower  angle 
of  attack 

Possible  higher 
drag  and  hinge 
moment 

H 

0-140  K 

20-30 

4.  Reentry  body 
shape  with  moveable 
wedge  controls  * 

a.  Simple  shape 

b.  Good  roll 
control 

Roll  control  sur¬ 
faces  separate 
from  maneuvering 
controls 

5-10 

7 

? 

5.  Lifting  body  * 

a.  Lower  angle 
of  attack 

b.  Control  force 

In  direction  of 
maneuver 

Complicated  aero¬ 
dynamic  shape 

3-6 

0-90  K 

10 

*  Comment.  Must  use  bank-to-turn  guidance. 


3.  Instruments 

.  .  •  Accelerometers  -  scale  factor  stability, 

bias  stability,  non-orthogonality,  **g2"  .  SCale 
factor,  third-order  scale  factor,  cross-axis  sensi¬ 
tivity,  cross-coupling,  scale  factor  asymmetry,  and 
rectification  error. 

•  Gyroscopes  -  scale  factor  stability; 
bias  stability;  non-orthogonality;  anlsoelastlc 
drift;  drift  rates  In  pitch,  yaw,  and  roll  due  to 
input  axis  and  spin  axis  accelerations  as  well  as 
those  (rates)  Independent  of  acceleration,  due  to 
torquer  nonlinearities,  and  due  to  electronic 
noise;  and  mass  unbalance. 

•  Porro-prlsm  azimuth  alignment. 

•  Laser  Inertial  measurement  unit  mis¬ 
alignment  with  respect  to  missile  body  axes. 

•  Accelerometer  triad  origin  displacement. 

•  Uplink/downlink  bias  and  calibration 

errors . 

•  Optical  correlator  errors . 

•  Radar  -  range  track  and  angle  track 
noise  and  accuracy,  ground  clutter  noise,  and  tar¬ 
get  glint  noise. 

t  Eternal.  Wind  magnitude  and  direction 

(initial  azimuth,  elevation,  and  roll  alignment  of 


the  missile;  target  velocity  and  Illumination 
jitter;  semiactive  laser  pointing  accuracy  and  beam 
divergence;  gravity  bias;  atmospheric  effects 
(e.g.,  upon  radio  range). 

5.  Subsystems 

•  Common  nonlinearities  -  saturation, 
coulomb  (and  other)  friction,  backlash,  and  bang- 
bang  with  dead  zone. 

•  Computer  -  quantization,  truncation,  and 
fixed  word  length. 

•  Seeker  -  boreslght  error  (In  pitch  and 
yaw)  due  to  servo  noise;  channel  crosscoupling; 
coupling  between  the  seeker  head  and  the  airframe; 
effects  of  the  radome  and  irdome  on  angle  linear¬ 
ity;  angle  bias  (the  electrical  equivalent  of  mech¬ 
anical  3SE);  gain  stability;  angle  noise;  and  bore- 
sight  error  in  pitch  and  yaw  due  to  clutter,  re¬ 
ceiver,  and  jamming  noise. 

•  Autopilot  -  bias  errors,  time  delays, 
and  gain  stability. 

•  Guidance  -  errors  In  initial  position, 
velocity,  and  acceleration. 

As  the  program  progresses,  these  error  sources 
will  be  analyzed  further  to  determine  which  might 
be  amenable  to  cancellation  by  appropriate  disturb¬ 
ance  accommodation  design  theory  and  techniques. 


4 


E.  Guidance  Laws 

This  is  a  major  thrust  area  within  the  program. 
As  indicated  above,  an  extensive  literature  search 
has  been  completed  and  documented.  This  was  follow¬ 
ed  by  placing  guidance  laws  in  five  categories  and 
describing  each  mathematically.  The  implementation 
and  predicted  performance  of  each  category  has  been 
initially  reported  in  Reference  1  and  subsequently 
updated  in  Reference  2.  A  summary  of  the  latter 
is  shown  in  Table  2. 

Investigations  continue  into  areas  of  optimal 
guidance,  in  particular,  emphasizing  digital  as¬ 
pects  in  anticipation  of  the  expected  use  of  on¬ 
board  digital  controllers.  Constant  effort  Is  made 
to  reduce  Implementation  complexity,  where  complex¬ 
ity  is  defined  as  requirements  for  hardware  and 
software.  Other  innovative  techniques  are  under 
investigation.  They  include  assessing  the  potential 
application  of  Singular  Perturbation  Theory,  Distur¬ 
bance  Accommodating  Theory,  and  means  of  determining 
or  predicting  the  very  Important  (for  optimal 


applications)  quantity,  time-to-go,  i.e.,  remaining 
time  Of  flight  at  any  instant  (see  Table  3). 

1 .  Optimal  Control .  A  conventional  implementa¬ 
tion  of  Linear-Quadratic  Optimal  Control  is  por¬ 
trayed  graphically  in  Figure  2.  The  optimal  con¬ 
trol  authority,  ujg,  is  selected  to  be 

u°Q  .  -R-ViCxi  .  0) 

The  matric  Riccati  equation  is  solved  to  obtain  the 
control  gain,  ( t ) ,  while  minimizing  a  quadratic 

performance  index. 

2.  Disturbance  Acconmodati nq  Control.  The 
Di sturbance  Ac commoda ting  Control  ( DA C )  t h eo ry  is 
described  for  the  continuous-time  domain  in  Refer¬ 
ence  6.  The  general  nature  of, the  DAC  controller 
is  to  generate  a  real-time  on-line  estimate  of  the 
actual  (instantaneous)  disturbance  waveform  and 
create  a  special  control  action  that  exactly 


Table  2.  Conventional  Guidance 


Approach 

Advantages 

Limitations 

1.  Attitude  Pursuit 

a.  Simplest  implementation 

b.  Fixed  targets 

Sensitive  to  target  velocity, 
disturbances 

2.  Velocity  Pursuit 

a.  Simple  implementation 

b.  Non-maneuvering  targets 

Sensitive  to  target  accelera¬ 
tion,  disturbances 

3.  Proportional  Naviga- 
tion  Guidance  (PNG) 

a.  Simple  implementation 

b.  Maneuvering  targets 

Sensitive  to  high  end  game 
maneuvers 

Table  3.  New  Guidance  Approaches 


Approach 

Advantages 

Limitations 

1.  Linear  Quadratic  (LQ) 
Regulator 

Better  than  PNG  against 
maneuvering  target 

a.  Tqq  estimate  required 

b.  Must  compute  time-varying 
control  gain  (Kx) 

c.  Disturbances  ignored 

2.  Linear  Quadratic 

Gaussian  (LQG) 

Regulator 

Better  than  LQ  against 
noise- type  disturbances 

a.  Tqq  estimate  required 

b.  Must  compute  time-varying 
control  gain  (i^) 

c.  Must  compute  Kalman  gain 
for  estimator 

3.  Disturbance-Utiliz¬ 
ing  Control  (DUC) 

Better  than  LQ  or  LQG 
against  waveform- type 
disturbances 

a.  Tg0  estimate  required 

b.  Must  compute  2  time- vary¬ 
ing  control  gains  (Kx,  Kxz) 

4.  Singular  Perturbations 


Computational  efficiency 


An  approximation  to  optimal 
control 


cancel s-out  the  disturbance  effect  on  the  missile. 
The  DAC  theory  will  be  extended  Into  a  discrete- 
time  domain.  To  date,  the  class  of  systems  and 
disturbances  amenable  to  a  discrete-time  version  of 
DAC  have  been  defined.  Work  is  underway  on  a  des¬ 
cription  of  the  state-reconstructor  that  will  be 
associated  with  a  digital  DAC  controller.  Also 
under  consideration  Is  the  possibility  of  applying 
DAC  to  various  missile  subsystem  or  component  out¬ 
puts,  such  as  sensors,  whose  normal  outputs  have 
been  modified  by  the  Influence  of  the  disturbances 
upon  the  sensors. 

An  Innovative  modification  of  DAC  is  Dis¬ 
turbance-Utilizing  Optimal  Control  (OUC)7.  In  this 
case  waveform-type  disturbances  are  exploited  opti¬ 
mally.  Examples  of  such  disturbances  are  drag, 
target  maneuvers,  wind  gusts,  any  effects  of  the 
gravitational  field.  A  graphical  portrayal  of  DUC 
implementation  Is  shown  In  Figure  3.  In  this  case, 
the  optimal  control  authority  Is  specified  as 


DUC 


-R-V(Kxx 


Kxz‘> 


(2) 


where  Kjt)  Is  found  by  solving  a  matrlc  Riccatl 
equation,  and  K  (t)  Is  found  by  solving  a  linear 
differential  eqQatlon.  To  date,  computer  simula¬ 
tions  have  shown  DUC  to  be  quite  effective  when 
compared  to  the  performance  achieved  by  using  con¬ 
ventional  LQ  controllers. 


3,  Singular  Perturbations  .  Appl  ication  of 
Singular7  Perturbation  Theory  to  missile  control  may 
be  attractive  if  It  is  deemed  necessary  for  the  « 
control  law  to  account  for  high  order  model  terms. 
The  standard  approach  is  to  approximate  the  model 
with  relatively  lower  order  equations.  However, 
the  neglected  higher  order  terms  may  be  dynamically 
significant.  3ecause  their  inclusion  might  create 
computational  problems,  a  possible  alternative 
might  be  the  application  of  Singular  Perturbation 
Theory.  The  application  is  particularly  amenable 
to  controller  design  where  the  open-loop  plant  has 
a  wide  eigenvalue  dispersion,  slow  and  fast  modes, 
or  parasitic  parameters. 

F.  Sensor  Characteristics 


Various  existing  and  predicted  sensors  have 
been  characterized.  Trade  studies  to  aid  In  their 
selection  have  been  identified.  It  now  appears 
that  a  sensor  using  some  form  of  pattern  recogni¬ 
tion  may  be  required.  The  problem  is  to  find  or 
begin  developing  a  sensor  that  can  provide  a  guid¬ 
ance  signal  with  a  superior  signal -to-noise  ratio 
from  data  that  has  been  deliberately  modified,  for 
example,  by  high  powered  jamming  equipment. 

G.  Computer-Aided  Design  Tools 


A  number  of  computer  programs  are  available 
(such  as  root  locus  and  other  frequency  domain 


Disturbance 


Estimate 


Figure  2.  Conventional  Linear  Quadratic  Optimal  Control. 


Disturbance 

w(t) 


Figure  3.  Disturbance-Utilizing  Optimal  Control. 
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techniques)  to  aid  in  design  of  control  systems. 
Most  of  them  are  in  the  continuous-time  domain* 
although  there  does  exist  a  z-transform  manipula¬ 
tion  program.  An  executive  routine  currently  Is 
being  coded  to  manipulate  these  programs  effi¬ 
ciently. 

H.  Microcomputer  State-of-the-Art 

A  continuing  assessment  of  microcomputer  state- 
of-the-art  has  been  instigated  within  this  program. 
Microcomputers  will  become  an  integral  part  of 
advanced  G&C  systems*  with  their  small  size,  light 
weight,  and  relatively  low  cost.  Presently,  micro¬ 
computers  software  capabilities  are  being  Investi¬ 
gated  with  a  Tektronix  8002  Microprocessor  Labo¬ 
ratory.  This  investigation  will  be  followed  by  an 
analysis  of  the  computational  capability  to  meet 
the  requirements  imposed  by  an  advanced  G&C  system. 

I.  Digital  Design  Tool  Development 

The  probability  that  an  advanced  G&C  system 
will  be  implemented  digitally  seems  to  be  nearly 
assured.  It  is  for  this  reason  that  the  foregoing 
work  has  been  oriented  strongly  in  the  direction 
of  digital  Implementation  and  its  consequences.  A 
review  of  analysis  and  design  tools  available  to 
the  G&C  system  designer  Is  being  conducted.  Where 
the  need  for  a  new  tool  appears  to  be  warranted,  It 
will  be  developed  within  this  program  when  possible. 
Emphasis  is  placed  on  both  simplicity  of  applica¬ 
tion  and  being  able  to  draw  on  practising  engineers 
existing  engineering  training  and  experience. 

One  example  of  such  a  tool  is  "SAM"  (acronym 
for  Systematic  Analytical  Method).  SAM  provides  an 
alternative  to  the  use  of  signal  flow  graphs  and 
the  application  of  Mason's  Gain  Rule  to  determine 
selected  states  of  the  missile.9  The  technique  Is 
particularly  useful  for  analysis  of  complicated 
samoled-data  control  systems.  An  advantage  of 
using  SAM  is  that  the  cumbersome  application  of 
Mason's  Gain  Formula  can  be  avoided.  Further,  the 
entire  method  of  constructing  signal  flow  graphs 
may  be  circumvented.  Since  only  the  equations  des¬ 
cribing  the  system  are  needed  for  SAM,  even  the 
customary  block  diagram  is  not  needed.  The  tech¬ 
nique  is  analytical  in  nature  and  makes  use  of  a 
systematic  manipulation  of  the  system  algebraic 
equations.  These  manipulations  follow  prescribed 
rules  set  forth  in  the  technique. 

A  second  simplified  technique  is  the  determina¬ 
tion  of  digital  control  system  response  by  cross  - 
multiplication.!0  This  technique  permits  the 
analyst  to  obtain  the  response  (at  the  sampling  in¬ 
stants)  of  any  system  state  from  Its  closed-loop 
transfer  function  expressed  in  the  complex  z-domain. 
If  it  is  desired  to  know  the  response  between 
sampling  instants,  either  the  submul ti pi e  method  or 
the  modified  z-transform  method  may  be  adapted  to 
the  cross-multiplication  technique. 

A  third,  more  sophisticated  technique  is  the 
Parameter  Space  Method.  It  is  being  developed  for 
determining  the  stability  and  dynamic  characteris¬ 
tics  of  a  digital  control  system  in  terms  of  sever¬ 
al  selected  system  parameters The  method  re¬ 
quires  that  the  system  characteristic  eeuation  be 
available  in  the  complex  z-domain.  Although  not 
necessary,  its  application  is  facilitated  by  aug¬ 
menting  the  analytical  results  with  graphical  por¬ 
trayals  in  a  selected  multi  parameter  space.  The 


method  is  based  on  analysis  and  synthesis  methods 
for  linear  and  nonlinear  control  system  design 
which  are  amply  described  in  Siljak's  excellent 
monograph  on  the  subject  essence,  the  para¬ 

meter  space  method  permits  the  designer  to  evaluate 
graphically  the  effects  of  the  locations  of  the 
roots  of  the  characteristic  equation.  Hence,  he 
may  design  the  control  system  in  terms  of  his  se¬ 
lected  performance  criteria;  e.g.,  absolute  stabil¬ 
ity,  damping  ratio,  settling  time.  He  is  able  to 
see  the  effects  on  the  characteristic  equation 
roots  (and  hence  on  system  dynamics)  of  changing 
several  adjustable  parameters.  The  method  has  been 
extended  to  portray  the  effect  of  varying  the 
sampling  period,  thereby  permitting  one  to  observe 
the  effect  of  the  choice  of  various  values  assigned 
to  the  sampling  period  on  absolute  and  relative 
stability.  Also,  simple  recursive  formulas  have 
been  derived  so  that  the  resulting  formulation  is 
deliberately  cast  in  a  form  particularly  amenable 
to  solution  by  a  digital  computer  or  a  desk  cal¬ 
culator,  emphasizing  the  interplay  between  analysis 
and  computing  machines. 

J.  Documentation 

A  major  portion  of  the  progress  reported  above 
has  been  documented  in  a  comprehensive  US  Army 
Missile  Corrmand  report. 13 

V,  Future  Plans 

Plans  for  the  next  fiscal  year  revolve  primar¬ 
ily  about  impl ementating  Task  3  (define  candidate 
G&C  systems)  of  the  development  program  as  well  as 
completing  any  partially-completed  portions  of  Task 
2  now  underway.  Close  communication  will  be  main¬ 
tained  with  the  intelligence  community  to  become 
aware  of  any  changes  to  the  presently  predicted 
targets  as  contrasted  to  the  attribution  mentioned 
above.  Oetailed  analytical  models  of  the  missile 
plant(s),  effectors,  and  autopilot(s)  will  be  devel¬ 
oped.  The  possibility  of  modifying  plant  charac¬ 
teristics  by  making  innovative  use  of  aerodynamics 
is  to  be  investigated.  Projected  advances  in  the 
field  of  propulsion  will  also  be  investigated.  The 
comparison  of  various  digitally-implemented  guid¬ 
ance  laws  will  be  continued,  including  those  that 
incorporate  DAC.  Applications  (other  than  to  guid¬ 
ance  laws)  of  DAC  theory  to  improve  system  perform¬ 
ance  will  be  investigated.  Development  of  a  mod¬ 
ular  guidance  simulation  to  implement  this  investi¬ 
gation  has  already  begun.  Trade  studies  concerning 
identified  existing  and  future  sensors  will  be  con¬ 
ducted,  with  the  possibility  of  developing  a  new 
sensor  with  capabilities  not  yet  in  existence.  Ad¬ 
vances  in  microcomputer  state-of-the-art  will  be 
watched  closely.  Finally,  tools  to  aid  the  G&C 
system  designer  to  handle  digital  impl  orientation 
will  continue  to  be  both  assessed  and  developed. 

VI.  Conclusions 

As  indicated  within  this  paper,  there  is  a  cley 
need  for  the  development  of  an  advanced  G&C  system 
for  US  Army  future  tactical  missiles.  This  need  is 
dictated  by  the  predicted  targets  of  the  future, 
the  anticipated  battlefield  of  the  future,  and  the 
characteristics  of  tactical  missiles  that  are  ei¬ 
ther  in  the  inventory  or  under  development  at  this 
time.  The  nature  of  the  future  G&C  system  will  be 
digital  so  that  the  overall  missile  system  may  be 
availed  of  present  and  predicted  advantages  that 
are  implicit  in  digital  controllers.  The  form  of 
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the  guidance  law  probably  will  be  optimal,  since 
the  performance  criteria  that  must  be  minimized 
must  take  more  Into  account  than  miss  distance. 
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ABSTRACT 

The  engagement  of  maneuvering  land  vehicles  with  gun  systems  place  extreme  performance  re¬ 
quirements  on  the  fire  control  system  designs.  The  effectiveness  of  a  gun  fire  control  sys¬ 
tem  depends  on  the  capability  to  provide  an  accurate  fire  control  solution,  i.e.  predict  the 
future  position  of  the  target  a  projectile  time-of-f light  later.  Non  linear  prediction  is 
shown  to  not  only  improve  performance  but  to  also  increase  available  time  for  firing  against 
maneuvering  targets.  Sub-optimal,  multi-variable,  adaptive  estimation  approaches  are  shown 
to  improve  the  effectiveness  of  predictive  fire  control  systems. 

Sensitivity  analyses  are  presented  that  relate  system  induced  errors  and  target  motion  in¬ 
duced  errors  to  tracking  noise  and  predictor  order.  Relationships  between  system  stability 
and  performance  for  two  basic  types  of  fire  control  systems  are  presented. 

INTRODUCTION 

This  paper  discusses  the  fire  control  system  problem,  the  nature  of  land  vehicle  mobility  and 
agility  and  the  ability  of  predictive  fire  control  systems  to  effectively  engage  maneuvering 
vehicles.  Existing  performance  specification  do  not  satisfactorily  describe  the  level  of 
maneuverability  expected  in  a  tactical  situation.  Rather,  present  specifications  define  per¬ 
formance  requirements  for  fixed  vehicle  speed  and  heading  movement  which  has  resulted  in  the 
development  of  fire  control  system  designs  that  are  significantly  degraded  in  a  maneuvering 
target  environment.  The  problem  is  addressed,  in  general,  for  the  four  cases  of  firing 
vehicle- target  vehicle  movements.  The  processes  required  in  the  fire  control  solution  are 
identified  and  the  sensitivity  of  system  performance  to  the  propagation  of  tracking  errors 
is  discussed.  The  stability  and  performance  characteristics  of  two  generic  fire  control 
system  configurations  are  analyzed  in  some  detail. 

GUN  FIRE  CONTROL  SYSTEM  PROELEM 

The  purpose  of  gun  fire  control  systems  is  to  have  a  projectile,  that  has  been  fired  a  time 
of  flight  previously,  impact  the  target  that  was  sighted  a  time  of  flight  earlier.  The 
critical  motion  parameters  that  degrade  the  performance  of  predictive  fire  control  systems 
have  been  identified  as  cyclic  oscillations  exhibiting  frequencies  that  are  within  the 
motion  capabilities  of  tactical  land  vehicles  (1).  Tracking  error,  defined  as  the 
difference  between  target  and  reticle  position,  does  not  in  itself  cause  the  performance 
degradation.  The  inability  of  the  fire  control  system  to  determine  the  motion  deriv¬ 
atives  of  the  line-of-sight  (LOS)  to  the  target,  and  predict  the  future  position  of  the 
target  are  the  two  main  factors  that  cause  fire  control  systems  degradation. 

The  error  in  the  ability  of  a  fire  control  system  to  cause  the  projectile  to  intercept  the 
target  a  time  of  flight  later  is  referred  to  as  total  gun  pointing  (TGP)  error.  TGP 
error  is  defined  as  the  offset  between  the  actual  gun  pointing  direction  at  round  exit  and 
the  location  of  the  target  centroid  at  round  impact.  The  TGP  error  is  the  sum  of  the  propa¬ 
gated  system  induced  (SI)  errors  and  target  induced  (TI)  errors  (i.e.  TGP  error  »  SI  errors 
+  TI  errors).  The  SI  errors,  considered  in  this  study,  are  the  tracking  error  (difference 
between  the  tracker  LOS  and  true  LOS  to  the  target)  at  the  time  of  firing  (lay  error)  and 
the  estimation  errors  (difference  between  estimated  LOS  states  and  true  LOS  states).  The 
SI  errors  propagated  through  a  projectile  time  of  flight  result  in  a  kinematic  lead  error. 

The  TI  error  is  caused  by  the  target  motion  during  the  time  of  flight  of  the  projectile. 

It  is  dependent  on  the  order  of  the  lead  solution  in  the  fire  control  system.  For  a  first 
order  lead  system  the  TI  error  is  the  difference  between  the  actual  LOS  movement  during  a 
projectile  time  of  flight  and  the  propagated  LOS  movement  assuming  perfect  LOS  rate  at  the 
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tima  of  fire.  The  first  order  predictor  system  TX  error  ignores  the  presence  of  actual  target 
acceleration  at  time  of  firing  and  during  projectile  flight  time.  For  a  second  order  lead 
system  the  TX  error  is  the  difference  between  the  actual  LOS  motion  during  a  projectile  tine 
of  flight  and  the  propagated  LOS  movement  assuming  perfect  LOS  rate  and  acceleration  at  the 
time  of  firing.  The  second  order  predictor  system  TI  error  accounts  for  target  acceleration 
at  time  of  firing  but  ignores  the  actual  target  acceleration  during  projectile  flight  time. 
This  distribution  of  errors  is  shown  in  Figure  1.  The  ballsitic  flight  characteristics  of 
the  projectile  are  ignored. 

The  fire  control  solution  occurs  during  a  short  time  interval  which  is  related  to  the  time 
of  flight  of  the  projectile.  The  motion  conditions  of  both  the  firer  and  the  target  are 
needed  tq  understand  and  solve  the  fire  control  system  problem.  Four  motion  conditions 
exist:  stationary  firer-stationary  target,  stationary  firer-moving  target,  n»ving  firer- 
stationary  target,  and  moving  firer-moving  target.  The  stationary  firer-stationary  target 
is  the  least  dynamic  situation  and  is  the  least  complex  case,  and  the  moving  firer-moving 
target  is  the  most  complex  case.  For  each  of  the  cases,  the  LOS  between  the  firer  and 
target  is  the  key  to  which  of  the  four  fire  control  processes  are  being  called  upon  in  a 
demanding  manner. 

FIRE  CONTROL  SYSTEM  PROCESSES 


A  fire  control  system  may  be  broken  down  into  four  distinct  processes.  Each  of  these  pro¬ 
cesses  are  present  in  all  types  of  fire  control  systems.  They  are:  tracking,  estimation, 
prediction,  and  gun  pointing.  In  specific  designs  these  four  processes  are  accomplished  in 
different  manners. 

The  tracking  process  is  Important  in  all  four  cases.  For  the  moving  firer  cases,  tracking 
becomes  more  critical  because  the  base  motion  of  the  firer  must  be  comoeruated  and  it  nav 
be  affected  in  a  secondary  manner  by  target  motion.  Tracking  is  usually  accomplished  man¬ 
ually  and  is  concerned  with  the  alignment  of  the  sight  reticle  with  the  target.  The  gunner 
is  involved  directly  at  this  stage  and  accuracy  of  tracking  will  be  a  characterization  of 
the  ability  of  any  given  gunner  to  perform  the  task.  Test  data  obtained  from  experimental 
investigations  can  be  used  to  determine  tracking  error  means,  standard  deviations,  and 
correlation  time  constants  useful  for  building  models  of  the  tracking  errors. 

The  estimation  process  is  the  intermediate  stage  between  the  tracking  process  and  the  pre¬ 
diction  process  and  its  configuration  is  dependent  upon  the  order  of  the  prediction  process. 
Estimation  is  the  process  of  filtering  the  tracking  data  to  provide  the  necessary  target 
motion  information  required  in  the  prediction  process.  The  accuracy  of  the  tracking  data 
will  influence  the  performance  of  the  estimation  process.  The  system  error  induced  by  Che 
estimation  process  decreases  with  improvement  in  tracking  accuracy. 

Prediction  of  target  future  position  to  obtain  intercept  between  projectile  and  target  is 
dependent  upon  an  estimate  of  the  present  motion  of  the  target  and  time  of  flight  of  the 
projectile.  The  output  of  the  estimator  is  not  a  complete  description  of  the  present  motion 
of  the  threat,  therefore,  the  predictor  does  not  have  the  necessary  information  to  calculate 
the  threat's  future  position  exactly.  If  restrictions  are  placed  on  the  allowable  threat 
motions,  then  the  predictor's  ability  to  determine  its  future  position  is  improved.  Over¬ 
simplification  of  allowable  threat  motions  has  placed  unrealistically  simplified  requirements 
on  the  operation  of  the  estimation  and  prediction  processes.  Realistic  threat  motions  are 
determined  by  the  mobility  capabilities  of  tactical  vehicles.  In  Che  past,  the  majority  of 
threats  that  have  been  studied  have  been  nonaccelerating.  The  requirements  of  an  estimator 
and  a  predictor  for  this  type  of  motion  are  to  combine  the  apparent  threat  velocity  estimate 
and  projectile  time  of  flight  for  the  lead  solution.  The  required  lead  is  constant  and  can 
be  realized  after  some  settling  time.  The  existence  of  acceierating  targets  requires  the 
estimator  and  predictor  to  develop  constantly  changing  lead  angles,  hence,  the  need  for  non¬ 
linear  prediction. 

An  important  point  to  observe  is  that  for  the  stationary  firer-moving  target  case,  the 
prediction  process  is  required  to  provide  gun  conanand  orders  that  orient  the  gun  to  account 
for  target  motion  during  the  projectile's  time  of  flight,  whereas  In  the  moving  firer- 
stationary  target  case  this  prediction  process  is  not  required  because  the  LOS  existing 
between  the  firer  and  target  at  instant  of  firing  does  not  move  during  the  projectile's 
time  of  flight.  For  the  moving  firer-moving  target,  LOS  also  moves  after  projectile  firing. 

The  gun  pointing  process  is  required  to  align  and  stabilize  the  gun  along  the  predicted  LOS 
to  the  target.  The  stabilization  and  response  of  the  gun  pointing  loop  is  a  major  concern 
for  fire  control  system  performance  against  maneuvering  targets.  Stabilization  of  the  gun 
pointing  process  could  have  an  adverse  effect  on  overall  system  performance.  The  moving 
firer  cases  will  stress  the  gun  pointing  process  most  severely  but  it  is  possible  that  the 
gun  pointing  process  will  be  equally  stressed  for  the  stationary  firer-moving  target  case  with 
non  linear  prediction. 
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FIRE  CONTROL  SYSTEM  CONFIGURATIONS 

The  three  currently  used  fire  control  configurations  are  known  as  manual  or  iron  sight,  dis¬ 
turbed  reticle  and  stabilized  sight-director  systems.  A  fourth  mehtod  called  closed  loop 
refers  to  projectile  spotting  to  adjust  the  fire  control  solution  and  is  not  considered  in 
this  discussion.  The  manual  fire  control  system  uses  the  brute  force  approach  and  concentrates 
on  stabilizing  the  gun  position  exclusively.  In  this  system  the  lead  is  introduced  manually, 
therefore,  there  is  no  automation  of  the  fire  control  estimation  and  prediction  process.  The 
disturbed  reticle  system  stabilizes  the  gun  position  and  disturbs  the  position  of  the  track¬ 
ing  reticle  from  the  gun  line  position.  In  this  scheme  the  tracking,  estimation  and  predic¬ 
tion  processes  are  inseparable  and  the  fire  control  solution  is  automated.  The  rejection  of 
firer  vehicle  base  motion  is  difficult  to  accomplish  in  this  type  of  system.  The  last  system 
to  be  considered  is  the  stabilized  sight-director  system.  The  tracking  process  is  accomplish¬ 
ed  by  a  tracker  which  is  rotated  from  a  stabilized  base  that  has  signals  applied  to  isolate 
firer  vehicle  base  movement.  The  resulting  LOS  orientation  is  referenced  to  interial  space, 
as  contrasted  to  the  gun  line  for  disturbed  reticle  systems.  The  estimation  process  is  the 
intermediate  link  between  the  tracking  process  and  the  prediction  process.  The  prediction 
process  uses  the  estimation  process  outputs  combined  with  projectile  time  of  flight  to  deter¬ 
mine  the  gun  pointing  commands.  The  gun  pointing  process  uses  the  estimated  LOS  to  the  target 
summed  with  the  calculated  lead  offset  of  the  gun  from  the  tracker  LOS  to  position  the  gun 
line. 


How  well  a  fire  control  system  configuration  performs  is  a  function  of  target  movement, 
firer  movement  and  fire  control  system  design.  The  analytical  methodology  required  to  study 
this  problem  should  be  constrained  to  real  time  solution  mechanisms.  Another  way  to  say  this 
is:  post  data  analysis  techniques  using  data  obtained  from  field  tests  will  not  provide  the 
insight  that  is  required  to  obtain  an  understanding  of  the  relative  performance  of  different 
fire  control  systems.  Probability  of  hit  information  is  useful  for  an  assessment  of  systems 
that  have  been  fielded  but  is  not  applicable  for  trade  off  studies  of  the  type  required  in 
this  study.  Analytical  methodologies  such  as  servo  mechanism  synthesis  and  modem  filtering 
technology  are  required  to  study  this  problem. 

MANEUVERING  TARGET  DESCRIPTION 

A  quantitative  description  -of  the  threat  is  required  to  evaluate  the  performance  of  fire 
control  systems  operating  against  maneuvering  targets.  To  develop  this  description,  it  is 
necessary  to  consider  the  mobility  and  agility  characteristics  of  threat  vehicles  in  a  real¬ 
istic  combat  environment.  A  thorough  description  of  anticipated  maneuvering  seems  to  defy 
identification  because  threat  maneuvers  constitute  a  large  set  of  possibilities  even  when 
constrained  by  tactical  doctrine,  driver  policy,  terrain  and  vehicle  capabilities.  Two 
approaches,  analytical  and  empirical,  are  available  for  consideration  in  the  attempt  to 
identify  the  maneuver  characteristics  of  land  vehicles.  An  analytic  approach  would  view  each 
maneuver  as  being  composed  of  elements  from  an  idealized  group  of  movements.  An  empirical 
approach  would  view  the  maneuvers  as  having  actually  occurred  during  limited  tests  for  diff¬ 
erent  types  of  maneuvering  vehicles.  Neither  of  these  approaches  provide  a  complete  maneuver 
description,  but  a  combination  of  these  two  approaches  offers  some  advantages  and  is  the 
rationale  adopted.  The  analytic  approach  will  partially  overcome  the  incompleteness  of  the 
empirical  data  base  while  the  empirical  data  will  offset  the  mathematical  idealizations  of 
the  analytic  methodology. 

Empirical  Approach 

When  U3ing  empirical  data  to  demonstrate  the  performance  of  a  gun  fire  control'  system,  base¬ 
line  performance  can  be  determined  with  no  concerns  arising  from  idealization  of  the  man¬ 
euvers.  Since  the  number  of  maneuvers  will  be  rather  small,  they  neither  provide  sufficient 
information  about  the  robustness  of  a  fire  control  design  mehtodology  nor  the  pathology  when 
the  fire  control  system  begins  to  degrade.  When  demonstrating  the  performance  of  a  fire 
control  system  against  experimental  data,  caution  must  be  exercised  to  assure  that  the 
empirical  data  is  properly  inputted  to  the  fire  control  system  model.  Matching  of  the  data 
rates  and  noise  levels  often  requires  some  preprocessing  of  experimental  data  to  prepare  it 
for  use  in  simulation  studies. 

Analytic  Approach 

As  a  supplement  to  the  empirical  approach  the  analytic  approach  is  used  to  investigate  sen¬ 
sitivity  effects  for  a  larger  group  of  movements.  Simulating  new  or  pathological  maneuvers 
require  that  the  analytic  capability  superimpose  maneuvers  arising  from  random  disturbances 
and  intentional,  voluntary  vehicle  driver  commands .  The  random  disturbances  may  be  represented 
in  terms  of  time  histories  or  power  spectral  densities.  The  time  history  approach  is  based 
on  the  development  of  a  mathematical  model  of  vehicle  movement  influenced  by  terrain  affects 
and  arbitrary  driving  habits  of  individual  drivers.  It  is  assumed  that  for  no  random  effects 
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caused  by  terrain  irregularities  or  driver  input,  the  vehicle  would  follow  a  straight  line- 
conatant  speed  path.  Maneuvers  are  viewed  as  perturbations  on  this  straight  line-constant 
speed  path.  Apparent  acceleration,  a  (t),  which  is  correlated  in  time,  accounts  for  the 
vehicle's  deviation  from  a  straight  line  path.  Maneuver  capability  is  expressed  by  three 
quantities:  the  variance,  or  magnitude  of  a  (t),  the  cyclic  maneuver  frequency  and  the  time 
constant  of  the  maneuver. 

Intentional,  voluntary  vehicle  driver  cooaaanded  motion  of  land  vehicles  over  terrain  is  a 
complicated  subject  in  itself  and  will  not  fee  investigated  in  this  study.  It  is  recognized 
however,  that  an  interaction  between  vehicle  horsepower,  weight,  suspension,  and  locomotion 
concepts  do  combine  with  terrain  over  which  it  is  moving  to  provide  different  levels  of 
mobility  with  respect  to  a  fixed  reference  frame.  Therefore,  different  vehicle  designs  will 
have  different  mobility  levels  defined  in  terms  of  motion  and  derivatives  of  motion.  Agil¬ 
ity  is  closely  related  to  mobility  and  yet  it  is  a  slightly  different  description  of  inten¬ 
tional  vehicle  motion.  Where  mobility  describes  the  movement  of  a  vehicle  from  one  location 
to  another  location  in  a  given  period  of  time,  agility  describes  the  vehicle's  ability  to 
alter  its  mean  path. 

SENSITIVITY  OF  FIRE  CONTROL  PROCESSES 

Degradation  in  gun  pointing  accuracy  results  from  two  major  error  sources,  system  and 
target  induced  errors.  The  target  induced  errors  are  caused  by  the  motion  of  the  target 
during  the  time-of-flight  of  the  projectile*  Since  the  target  has  the  capability  to  man¬ 
euver  within  constraints  of  the  terrain,  vehicle  characteristics  and  driver  policy  during  a 
projectile's  time-of-flight,  there  is  no  such  thing  as  a  correct  (perfect)  lead  solution. 

The  lead  solution  is  based  on  the  projected  target  position  using  the  present  target  states 
and  projectile  time-of-flight.  Therefore,  the  target  induced  error,  in  general,  cannot  be 
reduced  to  2ero  for  a  maneuvering  target.  However,  it  can  easily  be  shown  that  the  pre¬ 
diction  process  is  capable  of  reducing  the  gun  pointing  error  due  to  target  motion. 

The  system  induced  errors  are  made  up  of  bias  and  random  errors  emanating  from  specific 
components  and  subsystems.  The  propagation  of  these  errors  degrade  the  performance  of  the 
fire  control  system.  The  system  induced  errors  of  major  concern  are  those  occurring  in  the 
tracking  process.  Sensitivity  analyses  have  bean  performed  to  evaluate  the  degradation  of 
gun  pointing  commands  to  tracking  process  errors. 

The  analysis  considers  the  fire  control  system  processes  to  be  interfaced  in  tandem  with 
no  feedback  of  outputs  to  a  previous  process.  The  analysis  is  further  limited  to  a  segment 
of  a  maneuvering  target  path  which  was  generated  by  a  maneuvering  target  path  simulation 
program.  This  analytically  generated  path  provides  an  exact  time  history  of  the  target  states 
(position,  velocity  and  acceleration). 

The  tracking  process  is  modelled  by  summing  random  errors  of  known  variance  with  the  output 
of  a  perfect  LOS  sensor.  The  output  of  the  tracking  process  is,  by  definition,  in  LOS 
coordinates,  however,  cartesian  coordinates  are  used,  by  choice  and  not  a  limitation  of  the 
methodology,  in  the  estimation  processes.  For  simplication,  the  transformation  from  LOS  to 
cartesian  coordinates  is  accomplished  prior  to  adding  tracking  noise. 

A  sub-optimal,  adaptive  Kalman  filter  (KF)  is  used  for  the  estimation  process  in  the  generic 
fire  control  system  under  consideration.  The  noisy  tracking  process  signal  is  processed  by 
the  KF  to  provide  a  "best"  estimate  of  the  target  states  (position,  velocity  and  acceleration). 
The  estimation  errors  are  minimized  by  providing  the  filter  with  the  correct  variance  of  the 
observation  noise.  In  practice,  this  perfect  match  of  noise  variance  is  not  achieveable  but 
can  be  approached  with  detailed  error  analysis  of  the  tracking  process  or  with  software 
methodology  to  estimate  the  noise.  The  latter  is  probably  desirable  and  necessary  because 
the  variance  of  the  tracking  process  error  is  not  time  invariant  in  a  combat  environment. 

The  KF  is  the  generic  fire  control  system.  The  KF  equations  and  theory  are  well  known  and 
are  presented  elsewhere  (2,3).  However,  the  adaptive  feature  of  the  designed  KF,  which 
requires  on-line  computation  of  the  filter's  gain,  is  outlined  (4).  The  adaptive,  time 
varying  gain  is  obtained  by  changing  the  variance  of  the  uncertainty  of  the  embedded  target 
dynamics,  as  a  function  of  the  estimated  path  geometry.  The  forcing  function  for  the  target 
dynamics  is  modelled  as  a  random  (Gaussian  noise)  rate  of  change  of  acceleration.  The 
variance  of  y  is  defined  in  the  body  coordinates  of  the  target  as  constant,  diagonal  elements 
of  the  Q  matrix.  The  Q  matrix  is  rotated  as  the  target  maneuvers  co  provide  a  time  varying 
Q  matrix  in  the  filter’s  coordinate  system. 

The  sensitivity  of  the  estimates  to  the  tracking  process  noise  is  evaluated  for  a  typical 
maneuvering  target  path.  The  ground  track  of  the  maneuver  i^  shown  in  Figure  2.  The  max¬ 
imum  speed  and  lateral  acceleration  are  10  m/sec  and  2  m/sec~,  respectively.  Figure  3  shows 
the  degradation  in  velocity  estimates  as  the  standard  deviation  of  the  tracking  process  noise 
on  the  assumed  position  observation  is  increased  from  0.05  meter  to  1.0  meter.  The  degrada- 
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tion  in  the  estimates  of  lateral  acceleration  for  the  same  noise  levels  is  shown  in  Figure  4. 
A  comparison  of  these  two  figures  shows  that  the  velocity  estimates  are  not  as  sensitive  to 
the  propagation  of  tracking  noise  as  the  acceleration  estimates. 

The  prediction  process  provides  the  command  for  pointing  the  gun  to  the  predicted  target 
position.  The  estimated  future  position  of  the  target  depends  on  the  order  of  the  prediction 
process.^  Ideally,  one  would  like  to  forecast  the  target  position  so  that  a  projectile  fired 
a  time-of- flight  earlier  would  arrive  at  a  point  in  space  simultaneously  with  the  target. 
Unfortunately,  only  the  present  states,  which  are  never  known  exactly,  are  available  for  use 
in  computing  future  target  position. 

With  knowledge  of  the  true  future  position  of  the  target,  available  from  the  target  motion 
model,  the  degradation  in  the  gun  pointing  commands  can  be  evaluated  for  different  tracking 
errors.  Target  induced  errors  and  the  propagation  of  the  tracking  process  noise  are  analyzed 
to  evaluate  their  effect  on  gun  pointing  commands. 

The  target  induced  errors  are  functions  of  target  maneuver  characteristics,  projectile 
time  of  flight  and  prediction  order.  For  a  given  prediction  order  and  with  perfect  know¬ 
ledge  of  the  present  target  states  and  time  of  flight,  the  resulting  target  induced  errors 
are  lower  bound  prediction  errors.  Effects  of  time  of  flight  and  order  of  prediction  are 
shown  in  Figure  5  for  a  maneuvering  target  whose  maximum  speed  and  lateral  acceleration  is 
2 

10  m/sec  and  3.5  m/sec  .  Prediction  errors  are  improved  for  decreases  in  time-of- f light  and 
higher  order  of  prediction. 


First  order  prediction  is  linear  and  requires  only  accurate  estimates  of  velocity  to  approach 
the  lower  bounds  of  prediction  error.  Second  order  prediction  requires  not  only  accurate 
velocity  but  also  acceleration  estimates  to  minimize,  the  prediction  errors.  Figure  6  shows 
the  standard  deviation  of  prediction  error  for  the  target  maneuver  shown  in  Fig,.  2  as  a  func¬ 
tion  of  time-of-f light  and  variances  of  tracking  process  noise  for  first  order  prediction. 
These  results  indicate  that  the  degradation  in  prediction  error  is  minimized  as  the  quality 
or  tracking  improves.  However,  the  existence  of  the  lower  bound  curve  for  second  order  pre¬ 
diction  provides  additional  improvement,  not  realized  by  first  order  prediction.  Assuming 
position  observations  (input  to  the  KF)  with  a  lo  noise  of  1.0  meter,  Figures  6  and  7  show 
that  there  is  no  large  difference  between  first  and. second  order  prediction.  However,  second 
order  prediction  with  a  reduction  in  the  tracking  process  error  to  0.05  meter  (25  microradians 
at  a  range  of  2000  meters)  provides  a  significant  improvement  in  the  lead  solution.  Unlike 
first  order  prediction,  second  order  prediction  is  not  only  more  sensitive  to  the  tracking 
process  noise  but  also  to  the  observation  state.  Figure  7  shows  that  improvements  are  re- 
t/  ,  0h3*™az,±0TLS  are  rates  rather  than  position.  If  tracking  accuracies  of  0.04 

m/aec  (20  microradian/ sec  at  2000  meters)  are  achieved,  the  prediction  error  is  within  about 
ten  percent  of  the  lower  bound  for  second  order  prediction. 

lead  errors  discussed  above  are  the  differences  between  the  predicted  and  actual  target 
?n!fCt0n3/0r^anreStiiaat:6d  ti^e-of- flight .  Targets  are  not  point  sources  and  a  more  nean- 
ingiul  criteria  for  evaluating  the  system  processes  is  the  percent  time  on  target  for  a 

rn^*8lan?a8!ment  Assuming  a  tarSeC  size  of  2.3  meters  X  2.3  meters,  independent  of 

^  5en^tl0n’  C^e  ?ercent  t1®*5  target  for  the  same  tracking  accuracy  in  Figure  7  is 
depicted  in  Figure  8  for  times  of  flight  between  1.0  seconds  and  2.5  seconds, 
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f^re  contro1  configurations  in  existence:  manual,  disturbed  reticle  and 

casses^are  iden£lfied  of  how  thi  fire  control  pro!  ““ 

th!  nech“±zed'  ^  agisting  operational  systems  utilize  the  human  operator  to  null 

T  the  °bserved  tar8et  and  the  «ticle  position.  !£Td!gre!“f  p!rc£i- 

and  thatch?’  C  6  track1^8,  estimation  and  prediction  processes,  are  performed  by  the  man 

reticle  is  dlst^bed’i  ^  ?0sitlon  serve  as  the  reference  from  which  the 

the  distrubed  in  the  disturbed  reticle  system.  Involvement  of  the  human  gunner  in 

the  .urret  loop  for  the  distrubed  reticle  system  and  his  absence  from  the  turret  loop  for  the 
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stabilized  sight-director  system  is  a  distinguishing  feature  of  the  systems.  The  tracking 
process  is  therefore  more  isolated  from  the  estimation,  prediction,  and  gun  pointing  processes 
in  the  stabilized  sight-director  system. 

Disturbed  Reticle  Fire  Control  System 

One  fire  control  configuration  in  current  use  is  the  disturbed  reticle  concept.  The 
following  discussion  is  intended  to  describe  in  detail  the  functioning  of  the  disturbed 
reticle  fire  control  system  and  identify  the  four  processes,  showing  how  each  is  related  to 
the  other.  Figure  9  describes  the  signal  flow  in  such  a  system  and  the  four  major  processes 
have  been  identified  in  terms  of  where  in  the  system  each  is  accomplished. 

The  input  to  the  system  is  the  LOS  from  the  target  to  the  reticle  of  the  tracking  system. 

The  human  operator,  present  in  most  current  systems,  moves  the  handle  bar  controller  to 
align  the  reticle  of  the  tracking  system  to  be  coincident  with  the  target.  The  ability  of 
any  human  controller  to  accomplish  this  task  defines  the  quality  of  the  tracking  process. 
Handle  bar  controller  output,  which  is  directly  related  to  the  LOS  rate,  is  used  to  drive  two 
independent  subsystems.  The  first  is  the  turret  servo  which  is  commanded  to  rotate  at  a  rate 
directly  proportional  to  the  handle  bar  controller  deflection.  The  second  subsystem  driven 
by  the  handle  bar  controller  is  a  lead  screw  servo  and  reticle  system.  The  displacement  of 
the  lead  screw  servo  is  directly  proportional  to  the  filtered  handle  bar  controller  deflection 
multiplied  by  the  projectile  time  of  flight.  The  lead  screw  displacement  is  used  to  position 
the  reticle  of  the  tracking  system. 

These  are  two  distinct  feedback  signal  paths  in  the  disturbed  reticle  configuration  and  the 
human  is  a  series  subsystem  in  both  paths.  Another  important  observation  is  to  note  that 
the  signal  loop  made  by  the  turret  servo-man-handle  bar  controller  is  a  degenerative  feedback 
loop  because  of  the  negative  summing  junction.  The  signal  loop  made  by  the  filter-time  of 
flight  lead  servo-reticle  servo-man  handle  bar  controller  is  a  regenerative  feedback  loop 
because  or  two  negative  summing  junctions.  During  normal  operation  of  the  disturbed  reticle 
system,  the  performance  of  these  two  feedback  paths  give  rise  to  a  dynamical  system  that 
exhibits  some  undesirable  performance  characteristics.  Without  further  crossfeed  compensa¬ 
tion,  the  closed  loop  performance  of  the  disturbed  reticle  system  is  at  best  marginaliy  stable 
and  at  worst  unstable.  To  overcome  this  condition,  compensation  signal  paths  are  added.  The 
basic  compensation  is  a  tachometer  generator  signal  from  the  lead  screw  servo  which  is  com¬ 
bined  with  the  turret  servo  error  signal.  This  composite  signal  is  fed  to  the  turret  servo 
and  the  reticle  servo  to  compensate  for  the  dynamical  mismatch  that  occurs  in  the  reticle  and 
turret  servos.  However,  there  is  no  such  thing  as  a  perfect  compensation  and  the  undesirable 
performance  characteristics  alluded  to  earlier  can  never  be  completely  nullified,  not  to  men¬ 
tion  the  potentially  precarious  situation  that  might  occur  if  any  failure  or  shift  occurs  in 
the  compensation  paths. 

The  important  thing  to  observe  about  the  root  locations  in  the  figures  10  and  11  is  that 
there  are  numerator  roots  in  the  right  half  of  the  S  plant.  This  arises  from  the  basic 
disturbed  reticle  configuration  and  must  be  considered  a  fixed  element  phenomenon  in  this 
type  of  system.  The  poles  or  denominator  roots  describe  the  system  operating  point  for  a 
system  gain  of  zero.  The  zeros  or  numerator  roots  describe  the  system  operating  point  for  a 
system  gain  of  infinity.  The  dotted  trajectories  connecting  these  two  extremes  are  a  pictor¬ 
ial  description  of  the  operating  point  loci  for  all  intermediate  gains.  These  systems' exhibit 
conditional  stability  because  of  the  presence  of  positive  feedback  in  the  equivalent  transfer 
function  between  B  and  A.  These  are  different  closures  than  exist  for  a  negative  feedback 
that  occurs  when  both  the  reticle  and  turret  crossfeeds  are  present  as  shown  in  Figure  12. 

The  existence  of  these  simultaneous  crossfeeds  from  the  lead  screw  servo  and  turret  servo 
error  to  the  turret  servo  and  reticle  servo  tend  to  offset  the  non-minimum  phase  root  con¬ 
dition  shown  in  Figures  10  and  11. 

In  summary,  it  is  the  location  of  the  operating  points  that  determine  the  system  stability 
characteristics.  The  frequency  content  of  the  tracking  error  is  directly  related  to  the 
operating  points,  but  equally  important  is  the  magnitude  of  the  tracking  error  which  is  in¬ 
fluenced  by  the  location  of  the  numerator  roots  of  the  closed  loop  transfer  function.  These 
effects  are  interrelated,  but  the  fundamental  underlying  requirement  is  to  achieve  an  adequate 
stability  margin  of  the  closed  loop  system.  This  stability  consideration  is  important  for 
fire  control  system  performance  and  the  designers  must  take  these  factors  into  account.  The 
end  result  is  system  performance  which  nay  be  acceptable  or  not  acceptable. 

It  can  be  asked  why  so  much  concern  about  this  situation  because  disturbed  reticle  systems 
have  performed  satisfactorily  in  the  past.  Perhaps  this  is  so,  but  with  the  introduction 
of  maneuvering  targets,  the  performance  of  this  type  of  system  mav  be  adversely  affected. 

When  the  target  LOS,  0^,  shown  in  Figure  9  moves  at  a  constant  rate,  the  human  operator  is 
required  to  move  the  handlebar  controller  a  nominal  fixed  amount.  The  turret  servo  develops 
a  tixed  nominal  rate  and  the  lead  screw  servo  assumes  a  fixed  nominal  position.  It  then 
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becomes  the  task  of  the  human  to  perturbate  the  handle  bar  controller  about  this  nominal 
position  in  order  to  minimize  the  tracking  error.  When  the  target  LOS  rate  is  not  constant, 
which  is  the  situation  for  maneuvering  targets,  the  handle  bar  controller  must  be  moved 
consistent  with  the  changing  target  LOS  rate.  The  nominal  handle  bar  controller  position 
is  not  the  only  difference  in  the  system  operation  for  maneuvering  targets.  The  turret 
servo  accelerates  and  decelerates  and  the  lead  screw  servo  is  constantly  being  driven  to  a 
new  position.  The  position  of  the  reticle  is  a  result  of  these  two  signals  paths.  The 
dynamic  performance  mismatches  are  quaranteed  to  be  greater  than  for  the  non-maneuvering 
targets  and  the  tracking  performance  will  be  degraded.  This  degradation  occurs  from  the 
inability  of  the  closed  loop  system  to  accurately  null  the  constantly  changing  target  LOS 
rate.  The  extent  of  this  degradation  may  not  be  immediately  obvious  to  the  casual  analyst, 
but  the  oscillatory  nature  of  this  degradation  will  be  observed  once  a  sufficiently  close 
survey  of  the  tracking  error  is  made.  It  is  imperative  that  the  resulting  stability  margin 
of  the  closed  man-machine  system  be  large  to  insure  acceptable  performance  against  maneuv¬ 
ering  targets. 

Recent  work  has  shown  that  tactical  targets  can  execute  maneuvers  of  such  a  nature  that 
when  projectile  times  of  flight  of  1. 5-2.0  sec  are  considered,  target  induced  motion  after 
projectile  firing  will  cause  excessive  miss  distances  when  linear  predictor  fire  control 
systems  are  assumed  and  more  over  these  miss  distances  can  be  significantly  reduced  when  non¬ 
linear  or  higher  order  predictor  fire  control  systems  are  employed.  These  observations  indi¬ 
cate  lower  boundary  miss  distances  are  possible  for  non  linear  lead  systems.  When  this  sit¬ 
uation  is  presented  to  the  fire  control  designer,  his  inclination  will  be  to  consider  the 
possibility  of  including  nonlinear  prediction  in  the  fire  control  system.  In  the  disturbed 
reticle  configuration  shown  in  Figure  9,  this  may  be  a  design  impossibility  because  of  the 
level  of  tracking  performance  obtainable  from  the  operation  of  the  disturbed  reticle  systems. 
To  be  more  specific,  the  tracking  error  to  develop  usable  target  accelerations  for  nonlinear 
prediction  is  required  to  be  smaller  than  the  tracking  error  for  first  order  prediction. 

The  trade-off  between  the  propagated  system  induced  tracking  errors  for  the  nonlinear  esti¬ 
mation  process  must  be  offset  by  the  target  induced  prediction  error  improvements  realized 
by  the  higher  order  prediction.  The  key  ingredient  for  this  situation  to  exist  in  a  fire 
control  system  is  to  have  high  quality  tracking  errors. 

If  the  human  tracker  is  replaced  by  an  automatic  tracker  the  performance  limitations  imposed 
by  the  loop  structures  in  a  disturbed  reticle  system  may  negate  the  potential  improvement 
attainable  from  the  improved  tracking.  It  is  the  coupled  nature  of  the  tracking,  estimation, 
prediction,  and  stabilization  process  occurring  in  the  disturbed  reticle  configuration  that 
restrict  its  growth  to  better  fire  control  system  performance,  especially  against  maneuvering 
targets. 

Stabilized  Sight-Director  Fire  Control  System 

A  stabilized  sight-director  fire  control  system,  shown  in  Figure  13  is  actually  two  distinct 
systems  that  are  brought  together  to  accomplish  the  tracking,  estimation  and  prediction  pro¬ 
cesses  of  a  fire  control  system.  Stabilization  of  the  tracking  system  is  independent  from 
stabilization  of  the  turret.  The  stabilized  sight  is  decoupled  from  turret  and  hull  motion 
by  the  reverse  torquing  of  the  outer  gimbal  of  the  tracker  to  account # for  distrubances  of 
the  tracker  base  which  is  mounted  on  the  turret.  This  decoupling  enhances  the  ability  of  the 
tracker  to  maintain  coincidence  between  the  sight  reticle  and  the  target  LOS.  The  stabilized 
reticle  position  can  utilize  both  position  and  rate  feedback  to  augment  the  stability  of  the 
sight.  The  orientation  of  the  sight  reticle  is  therefore  an  independent  process  from  the 
turret  motion.  • 

Position  and  rate  of  the  LOS  are  fed  to  a  filter  or  estimation  process  to  determine  the 
necessary  information  about  the  LOS  to  the  target  that  will  be  needed  to  offset  the  turret 
servo  from  the  stabilized  tracker.  Multi-variable,  sub-optimal  technology  can  be  applied  to 
further  Improve  the  quality  of  tracking  that  can  be  realized  from  the  stabilized  sight- 
tracker.  Therefore  either  linear  or  non  linear  estimates  of  LOS  movement  can  be  considered 
as  possibilities.  If  LOS  accelerations  are  to  be  estimated,  the  appropriate  modeling  of 
target  dynamics  and  tracker  uncertainties  will  be  required  to  insure  that  the  degree  of 
suboptimality  is  not  excessive.  One  very  significant  plus  that  couples  the  estimation  and 
tracking  process  in  a  favorable  manner  is  the  utilization  of  3ight  line  rate  aiding  feedback 
to  the  tracker  obtained  from  estimation  of  the  target  rates  and  acceleration.  This  concept 
relaxes  the  task  of  the  human  tracker  or  auto- tracker  and  will  improve  the  minimization  of 
tracking  error. 

Output  of  the  target  state  estimator  is  used  in  two  separate  paths.  The  first  path  uses 

5t  and  9^  to  drive  the  turret  servo  as  a  director  to  fallow  the  tracker  LOS.  The  second 
signal  path  combines  target  state  estimates  with  projectile  time  of  flight  and  offsets 
the  gun  from  the  tracker  LOS  by  the  appropriate  value  to  permit  intercept  of  projectile  and 
target  a  time  of  flight  later. 
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Performance  of  the  stabilized  sight-director  system  should  not  be  compromised  by  maneuvering 
targets  to  the  extent  that  the  disturbed  reticle  system  is  compromised.  The  basic  reason 
for  this  is  that  the  cracking  system  is  essentially  decoupled  from  the  lead  prediction  sys¬ 
tem.  However,  there  are  some  inherent  stabilization  problems  that  can  occur  in  this  config¬ 
uration  and  they  are  accentuated  by  the  temptation  to  obtain  high  preformance  of  the  gun 
pointing  process.  The  argument  goes  as  follows:  with  increased  tracker  performance,  the 
gun  stabilization  servo  can  be  made  to  perform  more  rapidly,  thereby  increasing  the  overall 
capability  of  the  system.  However,  with  increased  performance  being  required  of  the  turret 
servo  to  follow  the  turret  commands,  the  stability  of  the  turret  servo  may  be  compromised 
because  of  the  high  gains  in  the  director-follower  loop.  Experience  with  similar  types  of 
systems  has  shown  that  because  of  non-rigid  gun  tube  and  hull  structures,  the  follower  loop 
system  must  be  phase  stabilized  and  not  gain  stabilized,  as  is  the  case  for  less  responsive 
systems  such  as  disturbed  reticle  systems.  This  requires  sophisticated  compensation 
circuits  to  overcome  system  instabilities. 

The  stabilized  sight  is  identified  between  the  target  input  and  the  sight  output  in  Figure  13 
It  in  turn  drives  the  gun  turret  servos  which  are  used  to  position  the  base  of  the  stabilized 
sight.  The  signal  flow  diagram  and  root  loci  for  the  stabilized  sight-director  system  is 
shown  in  Figure  14.  This  is  the  same  basic  root  locus  obtained  in  the  disturbed  reticle 
system  when  the  cross  feeds  were  included.  The  dotted  lines  show  the  movement  of  the 
stability  as  the  gain  is  increased.  The  addition  of  series  compensation  Circuits  in  the 

,  ,  ,  T,  S+l  ,  which  can  easily  be  added  in  a  straight-forward 

tracker  transfer  runction;  such  as  1  *  7 

T  S+l 

manner  will  alter  the  shape  of  the  loci  to  obtain  an  optimized  operating  point,  which  would 
be  difficult  in  the  disturbed  reticle  system.  The  fundamental  purpose  of  the  tracking  pro¬ 
cess  is  to  align  3g  with  3,.  Simultaneously  any  disturbances  on  the  stabilized  sight  are 
compensated  by  orientation ^of  the  sight  base  thereby  simplifying  the  tracking  task. 

CONCLUSIONS 

The  inherent  ability  of  a  stabilized  sight-director  fire  control  system  to  decouple  the  track 
ing  estimation,  prediction  and  gun  pointing  processes  may  be  exploited  to  improve  effective¬ 
ness  when  engaging  maneuvering  targets.  Accurate  tracking  is  necessary  for  non  linear 
prediction  and  multi- variable,  sub-optimal  design  technology  is  required  to  achieve  the 
needed  accuracy  of  the  target  state  estimates  for  mechanizing  nonlinear  prediction.  Further 
studies  are  required  to  identify  the  specific  details  of  the  resulting  system  design.  A 
complementary  methodology  employing  stability  and  performance  analyses  will  assist  in  this 
quest. 
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Figure  1.  Fire  Control  System  Error  Sources. 


Figure  2.  Typical  Maneuvering  Target  Path. 


Figure  4.  Sensitivity  of  Acceleration  Estimates  to  Observation  Noise. 
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Figurt  7  Second  Ord er  Prediction  Error*. 
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Figure  9  Disturbed  Reticle  Fire  Control  System. 
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Figure  10.  Disturbed  Reticle  Root  Locus, 
Fixed  Elements, 
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Figure  11.  Disturbed  Relicle  Roof  Locus, 
Turret  Servo  Crossfeed. 


Figure  12.  Disturbed  Reticle  Locus  with  Turret 
Servo  and  Reticle  Gossfeeds. 
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Figure  13.  Stabilized  Sight -Director  Fire  Control  System. 


Figure  14.  Stabilized  Sight  -  Director  Root  Locus. 
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DETECTION  AND  INITIATION  OF  FIRING  COMMAND  FOR  AN  ACCELERATION  PREDICTOR 
FIRE  CONTROL  SYSTEM  ENGAGING  MANEUVERING  TARGETS 

H.  H.  Burke  -  AMSAA 


Maneuvering  targets  degrade  the  performance  of  all  state  of  the  art  predic¬ 
tive  fire  control  systems.  Anti-aircraft  guns,  and  tank  guns  are  the  two  main 
weapon  systems  whose  effectiveness  is  lessened  by  the  presence  of  maneuvering 
targets.  The  application  of  the  methodology  described  in  this  paper  will 
significantly  increase  the  performance  for  these  systems,  especially  when  the 
jinking  or  agile  movement  of  the  target  increases. 

A  gun  fire  control  system’s  function  is  to  offset  the  gun  line  from  the 
target  line-of-sight  to  cause  a  projectile  to  intercept  the  target  a  time-of- 
flight  after  firing  the  gun.  Two  classifications  of  target  motion  occur;  non¬ 
maneuvering  and  maneuvering.  The  projectile-target  closest  approach  is  a 
measure  of  miss  distance  and  the  performance  of  the  gun  fire  control  system. 
Non-maneuvering  targets  require  a  constant  offset  between  the  target  line-of- 
sight  and  the  gun  line.  The  magnitude  of  the  required  offset  is  the  product  of 
the  target  velocity  perpendicular  to  the  line-of-sight  and  the  projectile  time- 
of-f light.  Maneuvering  targets  require  a  time  varying  offset  having  a  magnitude 
related  to  the  non-maneuvering  target  offset  plus  an  additional  amount  related 
to  the  target  acceleration,  acceleration  rate,  etc.,  perpendicular  to  the  line- 
of-sight  combined  with  the  appropriate  functions  of  projectile  time-of-flight. 

The  two  most  familiar  fire  control  systems  use  (1)  target  velocity  and  (2)  target 
velocity  and  acceleration  combined  with  time-of-flight  to  determine  the  gun  line 
offset.  Recent  application  of  sub-optimal  estimation  methodology,  specifically 
Kalman  filtering  has  resulted  in  the  development  of  velocity  plus  acceleration 
controlled  offsets  that  have  significantly  increased  the  performance  of  gun  fire 
control  systems  when  maneuvering  targets  are  engaged.  A  penalty  is  associated 
with  the  estimation  of  accelerations.  Whereas,  it  is  relatively  simple  to 
estimate  velocities,  such  is  not  the  case  for  accelerations.  As  projectile 
flight  times  increase  the  benefits  derived  from  acceleration  estimates  are 
penalized  by  the  errors  introduced  in  the  offset  calculation  by  the  "noise". 
Current  tactical  mobility  indicates  that  maneuvering  targets  require  target 
acceleration  estimates  in  addition  to  target  velocity  estimates  to  calculate 
effective  gun  line  offsets.  As  projectile  time-of-flight  increases  the  benefit 
of  this  offset  philosophy  is  degraded. 

When  the  firing  rate  of  gun  is  relatively  low  with  respect  to  the  cyclic 
motion  of  the  maneuvering  target,  it. is  possible  to  correlate  the  firing  times 
with  the  zero  crossing  of  the  velocity  perpendicular  to  the  line-of-sight.  The 
closest  approach  of  the  projectile  to  the  target  a  time-of-flight  later  crosses 
from  a  positive  to  a  negative  miss  distance  in  this  region  and  both  velocity  and 
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acceleration  are  used  to  determine  the  gun  line  offset.  The  LOS  velocity  zero 
crossing  times  can  be  sensed  with  rate  gyros  on  fire  control  system.  The 
apparent  displacement,  velocity  and  acceleration  of  a  maneuvering  target  is 
viewed  from  a  fire  control  system  situated  at  some  distance.  The  LOS  rate  zero 
crossings  are  detected  by  a  rate  gyro,  causing  the  gun  to  fire  a  projectile. 

The  closest  approach  miss  distance  between  the  target  and  projectile  a  time-of- 
f light  later  indicate  the  zero  crossing  character  of  the  closest  approach  miss 
distance  in  this  region.  Burst  fire  will  be  more  effective  in  this  region. 
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DETECTION  AND  INITIATION  OF  FIRING  COMMANDS  FOR  AN  ACCELERATION  PREDICTOR 
FIRE  CONTROL  SYSTEM  ENGAGING  MANEUVERING  TARGETS 

H.  H.  Burke  -  AMSAA 

(  DEVELOPMENT  OF  A  FIRE  CONTROL  SYSTEM  CONCEPT  TO  IMPROVE  PERFORMANCE  AGAINST 

MANEUVERING  TARGETS ) 

An  investigation  of  the  performance  of  different  order  predictors 
for  gun  fire  control  systems  engaging  maneuvering  targets  was  conducted  ^in 
1977  and  reported  on  in  Reference  1.  The  intent  of  that  study  was  to 
demonstrate  the  theoretical  reduction  in  miss  distance  that  was  obtainable, 
assumming  perfect  tracking  of  a  maneuvering  target.  The  notion  that  fire 
control  systems  should  consider  second  order  or  acceleration  terms  in  the 
prediction  of  the  gun  line  was  introduced  as  a  possible  alternative  to 
the  present  first  order  or  velocity  term  used  in  current  fire- control 
systems.  In  1978  a  paper  presented  at  AORS  XVTI  (Ref  2)  extended  this 
work  and  discussed  progress  in  the  technical  area  of  analytical  describing 
maneuvering  targets  for  the  purpose  of  studying  fire  control  system  performance 
The  three  processes  taking  place  in  all  fire  control  systems  were  identified 
as  tracking,  estimation  and  prediction  (There  is  a  fourth  process,  pointing 
which  was  added  later).  The  three  basic  known  types  of  fire  control  systems; 
manual,  disturbed  reticle,  and  stabilized  sight -director  were  described 
and  compared  in  a  basic  manner.  Prediction  characteristics  of  first  and 
second  order,  systems  were  discussed  and  the  relative  magnitude  of  fire 
control  system  errors  from  the  tracking,  estimation,  and  prediction 
processes  were  described.  A  brief  mention  of  a  firing  doctrine  utilizing 
second  order  prediction  to  maximum  advantage  was  presented  (This  idea 


11) 

(2) 


AMSAA  TR  234,  Aug  78. 

Given  at  XVII  AORS,  Nov  78. 


was  not  exploited  until  recently  and  later  in  this  paper  it  will  be  shown 
that  the  fundamental  utilization  of  this  firing  doctrine  is  the  concept 
that  can  be  easily  mechanized  and  applied  to  tank  gun  fire  control  systems, 
with  results  that  offer  major  performance  improvement  over  existing  first 
order  prediction  concepts)  A  paper  describing  the  application  of  sub-optimal 
state  estimation  design  methodology  to  the  development  of  improved  fire 
control  systems  was  the  subject  of  Reference  3.  Some  of  the  earlier  work 
was  contained  in  this  paper  along  with  the  introduction  of  the  stability 
characteristics  of  the  two  main  fire  control  concepts;  disturbed  reticle 
and  stabilized  sight -director.  It  was  argued  that  the  inherent  performance 
of  a  stabilized  sight  was  superior  to  a  disturbed  reticle  because  the  track¬ 
ing,  estimation  and  prediction  processes  were  more  decoupled  and  therefore 
each  process  could  be  designed  to  function  better.  Some  results  presented 
in  Reference  3  indicated  that  the  so  called  "lower  bound"  performance  for 
second  order  predictions  was  superior  to  first  order  predictors  for  an 
analytically  generated  path  when  tracking  noise  was  introduced  in  the  form 
of  uncorrelated  noise.  This  advantage  decreased  and  for  larger  values  of 
tracking  noise  an  inversion  in  performance  between  the  two  prediction  orders 
occurred.  This  trend  was  used  to  argue  that  the  superior  tracking  per¬ 
formance  of  a  stabilized  sight -direct or  system  would  enhance  the  benefit  of 
second  order  prediction.  Another  method  of  deciding  on  the  benefits  of  the 
two  orders  of  predictors  was  offered;  time  on  target.  It  was  also  shown  that 
for  the  analytically  modeled  target  motion  being  studied,  that  the  'basic 
shape  of  the  prediction  error  for  first  and  second  order  systems  was  dif¬ 
ferent,  and  it  was  believed  that  this  was  the  fundamental  reason  that  the 
time  on  target  indication  showed  that  for  low  tracking  errors,  second  order 

^  Given  at  Univ  of  Pitt  Modeling  £  Simulation  Symposium,  Apr  79. 


prediction  was  superior  to  first  order  prediction.  A  separate  AMSAA  Interim 
Note  (Reference  4)  was  written  in  May  79.  The  purpose  of  this  paper  was 
to  more  fully  develop  the  argument  that  predictor  orders  were  all  improtant 
in  determining  the  distribution  of  errors  in  a  fire  control  system  when 
maneuvering  targets  are  engaged.  The  target  induced  errors  were  claimed 
to  be  much  larger  then  the  tracking  and  estimation  errors  in  themselves , 
but  it  was  admitted  that  the  tracking  and  estimation  errors  had  to  be 
small  to  utilize  the  second  order  prediction,  thereby  reducing  the  target 
motion  induced  errors.  It  was  believed  and  still  is  that  the  disturbed 
reticle  fire  control  configuration  can  not  meet  the  demands  placed  on  it 
that  second  order  prediction  requires.  This  will  have  to  be  studied  in 
detail  before  a  positive  conclusion  can  be  reached. 

Throughout  the  period  of  this  work  we  have  communicated  with  many 
different  groups,  both  within  and  outside  the  government.  One  of  the  items 
we  have  shared  with  them  is  a  data  tape  having  six  time  histories  of  man¬ 
euvering  targets.  Initially  our  work  centered  on  these  paths  but  in  order 
to  remove  the  instrumentation  noise  (the  data  gathering  system  had  no 
relationship  to  a  fire  control  tracking  process,  but  was  rather  a  ranging 
system)  the  paths  were  smoothed  to  surpress  the  high  frequency  noise, 
leaving  frequencies  no  greater  than  1/8  HZ  in  most  cases.  These  data  were 
then  used  to  determine  the  " lower  bound' 9  performance  measures  mentioned 
earlier.  Based  on  these  smoothed  paths,  the  analytically  generated 
maneuvering  target  generator  was  developed.  For  the  fire  control  studies 
conducted,  to  determine  the  relative  performance  of  first  and  second  order 
predictors,  the  analytical  generated  paths  were  used. 

^ AMSAA  Interim  Note  C-82,  May  79. 
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Some  other  groups,  using  the  data  types  supplied  them  chose  to  operate 
directly  on  the  data  following  the  removal  of  some  obvious  data  timing 
irregularities.  This  approach  was  more  in  keeping  with  real  world  data, 
in  that  there  were  certainly  frequencies  greater  than  1/8  HZ  in  the  data. 

The  filtering  algorithms  used  were  therefore  required  to  cope  with  correlated 
high  frequency  noise  superimposed  on  the  basic  path  motion  that  the  analytically 
generated  paths  contained.  The  results  obtained  from  these  studies  were 
not  as  favorable  toward  second  order  prediction  as  the  AMSAA  studies.  In  an 
attempt  to  sort  out  the  reasons  for  this  difference  in  reuslts  and  to 
determine  the  mechanism  that  causes  the  performance  trends  to  shift  when 
empirical  data,  such  as  the  empirical  data  tapes,  are  used  compared  to  the 
analytically  generated  paths ,  the  same  data  were  analyzed  by  us ,  using  the 
optimal  filtering  designs  we  had  developed  using  the  maneuvering  path 
generator.  Before  discussing  the  findings  of  this  effort  it  should  be 
mentioned  that  in  Reference  1,  for  the  "lower  bound"  studies,  it  was  found 
that  if  the  cyclic  period  of  the  target  motion  was  twice  as  rapid  as  the 
time -of -flight  of  the  projectile,  an  inversion  between  first  and  second 
order  predictors  occurred.  Three  of  the  paths  reported  on  in  Reference  1 
had  this  characteristic  and  the  results  are  described  in  Reference  1. 

It  was  argued  that  realistic  maneuvering  targets  would  not  exhibit 
frequencies  in  this  region  (/£  HZ).  The  empirical  data  contained  in  the 
subject  tapes  supplied  to  the  other  researchers  and  now  being  studied  by 
us  probably  exhibited  such  frequencies,  not  from  actual  vehicle  movement, 
but  from  instrumentation  imperf ections .  This  was  the  basic  reason  for 
development  of  the  analytical  generator,  but  one  of  its  short  comings  may 
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be  that  the  correlated  high  frequency  tracking  errors  are  not  realisticly 
represented.  A  counter  argument  to  this  idea  is  that  the  data  contained 
on  the  tapes  does  not  reflect  the  correlated  tracking  error  signal  for  a 
given  fire  control  system;  be  it  a  disturbed  reticle  or  stabilized 
sight-director  system. 

However  several  of  the  paths  studied  by  some  of  the  other  researchers 
have  recently  been  processed  by  AMSAA  using  the  same  estimation  model  oath 
developed  using  the  analytical  path  generator*  Some  extremely  interesting 
and  exciting  results  have  emerged.  Initially,  time-on— target  was  used  as 
a  figure  of  merit.  Indeed,  as  some  of  the  other  investigators  have  pointed 
out,  the  second  order  predictor  does  worse  than  the  first  order  in  some 
regions  and  better  in  other  regions  of  the  paths.  Where  is  the  cue  that 
controls  this  perf ormance  reversal?  Inspection  shows  that  as  the  cyclic 
movements  referred  to  in  Reference  1,  are  less  than  twice  the  flight  time 
the  second  order  system  out  performs  the  first  order  system  and  in  periods 
where  the  cyclic  motion  becomes  more  prononunced  at  higher  frequencies, 

(i.e.  higher  apparent  accelerations)  the  first  order  system  out  performs 
the  second  order  system.  This  explains  the  reversal  of  performance 
obtained  using  the  empirical  data.  The  smoothed  data  having  only  the 
basic  path  motion,  with  the  instrumentation  noise  removed,  did  not  exhibit 
a  reversal  (the  frequencies  in  the  data  were  such  that  the  apparent  periods 
were  always  greater  than  twice  the  time-of -flight. ) 

It  was  still  believed  that  in  spite  of  this  turn  of  events?seen  from 
analyzing  the  empirical  tapes,,  that  some  advantageous  trade-off  existed  for 
second  order  predictors  over  first  order  predictors  using  the  empirical 
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path  data.  This  was  based  on  the  fact  that  the  miss  distance  time  histories 
for  the  two  predictors  were  entirely  different  for  a  given  path.  This  trend 
is  demonstrated  graphically,  in  Figure  5  of  Reference  3  for  the  analytical 
path  generator  studies  and  is  repeated  for  all  of  the  empirical  data  analyzed 
from  the  data  tape  mentioned.  Examples  of  this  are  shown  in  Figure  1  and 
appendices  2  £  3.  As  discussed  in  Reference  3,  the  trend  shown  in  Figure  1 
is  not  as  clear  cut  when  larger  values  of  noise  are  introducted,  but  it 
does  still  exist  in  practice  theory  according  to  the  empirical  paths  we  hava  analyzed 
Thus  far  this  trend  does  persist  with  suprisingly  accurate  timing.  This 
situation  speaks  well  for  the  model  we  are  using  for  the  sub-optimal  estima¬ 
tor  and  for  the  fact  that  it  was  developed  using  an  analytical  path  generator 
model  rather  than  empirical  data.  The  main  effects  described  in  Figure  5 

m 
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of  Reference  3  and  Figure  1  of  this  paper  are  (1)  the  miss  distance  for  a 
second  order  predictor  passes  through  zero  a  time-of-flight  increment  after 
the  apparent  velocity  sensed  by  the  fire  control  system  passes  thru  zero 

(2)  the  miss  distance  for  a  first  order  fire  control  system  passes  thru  zero 

a  time-of-flight  after  the  apparent  velocity  passes  thru  its  maximum  value  and 

(3)  when  the  time-of-flight  exceeds  the  time  remaining  before  the  apparent 
acceleration  reverses  direction,  the  miss  distance  errors  for  a  second  order 
predictor  becomes  very  large. 

These  findings  have  been  found  to  exist  in  the  empirical  data  and  can 
be  utilized  to-  develop  a  fire  control  system  concept  that  will  greatly 
improve  performance  against  maneuvering  targets,  over  existing  fire  control 
systems  using  no  prediction  or  l'st  order  prediction.  The  concept 
has  the  inherent  built  in  ability  to  revert  to  first  order  automatically 
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and  become  second  order  when  required.  It  is  based  on  the  firing  doctrine 
concept  mentioned  in  Reference  2  which  was  ignored  by  us  for  a  long  time 
in  favor  of  the  time  on  the  target  advantage  that  second  predictor  provided 
when  analytical  paths  were  studied.  Noise  detracts  from  the  comparison. 

The  cue  or  sensing  device  required  is  a  rate  gyro  to  sense  LOS  rate.  As 
the  LOS  rate  crosses  zero,  the  fire  control  system  will  fire  and  the  projec¬ 
tile  will  impact  the  target  a  time  of  flight  later.  The  timing  accuracy  of 
these  "firing  windows"  is  very  accurate,  even  when  the  empirical  data  tapes 
are  used.  In  some  cases  there  is  less  time  interval  in  the  "firing  window" 
than  others,  but  for  the  majority  of  opportunities,  based  on  the  logic  of 
the  LOS  rate  passing  thru  zero,  the  time  intervals  vary  between  100  ms  and 
300  ms.  When  the  firing  rate  is  in  the  vicinity  of  1  round/8  sec  for  tank 
guns,  this  rationale  is  reasonable,  in  that  the  best  opportunity  to  hit  the 
maneuvering  target  is  determined  by  utilizing  these  "firing  windows". 

When  the  target  is  not  maneuvering  too  much,  the  accleration  levels  occurring 
at  these  regions  where  the  LOS  rate  passes  thru  zero  are  not  significant 
and  do  not  corrupt  the  basic  goodness  of  first  order  prediction.  Analysis 
of  the  empirical  data  seems  to  indicate  this  effect. 

It  is  believed  that  these  findings  justify  the  incorporation  of  second 
order  prediction  in  fire  control  systems,  and  that  instead  of  assuming  a 
"fire  hose"  or  arbitrary  firing  time  to  compare  different  predictors,  that 
a  firing  doctrine  similar  to  the  one  just  described  be  adopted.  This 
concept  degrades  gracefully  to  the  manual  offset  or  no  predictor  case 
ar  the  maneuver  ievel  lessens  and  both  velocity  and  acceleration  are  small. 
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The  inherently  better  tracking  capability  of  a  stabilized  sight -director 
system  will  provide  a  more  precise  sensing  of  the  "firing  window"  than  can 
be  obtained  with  a  disturbed  reticle  system.  Also  this  firing  doctrine  can 
be  used  to  advantage  for  a  rapid  burst  fire  gun  by  firing  the  burst  in  this 
"firing  window". 

For  cyclic  frequencies  (be  they  real  maneuvering  vehicle  movements  or 
tracking  line-of -sight  motion)  there  is  a  penalty  for  higher  frequencies  and 
extended  time-of-flight  as  shown  in  appendix  1.  This  has  implications  on  gunner 
sight  line  movement  frequency  characteristics  and  long  engagement  ranges  if  there 
is  validity  to  the  serpentine  or  cardoid  model  in  representing  real  world 
vehicle  maneuvering.  The  results  of  the  analysis  reported  in  appendices  2  £  3 
for  ATMT  empricial  paths  433  and  315  indicate  there  is  such  validity  and 
a  second  order  parabolic  predictor  can  predict  serpentine  vehicle  movement  in 
this  region  of  apparent  velocity  nulling  (region  1)  with  errors  described  in 
appendix  1.  When  the  first  order  predictor  outperforms  the  second  order 
predictor,  close  inspection  of  the  region  reveals  that  the  cyclic  frequencies 
are  larger  (-  X  HZ).  It  is  believed  that  suppression  of  these  tracking 
frequencies  will  be  a  minor  problem  to  cope  with  by  the  fire  control  designer. 

HAROLD  H.  BURKE 
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APPENDIX  1 


2fnd  ORDER  PREDICTOR  AND  SERPENTINE  MANEUVERING  VEHICLE 

MOVEMENTS 


The  maneuvering  path  generator  is  built  to  generate  movement  of  a  maneuvering 
vehicle  from  point  A  to  B  as  though  the  vehicle  moves  on  arcs  of  circles 


10% 


v 


If  this  movement  is  viewed  from  a  distance,  such  that  it  is  approaching 
the  FCS,  and  if  the  velocity  of  the  vehicle  is  constant  thru  arcs  Ql  and  a2. 


the  apparent  positions,  velocities  and  accelerations  perpendicular  to  the  LOS 
can  be  assumed  to  be  simple  harmonic  motion.  Points  1  and  2  are  the  places 
where  the  apparent  position  reverses  direction,  apparent  velocity  crosses 
zero  and  apparent  acceleration  peaks.  A  detailed  discussion  of  this  set 
of  conditions  is  contained  in  a  paper  delivered  at  AORS  XVII,  Nov  79.  This 
"serpentine"  math  model  has  been  in  a  sense  validated  by  the  fact  that  for 
the  AT MX  empirical  data  tapes,  at  least  in  the  regions  surrounding  points 
1  and  2,  these  apparent  position,  velocity  and  acceleration  trends  exist. 

The  next  question  to  address  is  how  well  does  a  second  order  or  velocity 
plus  acceleration  offset  predictor  work  for  this  region  about  1  S  2. 

Assuming  the  analytic  model  of  apparent  vehicle  motion  from  the  maximum 


position  point  at  t*0,  the  equation  relating  apparent  movement  to  time 
is  y  =  —a  +  acoswt 


Representing  cos  dc.  in  a  series,  we  have 


where  a  »  amplitude  or  radius  of  circle 
u  »  cyclic  frequency  of  movement 
t  =*  time 


cos  at  =  1 


(at)  2  (at)  4 

2!  41 


r  a(  (it)  2  a(  ht)  *  n 

Then  y  =  -a  +  ‘•a-———  +  — 4T“  J 


or 


y  *  -(«.*)  -  ♦  —  n- 


But  am2  s  magnitude  of  apparent  acceleration  at  point  1  (A  ) 


Then 


or 


Aapp  ^  Aapp  tj  2 


2  12 


Aapp 


t2  IT  -  2  7 

2  [  12  \ 


The  apparent  movement  of  y  for  the  serpentine  math  model  of  the  path 
generator  from  point  1  (t=0)  approximates  the  parabolic  second  order- 
predictor  concept,  ie: 
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yt0  *•  t  -App  r 


7to  “  0  t 
en  y  a  -App  - 


with  the  w*t2  tern  representing  the  error  between  the  simple  harmonic  and 
12 

parabolic  curves. 

To  obtain  some  insight  into  the  magnitude  of  this  error  the  following  table 
shows  the  error  term  for  various  cyclic  frequencies  and  times  of  flight 


t  (time  of  flight)  1  sec  1.5  sec  2.0  sec 

w  a  1/4  HZ  (1.5  r/s)  0.1875  0.421  0.73 

m  *  1/8  HZ  (.785  r/s)  0.05  0.078  0.21^ 

m  a  1/10  HZ  (.628  r/s)  0.033  0.049  0.13 

Plotted,  this  demonstrates  the  degrading  relationship  that  u (cyclic 
frequency)  has  on  the  ability  of  a  second  order  predictor  to  predict 
simple  harmonic  motion  in  this  localized  region  (point  1). 


APPENDIX  2 


Inclosure  la 
ATMT  Path  '433 

Tp  *  1.3  sec.;  Filter:  Q  »  8.,1.,  R  »  .01 


(1) 

Firing  Occurs  At 


2.9s 

l.atof 

4.4s 


Vv»3.2  M/S 

A 

Vy«4.3  m/s 

v  _ 

3a  taxi  Jf  a  39.a 

V„ 


Miss  Distance 


2'nd 

Order 

1st  Order 

(M) 

Ol) 

4.2 

3.09 

4.66 

4.3 

1.62 

3.98 

4.4 

0.19 

3.25 

4.3 

-0.48 

2.80 

4.6 

-0.17 

2.90 

4.7 

-0.28 

2.78 

(2) 


Firing  Occurs  At 

12.7 

12.8 

+0.196 

-1.94 

11.5s 

12.9 

-  .484 

-2.44 

l.Stof 

13.0 

-  .69 

-2.55 

13.0s 

13.1 

-  .81 

-2.59 

13.2 

-1.23 

-2.91 

13.3 

-1.69 

-3.26 

13.4 

-1.94 

-3.35 

Vx-7.5 

Vy-2.2 
9  a  16.3°' 
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AIMT  PATH  433 


(3  ) 

Miss  Distance 

Firing  Occurs  At 

2nd  Order 

1st  Order 

IX) 

l  M) 

16.2s 

1.5tof 

17.4  .338 

2.81 

17.7s 

17.5  -  .491 

2.64 

17.6  458 

2.57 

17.7  -  .253 

2.85 

17.8  +  .148 

3.17 

17.9  +  .766 

o .  o  2 

18.0  +  .846 

3.57 

VX  =  5*5 


Vy  =  3.5 
9  =32.5° 


(4  ) 

26.0 

-4.891 

-7.21 

Firing  Occurs  At 

26.1 

-4.692 

-7.20 

26.2 

-4.216 

-6. 95 

24.3s 

26.3 

-2.12 

-5.77 

1.5tof 

26.4 

-  .072 

—4 .  o  3 

2Q.3s 

to 

• 

to 

CO 

+2 . 06 

-3.17 

26.6 

+3.66 

-1.98 

V  =  6.6 

A 

VY  =  6.4 
3  =  44.1° 


(5  ) 

Firing  Occurs  At 


31.1 

2.11 

4.56 

30.0s 

31.2 

2;85 

4.96 

1 . 5tof 

31.3 

2.12 

4.  58 

31 .  os 

31.4 

1.08 

4.03 

31.5 

.25 

3.51 

31.6 

-  .56 

2.94 

VY  =  4.7 

31.7 

-  .98 

2.70 

X 

31.8 

-1.52 

2.34 

vv  =4.5 

31.9 

-2.36 

1.77 

9  3  43.8° 
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ATMT  PATH  433 


(  6  ) 

Firing  Occurs  At 
39.3s 
l.Stof 

40.  as 

(acceleration  reversal  during  tof) 

Vx  .  3.3 

VY  »  .121 

9  »  1.3° 


Miss  Distance 


2nd  Order 

1st  Order 

(M) 

U) 

40.2 

-1.95 

-3.25 

40.5 

-1.12 

-2.84 

40.4 

.45 

-1.96 

40.5 

1.50 

-1.28 

40.6 

1.63 

-  .996 

40.7 

1.81 

.732 

40.8 

2.42 

-  .239 

40.9 

3.80 

-  .681 

(-) 

Firing  Occurs  At 

40.7s 

1.4tof 

42.2s 


VY  *  5.4 

A 


4.2 


41.8  +  .434  1.42 

41.9  +■  .169  .98 

42.0  +  .197  8.54 

42.1  -  .041  .87 

42.2  -  .650  .74 

42.3  -1.07  .32 

42.4  -1.22  .00 

42.5  -1.18  -  .16 

42.6  -1.01  -  .19 

42.7  -1.04  -  .13 


(  3} 

Firing  Occurs .At 
57.8s 
1 . 5tof 
59.3s 


Vx  -  7.1 
VY  »  2.3 
9  »  17.9° 


58.9 

59.0 

-2.85 

-3.13 

59.1 

-2.50 

-3.04 

59.2 

-2.40 

-3.12 

59.3 

-1.81 

-2.86 

59.4 

-  .67 

-2.27 

59.5 

+1.09 

-1.22 

59.6 

+2.50 

-2.55 

59.7 

-3.26 

+  .343 
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ATMT  PATH  433 


(9  ) 

Firing  Occurs  At 
60.5s 
1. 5tof 
62.0s 


Miss  Distance 

2nd  Order  1st  Order 


IM) 

IM) 

61.6 

.687 

2.47 

61.7 

.17 

2.16 

61.3 

.46 

2.44 

61.9 

.13 

2.32 

62. Q 

— 

.26 

2.13 

62.1 

- 

.34 

2.07 

62.2 

- 

.11 

2.17 

62.3 

+ 

.26 

2.29 

62.4 

.75 

2.54 

VX  "  5*7 


V  * 
Y* 

a  - 


2.3 

18.9° 


(io) 

Firing  Occurs  At 


67.0s 

68.1 

-1.08 

-3.39 

1. 5tof 

68.2 

-1.56 

-4.19 

68.5 

68.3 

-2.87 

-4.98 

68.4 

-2.99 

-5.24 

68  •  5 

-3.45 

—5 . 53 

68  •  6 

-3.57 

—o .  7o 

68.7 

-3.46 

-5.34 

'  68.8 

CM 

. 

CM 

1 

-5.70 

VX  =  6,7 


VY  =  2.3 


19.4 


(11) 


Firing  Occurs  At 

77.8 

2.42 

4.12 

76.7s 

77.9 

1.83 

3.84 

1 . Stof 

78.0 

1.33 

3.51 

78.2S 

78.1 

1.19 

3.57 

78.2 

1.68 

3.85 

78.3 

1.78 

3.35 

78.4 

0.47 

3.17 

V  =8.2 

78.5 

-  .380 

2.51 

A 

78.5 

-1.15 

1.95 

21.3° 


9  = 
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ATMT  PATH  433 


( 12)  Miss  Distance 

Firing  Occurs  At  2nd  Order  1st  Order 


81.23 

MiMI 

(M) 

(a) 

l.Stof 

82.2 

—1.89 

-3.41 

82.7s 

82.3 

-1.71 

-3.36 

82.4 

-1.38 

-3.30 

82.5 

-  .44 

'  -2.78 

82.6 

+  .23 

-2.34 

82.7 

.14 

-2.50 

82.3 

-  .147 

-2.70 

82.9 

-  .53 

-3.26 

83.0 

—1.67 

-3.43 

Vx  *  6.9 

Vy  a  3.6 

6  =27.6° 

83.1 

-2.06 

-3.51 

(  13) 

Firing  Occurs  At 

85.4 

.521 

2.35 

84.4s 

85*5 

—  .560 

1.88 

l.Stof 

85.6 

-2.26 

.771 

85.9s 

85.7 

3.43' 

-  .019 

85.8 

3.03 

-  .174 

(Velocity  reversal  during  tof) 

85.9 

-2.55- 

.012 

86.0 

1.195 

.695 

86.1 

.044 

1.13 

Vv  =  7.7 

•  86.2 

1.028 

1.69 

X 

VY  =  1.4 

'9  -  10.3° 

86.3 

1.18 

1.74 

(14)  .  ■ 

Firing  Occurs  At 

-2.69 

91.6s 

92.8 

-  .74 

1 . 5tof 

92.9 

-  .88 

-2.75 

93.1s 

93.0 

-1.52 

-3.16 

93.1 

—1  •  82 

-3.35 

93.2 

-1.66 

-3.37 

93.3 

-1.12 

-3.14 

Vx  =  8.3 

l* :  s-.9f 

93.4 

-  .847 

-3.04 
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APPENDIX  3 


t*0 


METERS 


?CS  location 
x-59150  *" 
/■73250 


ATMt  Path  315,  Twister  Vehicle 


Inclosure  lb 


ATMT  315 

T_  *  1.5  sec 
r 


X 


4.8 


(1) 

Firing  Occurs  At 

Miss  Distance 

7.8s 

2nd  Order 

1st  Order 

* 

1. 5tof 

IM) 

IM) 

9  *  os 

8.9 

9.0  -2.241 

■  -  2.98 

9.1  -1.996 

-  2.88 

9.2  -1.821 

-  2.36 

9.3  -1.022 

—  2.53 

9.4  —  •  o  /  4 

-  2.31 

9 . 5  -  .247 

-  2.26 

9.6  -  .382 

-  2.42 

9./  —  • ooo 

-  2.58 

VY  =  2.9 
■3  =  31.  l' 


(2) 

Firing  Occurs  at 


15. 6s 

16.8 

1.5tof 

16.9 

-1.99 

.793 

17.1s 

17.0 

-2.67 

.258 

17.1 

-2.52 

.220 

17.2 

-1.72 

.613 

17.3 

-1.48 

•  683 

17.4 

-1.29 

.714 

17.5 

-  .601 

#969 

17.6 

+  .449 

1.43 

17.7 

+1.692 

2.021 

Vx  -  8.2 


=  1.4 
9  7° 

9  »  9*7 
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AT!CT  PATH  315 


(  3) 

Firing  Occurs  At 


24.3s 

1.5tof 

25.8 


Vx  =  7.2 


V„  =  3.2 

1  ■ 

3  =  24.0° 


Miss  Distance 


2nd 

Order 

1st  Order 

IMJ 

25.3 

25.4 

IM) 

25.5 

-2.259 

+3.46 

25.6 

-1.520 

+3.09- 

25.7 

-1.210 

+3.01 

25.8 

-1.042 

+2.96 

25.9 

-1.099 

+3.00 

26.0 

-  .812 

+2.86 

26.1 

-  .  866 

+2.892 

26.2 

-1.47 

+3.17 

(4  ) 

Firing  Occurs  At 

29.4s 
1 .  atof 
30.9 


V.c  =  9.6 


0 


30.4 

0  0 . 0 

+  .969 

1.165 

30.6 

*  +  . 884 

1.092 

30.7 

.557 

.892 

30.8 

+  .208 

.721 

30.9 

-  .611 

.284 

31.0 

-1 . 65 

-  .322 

31.1 

-2.218 

-  .707 

31.2 

-2.22 

-  .799 

31.3 

31.4 

0 1 . 0 

(  5  ) 


Firing  Occurs  At 

40.0 

-1.226 

-1.92 

38.9s 

40.1 

-1.077 

-1.86 

1.5tof 

40.2 

-1.452 

-2,06 

40.4s 

40.3 

-  .945 

-1.79 

40.4 

-  .313 

-1.39 

40.5 

+  .143 

-1.12 

Vv  =  7.2 

40.6 

+  .675 

-  .791 

X 

40.7 

+1.560 

.257 

V-  =  1.3 

40.3 

+2.75 

+  •  0  27 

3  =  10.2° 
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ATMT  PATH  315 


(  5  )  Miss 


Firing  Occurs  At 

2nd  Order 

44.2s 

(M) 

l.Stof 

45.3 

1.28 

45.7s 

45.4 

1.17 

45.5 

.647 

45.6 

•  •  2o  o 

45.7 

-  .966 

45.8 

-1.4 

45.9 

-1.53 

46.0 

-1.77 

V,c  =  1.5 

V,.  =  .25 

i 

0 

3  •  =  9.5 

(  7  ) 

Firing  Occurs  At 

59.1s 

60.2 

-2.3 

1 • 5tof 

60.3 

-2.7 

60.6 

60.4 

to 

. 

CM 

1 

60.5 

-2.1 

60.6 

-1.3 

60.7 

-  .411 

60.8 

+2.12 

Vx  =  4.7 

Vy  »  1.9 

,  _  22.0° 
c  — 

60.9 

+  .472 

(8  ) 

75.0 

.514 

Firing  Occurs  At 

75.1 

.324 

74.0s 

75.2 

1.16 

1. 5tof 

75.3 

1.51 

75.  os 

75.4 

1.24 

75.5 

.974 

75.6 

-  .089 

75.7 

-  .786 

Vv  =  8.0 

75.8 

-1.25 

X 

75.9 

-1.92 

VY  -  1.7 

Distance 

1st  Order 

[IV 

•  977 
.965 

•  598 
.239 

-  .142 

-  .455 

-  .620 


2*  365 
2.369 
2.254 
1.94 
1.46 
1.12 
.908 
.  655 


1.86 
1.74 
2.21 
2.47 
2.44 
2.37 
1.73 
1.31 
1.05 
.  55 


3 
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ATMT  PATH  315 


(9  ) 

Firing  Occurs  At 
76.6s 
1 . 5tof 
73.1s 


V,,  =  4.4 
3  =  35.4° 


Miss  Distance 

2nd  Order  1st  Order 


01) 

on 

77.6 

— 6 . 35 

6 . 30 

77.7 

-5.23 

6 . 00 

77.8 

—4. 1 5 

5.77 

77.9 

-2.29 

4.92 

78.0 

-  .740 

4,11 

78.1 

-  .025 

3*71 

78.2 

-  .244 

3.86 

78.3 

-  .656 

3.96 

78.4 

-  .665 

3.90 

78.5 

+  .020 

3.51 

( lo) 

Firing  Occurs  At 


81.1s 

82.0 

o .  42 

6 . 57 

1. otof 

82.1 

3.36 

6.32 

82.6s 

82.2 

2.41 

6.71 

82.3 

1.48 

6.21 

82.4 

-  .410 

5 . 33 

82.5 

-1.50 

4.63 

32.6 

-1.58 

3.89 

'  82.7 

-1.13 

4.04 

82.8 

-  .314 

4.08 

32.9 

33.0 

VY=-  5.0 

1  -  0 
8  s  45.0 


(ID 

Firing  Occurs  At 
84.3s 
1. 5tof 
85.8s 


V„  *  8.5 


V  =  5.0 

v  o 

*  -  30.2 


85.4 

-3.84 

85.5 

-2.011 

85.6 

-  .673 

85.7 

+1.711 

85.8 

+3.46 

85.9 

+4.71 

86.0 

+5.049 

36.1 

4.93 

86.2 

-7.302 
—6 ,251 
-5.40 
-4.01 

-  2.37 

-  1.90 

-  1.23 

-  .937 


3 
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ATMT  PATH  315 


(  12) 

Firing  Occurs  At 

87.5s 
1 . 5tof 
89.0 


Miss  Distance 

2nd  Order  1st  Order 


(M) 

(M) 

38.6 

3.628 

5.938 

38.7 

2.963 

5.695 

88.3 

1.981 

5.331 

88.9 

.850 

4.722 

89.0 

-  .678 

3.768 

89.1 

-1.49 

3.197 

89.2 

-1.88 

2.817 

89.3 

-2.13 

2.542 

89.4 

89.5 

-3.47 

1.630 

=  6.6 

Vy^4.0 


(  is) 


Firing  Occurs  At 

93.0 

l 

to 

• 

cn 

-5.86 

91.9s 

93.1 

-1.68 

-5.237 

1 . 5tof 

m 

93.2 

-1.29 

— o . Qoo 

93.4s 

93.3 

-  .743 

-4.603 

93.4 

+1.72 

-3.088 

93.5 

3.03 

-2.127 

93.6 

4.79 

-  .816 

93.7 

— =-r024 

*  93.8 

+  .789 

V  =  8.4 


(  14) 


Firing  Occurs  At 

97.1 

5.2 

4.98 

- 

96.1s 

97.2 

5.6 

5 . 46 

1.5tof 

97.3 

5.4 

5.49 

97.6s 

97.4 

4.6 

5.28. 

- 

97.5 

3.8 

5.02. 

97.6 

2.3 

4.23- 

97.7 

.934 

3.43 

V„  =  1.6 

97.8 

•  o  2o 

3.23 

97.9 

-  .678 

2.82 

VY  =  5.0 

98.0 

-1.67 

2.26 

n 

98.1 
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HYBRID  COMPUTER  SIMULATION  OF 
COMBAT  TANK  DRIVEN  RETICLE  FIRE  CONTROL 

By  John  N.  Groff 

US  Army  Materiel  Systems  Analysis  Activity 


1.  INTRODUCTION 

The  US  Army  Materiel  Systems  Analysis  Activity  (AMSAA)  and 
the  Ballistic  Research  Laboratory  (BRL)  have  jointly  developed  a  man-in- 
the-loop  hybrid  computer  simulation  of  a  combat  tank's  turret/weapon 
stabilization  drives  and  those  portions  of  the  fire  control  that  are 
necessary  for  lead  angle  generation.  The  model  that  has  been  imple¬ 
mented  on  the  BRL  EAI  690  Hybrid  Computer  is  linear  although  it  does 
possess  the  following  non-linear  features: 

•  Hand  Control  Dead  Zones 

•  Response  Rates  versus  Handle  Control  Deflection 

Currently,  the  simulation  is  still  undergoing  development  of 
the  gunner's  oscilloscope  display,  scaling  of  the  input/output  variables, 
and  programming  of  the  time  series  software  necessary  to  reduce  the  data. 
However,  qualitative  rather  than  quantitative  information  concerning 
gunner  tracking  performance  is  being  provided  by  the  simulation. 

2.  BACKGROUND 

In  1978  AMSAA  and  BRL  jointly  undertook  the  responsibility  of 
investigating  whether  proposed  modification  to  a  driven  reticle  continuous 
lead  insertion  system  had  a  significant  effect  in  reducing  weapon 
pointing  errors.  The  basic  problem  was  a  response  mismatch  between  the 
reticle  projection  unit  (RPU)  of  the  gunner's  sight  and  the  turret/weapon 
stabilization  drive  which  resulted  in  the  weapon  lagging  the  RPU.  The 
resulting  offset  appeared  in  the  sight  as  apparent  tracking  error  which 
the  gunner  attempted  to  null  out.  The  problem  tended  to  be  further  ag¬ 
gravated  if  evasive  or  sinusoidal  target  motion  was  introduced.  This 
tended  to  produce  an  oscillatory  instability  or  “rubberband"  tracking 
error  effect.  The  investigation  took  the  form  of  a  paper  analysis 
using  available  system  description  documentation.  The  results  of  this 
paper  analysis  are  presented  in  Section  3  and  served  as  the  basis  for 
development  of  the  AMSAA/BRL  Hybrid  Simulation.  Further  motivation  for 
the  hybrid  work  was  provided  by  AMSAA' s  delivery  accuracy  efforts. 

AMSAA  has  the  responsibility  of  determining  quasi-combat 
hitting  probabilities  for  armored  vehicles,  and  the  use  of  engineering 
simulation  techniques  have  proved  to  be  useful  in  performing  this 
work.  Currently,  engineering  simulations  for  the  M60A1,  M60A3,  and 
XM1  tanks  are  being  developed  in  which  all  major  components  of  these 
tanks'  fire  controls,  weapon  stabilization  drive  loops,  suspension 
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system  as  well  as  terrain  profile  and  human  characteristics  will  be 
modeled  and  integrated  into  a  single  simulation. 

Previous  AMSAA  man-in-the-loop  studies  have  indicated  that 
the  gunner  model  used  in  these  simulations  may  not  be  entirely  satis¬ 
factory  for  delivery  accuracy  work.  Part  of  the  problem  stems  from 
the  fact  that  the  McRurer/Krendel  gunner  model  which  will  be  used  does 
not  possess  adaptable  features,  i.e.,  the  ability  of  the  gunner  model 
to  make  internal  adjustments  in  gain,  bandwidth,  and  neuromuscul ar 
reaction  time  in  response  to  target  motion  and  characteristics  of  the 
tracking  system  in  which  the  gunner  is  operating.  Another  inadequacy 
of  the  present  model  is  that  it  is  a  purely  deterministic  model  lacking 
any  residual  randomness. 
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3.  RESULTS  OF  THE  PAPER  ANALYSIS 


Figure  1  depicts  the  block  diagram  of  the  driven  reticle 
system  that  is  being  used  in  this  study. 
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3.  CLOSED  LOOf  CONNGURATIONS: 

a.  FIRST  ORDER  LAG  DEMODULATOR  0,(3) 

b.  SECONO  ORDER  QUADRATIC  DEMOOUtATOR,  0,(3)  »  Q  (Co.«o> 


FIGURE  1.  DRIVEN  RETICLE/CLOSED  LOOP  CONFIGURATION 

A  McRurer/Krendal  man-model  was  interfaced  with  this  system  and  used 
in  the  analysis. 

The  general  human  transfer  function  used  is  of  the  form: 


-TS 

Kq  (S  +  u>2)  e 

G(s)  =  J - 

(s  +  a»l)  (s  +  (1)3) 


(3.1) 


Where,  t  is  a  neurological  time  delay;  K,  103  are  values 

that  the  human  adjusts  for  the  task  at  hand.  Typically,  he  adjusts 
these  to  achieve  a  loop  crossover  frequency  of  3  rad/sec  with  a 
phase  margin  of  35  to  45  degrees.  (1)3  is  associated  with  the 
human  bandwidth. 


For  this  investigation,  the  following  nominal  values  were 
used;  (dj  =  0.0  rad/sec,  =  0.667  rad/sec,  013  =  6.667  rad/sec, 
and  t  =  0.1  sec.  The  human  gain,  K,  was  established  on  the  basis  of 
subsequent  analysis. 

Specifically,  open  loop/root  locus  techniques  were  applied  to 
the  overall  system  for  the  purpose  of  determining  relative  stability  of 
this  system  with  the  gunner  model  included.  On  the  basis  of  these 
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analyses,  the  open  loop  gain  was  adjusted  to  provide  a  gain  margin  of 
3DB  at  the  crossover  frequency.  Figure  2  depicts  a  further  simplifi¬ 
cation  of  the  overall  system  shown  in  Figure  1.  Essentially,  the  TACH 
feed-ahead  loop  shown  in  Figure  1  has  been  eliminated  and  replaced  with 
two  parallel  branches;  an  RPU  loop,  and  a  Turret  drive/stabil ization 
loop.  In  addition,  the  common  gain  in  both  these  loops  has  been  com¬ 
bined  with  the  human  gain,  optical  sight  gain,  and  the  handle-bar  gain 
to  yield  a  total  system  gain,  Ks.  The  open  loop  transfer  function  for 
the  system  depicted  in  Figure  2  is  defined  by  Equation  3.2, 


8g-dl 


Figure  2.  SIMPLIFIED  AZIMUTH  BLOCK  DIAGRAM 


Where, 

G]_  (s)  =  Ks  (s  +  0.667)  EXP  (-ts)/  s(s  +  6.667) 


G2  (s)  «  (tf/(l  +  0.39s)  Q  (.3,  21.9))  ((1/  D  (s)  Q  (.45,  18.0J-1) 

/1/(1  +  0.23s)  First  Order  Demod  1 
0  (s)  •]  \ 

ll/Q  (.7,  100)  Second  Order  Demod/ 

S0i  (s)  -  Gx  (s)  (S2  (s)  +  G3  (s))  3.2 


The  subscript,  ol ,  denotes  open  loop.  The  open  loop  gain,  Ks,  has  been 
computed  for  the  cases  where  D(S)  was  either  a  first  order  or  second 
order  demod  filter,  i.e. , 


so 


D  (S) _ 

First  Order  Filter 
Second  Order  Filter 


Open  Loop  Gain 


6.3 
12.5 

With  these  gains,  the  respective  closed  loop  frequency 
responses  were  generated  for  the  systems  having  either  a  first  order  or 
second  order  demod  filter.  Figure  3  depicts  the  bode  frequency  response 
for  these  systems.  It  is  apparent  that  substantial  improvement  would 
be  realized  from  a  system  using  a  second  order  demod  filter  in  place  of 
the  first  order  filter. 
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Figure  3  Driven  Reticle /Closed  loop  Frequency  Response. 


In  order  to  quantify  this  improvement,  the  following  statistical 
techniques  employing  power  spectral  density  of  target  motion  were  used. 
Equation  3.3  represents  a  one  sided  power  spectral  density  model  of 
target  lateral  linear  acceleration,  y. 


r..  (») 

x 


2a2a  (w2  +  (a  2  +  a»c2) 
ai4  +  2(a  4  -  ucC)uc  +  (a  c  +  a)c^)2 


3.3 
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Where, 

a  and  oic  are  positive  constants 
1  -  a  *  1,2,  3  meters/sec^ 
u  =  Angular  frequency 

The  basic  technique  is  to  play  the  target  PSD  through  an  error 
transmissabil ity  relationship,  i.e., 


Qg  -  ®L  3  Gol  (s) 


Gd(s) 


1  +  Gol  (s) 

where  Gc-|  (s)  =  Closed  Loop  Transfer  Function 
Gol  (s)  =  Open  Loop  Transfer  Function 

Now, 

e  =  ©t  -  (0g  -  Q|_)  =  ®t  "  Gcl  (s)Qt 

■  -  (1  -  Gcl  (s))©t 

G01  (s) 


1  +  Gol(s) 
1 


0. 


1  +  G0-|  (s) 

Equation  3.8  may  be  rewritten  in  the  following  form 


1 


e  = 


0, 


1  +  Gol (s) 

For  sufficiently  small  0^, 

0t  -  X  • 

TT 

Equation  3.9  now  becomes 
1  1 


s7  l+Gol(s) 


3.4 


3.5 


3.6 


3.7 


3.8 


3.9 


3.10 
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3.11 


In  terns  of  PSD  relationships.  Equation  3.10  becomes, 

r  (u)  =  r..(<o)/  (1  +  G0i  (ito)  Go!  ( i <**) )  0,4 
e  X 

Integration  of  r  (<*>)  over  the  frequency  domain  yields  the  variance 
£ 

of  the  error,  i.e. , 

1  ' 

a  2  3  ___  j  r  (a»)dtu. 

ee  2ir  e  3.12 

Table  3.1  presents  the  error  one-sigma  values  for  the  two  different 
systems  as  a  function  of  lateral  acceleration.  A  1500  meter  range  was 
used. 


TABLE  3.1  CLOSED  LOOP  ERROR,  a  (R  ■  1500  METERS) 

££ 


Lateral 

First  Order 

Second  Order 

Accel erati  mi 

Demod  Filter, 

Demod  Filter 

meters /seer 

a 

a 

a 

X 

£ 

e 

1.0 

0. 280 

0.135 

2.0 

0.559 

0.270 

3.0 

0.839 

0.406 

4.  DESCRIPTION  OF  HYBRID  SIMULATION 
4.1  Oscilloscope  Display. 

The  AMSAA/BRL  simulation  presents  a  fixed  base  display  to  the 
gunner  whereas  the  actual  system  has  the  ability  to  slew  the  turret. 

In  terms  of  task  description,  the  actual  systems  might  be  viewed  as  a 
pursuit  tracking  task  while  the  simulation  is  compensatory.  Figure  4 
shows  the  oscilloscope  display  presented  to  the  gunner  which  consists 
of  a  point  target  and  a  cross  hair  reticle.  Both  the  point  target  and 
reticle  possess  the  ability  to  be  offset  from  the  center  of  the  scope 
display.  This  is  accomplished  by  commands  received  from  the  EAI680 
analog  consoles.  These  commands  may  be  generated  in  a  variety  of  dif¬ 
ferent  ways  depending  on  how  the  various  angular  offsets  depicted  in 
Figure  1  are  formed  into  different  display  commands. 
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FIGURE  4.  DRIVEN  RETICLE  DISPLAY 
Consider  the  residual  error  0r  which  is  presented  to  the  gunner, 

®R  3  ®t  +  ®L  “  ®g  4. 1 

In  terms  of  the  control  handle  output,  ©g  ,  Equation  4.1  takes  the  form; 


©O  -  ®+.  G|  (s)  Gjyjy  (s)  0 


R  =  t 


Where, 


>g  -  Gg  (s)  Gtac(s)  D(s)  Grpu(s)  0  -  Gg  (s)  0 

c  c  c 

4.2 


0R  *  Angular  Offset  between  Point  Target  and  Reticle, 
0-t  =  Angular  Offset  of  Target,  Motion 


0a  =  Commanded  Gun  Rate 
yc 


GL(s)  =  Tf  Kf 

- - ,  where  tf  is  time  of  flight. 


GTac(s) 


TjS+1 

Ktacs 


KGear 

D(S)  ■  Demodulator  Transfer  Function, 
Grpu  (s)  3  Kr/Q  Ur,  ojr),  and 
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Equation  4.2  may  be  rewritten  in  terms  of  the  commanded  inputs,  i.e. , 

Qr  =  Ci  +  C2  -  C3  -  C4  4.3 

Where, 

Cl  =  ©t. 

C2  =  ©l  =  Gj_t  (s)  Grpu  (s)  0gC, 

f 

C3  =  Gg  (s)  Gtac  (s)  Gt_t  (s)  D(s)  0gc 

C4  =  Gg(s)  0gc. 

Let  Cs  represent  the  spot  command  and  Cc^  the  cross  hair  command.  Var¬ 
ious  ways  exist  of  summing  the  various  commands  (01,02,03,04)  into  the 
spot  command.  Cs  and  cross  hair  command  Ccj,.  Three  ways  of  summing  the 
commands  are  presented  in  this  paper. 

Method  1. 

Cs  =  Ci  +  C2  -  C3  -  C4 
Cch  *  0.0,  and 
Method  2. 

Cs  -  Ci 

Cch  =  C2  -  C3  -  C4. 

Method  3. 

C  ^  =  Ci  —  C  4 

Cch  3  C2  -  C3 

The  first  method  results  in  a  purely  compensatory  type  of  tracking  task 
presented  to  the  gunner.  This  proved  unsatisfactory,  since  none  of  the 
"rubberband"  phenomenon  was  exhibited.  The  gunner  simply  nulled  out 
the  residual  error. 

The  second  mechanization  presents  a  pursuit  tracking  task  to 
the  gunner  in  which  he  must  chase  the  point  target  over  the  display 


screen:  whereas,  in  the  actual  tank  the  pursuit  track  task  results  in 
the  turret  rotating  and  the  reticle  staying  close  to  the  center  of  the 
sight  display. 

Method  two  has  not  yet  been  mechanized  and  checked  out  and 
may  prove  unfeasible  for  certain  types  of  target  motion. 

Method  three  is  currently  mechanized  and  will  be  demon¬ 
strated.  This  mechanization,  while  compensatory,  does  simulate  the 
"rubberband"  phenomenon  and  can  accommodate  all  types  of  target  motion. 

4.2  Overview  of  Hybrid  Computer  Simulation. 

Figure  5  provides  the  reader  with  an  overview  of  the  AMSAA/ 
BRL  Hybrid  Computer  Simulation. 


FIGURE  5  OVERVIEW  OF  HYBRID  SIMULATION 


All  of  the  major  computing  components  have  been  identified  and  for  the 
most  part  their  functions  are  self-evident.  However,  some  classification 
concerning  the  model  that  has  been  programmed  and  the  planned  use  of 
the  output  data  is  needed.  It  is: 
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a*  The  turret  drive  and  fire  control  model  shown  in  Figure  1 
of  this  paper  has  been  programmed  on  the  EAI  680  analog  computer.  An 
elevation,  drive  model  has  also  been  implemented.  For  the  time  being  a 
very  simple  target  path  generator  has  also  been  programmed  on  the  680. 
The  following  equation  defines  the  evasive  part  of  the  target  path 
generator. 


Ot(s)  =  G(s)  Wn(s) 


4.4 


Where 


G(s)  =  K/(ts  +  1) 

Wn  =  White  Noise 
Q-t  *  Target  Angular  Velocity 

The  constants  K  and  x  have  been  chosen  so  as  to  provide  zero  to  peak 
amplitude  of  12  mils/sec  with  a  bandwidth  of  0  to  .2  Hz.  The  simula¬ 
tion  also  can  be  exercised  against  a  constant  velocity  target. 

b.  The  analytical  objectives  of  this  effort  are: 

(1)  Confirm  the  results  obtained  from  the  paper  study 
of  Section  3.0. 

(2)  Create  a  data  base  suitable  to  refine  the  gunner 
model  currently  being  used  in  digital  engineering 
simulation. 

Initially,  it  is  planned  to  develop  a  residual  error  model  which 
describes  the  random  portion  of  the  human.  Box  Jenkins  and  Parameter 
Identification  techniques  appear  suitable  for  this  purpose.  Ultimately, 
a  simplified  version  of  the  Kleinman  (Kalman  Filter/Predictor)  model 
would  be  the  desired  goal  of  this  work. 
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INTRODUCTION: 

Traditionally  fire  control  prediction  algorithms  have  been  based  upon  third- 
order  system  models  in  each  coordinate  (1) .  In  a  previous  paper  by  the  authors 
(2),  the  concept  of  using  higher  order  autoregressive (AR)  models  was  introduced. 
It  was  found  that  in  a  noiseless  environment  these  models  provided  robust 
predictors  which  could  significantly  improve  the  capabilities  of  an  anti-aircraft 
(AA)  artillery  weapon  system  against  a  large  class  of  aircraft  maneuvers  at 
extended  times  of  flight.  This  follow-up  study  examines  the  case  where  the 
observations  made  by  the  weapon  system  are  corrupted  by  noise.  The  achieved 
results  show  that  robust  higher  order  AR  models  still  yield  considerable 
improvement,  especially  in  the  filtering  of  the  sensor  signals.  In  addition, 
an  artificial  flight  path  was  constructed  by  finding  those  autoregressive 
coefficients  which  maximize  the  variance  of  the  estimated  aircraft  velocity. 

It  was  found  that  AR  predictors  modeled  after  real  aircraft  data  work  well 
against  this  "worse  case"  flight  path. 

TEST  DATA  AND  SIMULATION  PROGRAM: 


The  system  models  discussed  here  and  in  Reference  2  were  determined  using  time 
series  analysis  from  test  data.  This  data  was  collected  by  the  US  Navy  from 
flight  paths  flown  by  an  A-7E  high  performance  aircraft  simulating  bombing  a 
defended  ground  target.  The  recorded  data  was  manipulated  and  smoothed  to 
yield  consistent  position,  velocity,  acceleration  and  acceleration-dot (the 
first  derivative  of  acceleration)  data  at  0.1  second  intervals  in  the  XYZ 
coordinate  system.  The  following  descriptions  should  make  it  clear  that 
leaving  the  fourth  and  higher  derivatives  of  position  all  undefined  does  not 
cause  any  of  the  complexities  of  the  maneuvers  to  be  lost.  Although  the 
complete  data  set  consist  s  of  eleven  different  flight  paths,  results  for 
only  three  representative  ones  will  be  presented  here.  These  will  be  the 
same  passes  that  were  used  in  Reference  2  and  the  reader  is  encouraged  to 
refer  to  this  paper  for  a  schematic  figuration  of  these  passes. 

Flight  path  #1  is  representative  of  a  general  class  of  maneuvers  known  as  a 
"dive  toss".  That  is, the  payload  is  released  while  the  aircraft  is  rolling 
and  pitching  which  gives  the  effect  of  the  bomb  being  tossed.  Hence,  the 
weapon  release  point  can  be  some  distance  from  the  target;  in  this  case, 
the  aircraft  is  1400m  downrange.  The  tracking  data  for  this  pass  is  initiated 
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at  a  range  of  5180m — :4020m  downrange  and  3270m  in  elevation.  In  the  middle 
of  the  pass  there  is  a  sharp  3.4g  diving  turn  followed  by  ten  seconds  of 
rolling  back  and  forth  and  then  a  5.5g  maneuver  away  from  the  target.  Less 
than  25  seconds  elapse  between  the  initial  point  and  the  instant  of  minimum 
altitude  (560m) ,  which  occurs  just  after  the  aircraft  passes  the  target. 

Flight  path  #2  bears  some  similarity  to  pass  #1.  The  significant  differences 
are  that  the  rolling  back  and  forth  in  the  middle  of  the  pass  is  not  as 

pronounced  and  the  weapon  release  point  is  closer  to  the  target  since  this 

pass  is  representative  of  the  common  "dive"  maneuver. 

Flight  path  #3  is  called  a  "pop-up"  maneuver.  Here,  the  pass  is  initiated 
with  the  aircraft  5800m  downrange,  at  an  altitude  of  260m  and  pulling 
almost  6g's  as  it  starts  a  climb.  As  the  aircraft  climbs  over  the  next 
10  seconds  to  an  altitude  of  1440m,  it  rolls  over  so  that  at  the  peak  of 
its  climb  it  is  on  its  back.  It  continues  to  roll  as  it  begins  to  dive  at 

the  target.  It  performs  a  3g  turn  one  way  followed  quickly  by  a  4g  turn 

the  other  before  it  steadies  to  deliver  its  payload.  Then  it  executes  a 
6g  bank  away  from  the  target.  All  this  maneuvering  occurs  within  30  sec. 

To  compare  the  effectiveness  of  the  various  prediction  algorithms,  the  simulation 
program  discussed  in  Reference  2  is  used.  This  program  generates  a  two- 
dimensional  histogram  (prediction  times  vs  miss  distance) ,  which  will  be 
the  means  of  comparison.  The  ballistics  (i.e.  range  as  a  function  of  time) 
of  the  AA  rounds  are  assumed  to  be 

R  =  V  t/(l+0.129t)  (1) 

m 

where  V  (assumed  to  be  1175m/sec)  is  the  muzzle  velocity.  The  results 
presented  in  Reference  2  used  constant  velocity  ballistics) . 

PREDICTION  WITH  NOISELESS  OBSERVATIONS: 


Heretofore,  models  in  fire  control  predictors  have  been  rather  simple.  The 
current  field  peice  utilizes  a  linear  or  position-plus-rate-times-time 
algorithm.  Previous  studies  have  suggested  the  use  of  a  quadratic  algorithm  which 
is  the  above  plus  an  extra  term  comprised  of  acceleration  multiplied  by  one- 
half  of  the  prediction  time  squared.  A  variation  of  this  approach  to 
prediction  is  based  on  a  first-order  Markov  model  of  acceleration  in 
continuous  time ,  i.e., 

a(t)  -  -6/a(t)  +  u(t)  (2) 

Solving  this  general  equation  gives  the  following  discrete  predictors 

xn+k  =  Vl  +  TXn-l  +  1/u>2  (exP(’wT)  +6oT-1)  xn_i 
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(3) 


* 


A  *  2  #  * 

y  =  y  +  Ty  +  l/6o  (exp(-wT)  +  u/T-l)y 

n+k  n-1  n-1  n“l 

z  =  z  +  Tz  +  l/to2(exp(-WT)  +  CoT-1)  z 

n+k  n-1  n-1  n-1 


where  T  =  (k+l)4  =  (k+l)*0.1  second 


This  prediction  algorithm  will  be  used  as  the  benchmark  for  comparison  since 
it  demonstrated  a  prediction  capability  better  than  the  other  predictors  mentioned 
above.  The  results  obtained  using  this  predictor  in  the  simulation  program 
assuming  no  noise  in  the  observations  are  presented  in  Table  1.  Notice  that  for 
prediction  times  of  less  than  one  second  all  the  rounds  for  all  three  passes  score 
"hits”  (defined  as  a  miss  distance  less  than  5m) .  For  prediction  times  of  one 
to  two  seconds ,  the  fractions  of  hits  ranges  from  one- third  (passes  #1  and  #2) 
to  one-fifth  (pass  #3).  The  number  of  rounds  within  fifteen  meters  is  important 
too,  for  close  rounds  can  have  an  effect  on  the  pilot1 s  resolve  to  carry  out 
his  mission.  However,  Table  I  does  not  show  any  hits,  or  even  a  consistent 
number  of  close  rounds,  for  prediction  times  much  greater  than  two  seconds.  The 
similarity  of  the  effectiveness  of  this  algorithm  across  the  different  flight 
paths  should  be  noted. 

The  authors  in  Reference  2  proposed  that  acceleration-dot  can  be  modeled  as  a 
fifth-order  autoregressive  process  (see  difference  3  and  4),  i.e.. 


*„  -  ftV 1  +  ft  V  2  +  ftVs  *  M.-4  +  ft  Vs  * 


(4) 


where  the  residuals  (u  )  are  uncorrelated  and  zero-mean.  This  AR  model  can  be 
combined  with  the  stanHard  expansions,  (for  X-direction) , 


x  =  x  ,  +dv  +  (42/2)a  ,  +  (43/6)a  , 

n  n-1  n-i  n-1  n-1 

v  =  v  .  +Xla  +  (xi2/2)a 

n  n-1  n-1  n-1 


(5) 


a  =  a  +  da  ,  , 
n  n-1  n-1 


(where  &  =0.1  second),  to  yield  a  matrix  one-step  predictor 

A  B 


A 

s  = 
~n 


0  $ 


-  n-1 


(6) 


where  s  =  (x„  ,  V  an,  a^a^,  a^,  an_3,an_4)  • 


n 


n 


A  = 


1  &  &/2 

0  1  A 

0  0  1 
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r 


d3/6  0  0  0  0 

tf/ 2  o  0  0  0 

4  0  0  0  0 


and  £  is  the  zero  matrix.  Although  Akaike's  AIC  statistic  (5)  indicates  that  the 
optimal  model  order  is  larger  than  five,  it  was  found  that  models  with  lower  orders 
provide  a  significantly  poorer  fit  to  the  data  while  models  with  higher  orders  do 
not  give  substantial  improvements.  Therefore,  a  fifth-order  AR  model  is  used. 
Predictions  for  longer  times  of  flight  can  be  made  by  raising  the  partitioned 
matrix  in  Equation  6  to  the  appropriate  power.  For  position  prediction  there  are 
similar  matrices  for  the  Y  and  Z  directions. 


Implementation  of  this  prediction  algorithm  leads  to  the  results  in  Table  II.  The 
robustness  of  this  type  of  predictor  that  was  reported  in  Reference  2  can  be  seen 
in  these  results  since  a  single  set  of  coefficients 

h 
h 
b 

#4 

/*5 

was  used  for  all  three  directions  and  for  all  three  flight  paths  .  This  single 
set  was  derived  via  a  least-squares  criterion  from  the  X-direction  acceleration- 
dot  of  flight  path  #2.  Notice  that  now  almost  all  the  rounds  score  hits  for 
prediction  times  up  to  two  seconds,  twice  as  long  as  for  the  benchmark  predictor. 
Also,  there  is  a  reasonable  amount  of  success  for  prediction  times  of  two  or 
three  seconds  —  about  one-fifth  score  hits  for  passes  #2  and  #3  while  one-third 
doe  for  pass  #1,  plus  more  than  half  the  rounds  are  within  15m  for  all  three 
passes.  For  longer  prediction  times  the  results  differ  from  pass  to  pass. 

Pass  #1  still  has  a  good  number  of  close  rounds  while  pass  #3  has  almost  nothing. 
But  again,  for  prediction  times  up  to  three  seconds  at  least,  the  results  for 
this  algorithm  over  very  different  flight  profiles  demonstrate  a  consistent 
performance  which  is  a  great  improvement  over  the  benchmark  algorithm. 


=  4.029471 
=-6.536585 
=  5.313477 
=-2.148406 
=  0.340802 


* 
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PREDICTION  WITH  NOISE-CORRUPTED  OBSERVATIONS 


When  the  observations  are  corrupted  by  noise,  the  problems  become  two-fold;  first 
the  observations  must  be  filtered  to  give  estimates  of  the  various  quantities  at  the 
appropriate  time  instants,  and  second,  these  estimates  must  be  manipulated  to  yield 
reasonable  position  predictions.  Since  the  theory  says  that  a  Kalman  filter  is 
optimal  for  a  linear  system  (6)  and  since  the  assumption  that  acceleration-dot  is 
an  AR  process  allows  a  linear  definition  of  the  aircraft  motion,  a  Kalman  filter 
will  be  used  to  provide  the  state  estimates  in  each  direction. 

Given  a  model  of  the  aircraft* s  kinematic  motion  in  one  coordinate 

s^  =  Fs  -  +  u  ,  (7) 

-n  — n-l  -n 

where  the  model  residual  vectors  (un)  are  zero-mean  and  uncorrelated  with  covariance 
matrix  Q.  The  position  and  velocity  observations  are  given  by 

*=  Hs  +  v  ,  where  v,^  NfO  £  1  (8) 

— n  — n  -n  -n  U-  -nJ 


Then  the  Kalman  filter  can  be  written  as  follows: 


s  =  Fs  +  K  f  z  -  HFs  1 
n  —  n-l  “n  *-  -n  - n-ll 


where  the  Kalman  filter  gain  matrix  is 

K  =  P  H  '  (HP  h'  +  R  )_1 
-n  -n  —  — n~  -n 


(9) 


(10) 


and  the  error  covariance  matrix  is 


P 

-m 


*  E 


(I  -K  H) P  F 
—  “-n-l - n-i— 


+  £ 


(ID 


The  quantity  in  the  square  brackets  in  Equation  9  forms  the  innovations  process  which 
should  be  a  white  noise  process  if  the  filter  is  properly  tuned.  The  filtering 
algorithm  to  be  part  of  the  fire  control  simulation  program  will  include  three 
such  Kalman  filters  —  one  for  each  direction. 


The  sensors  employed  by  AA  weapons  systems  measure  ran<je  R,  the  elevation  angle  E, 
the  azimuth  angle  A,  range-rate  R  and  the  angle  rates  E  and  A.  Each  of  these 
measurements  will  have  some  independent  noise  in  it.  The  results  to  be  presented 
were  produced  assuming  gaussian  noises  with  the  following  one-sigma  values: 
range-2m,  elevation  and  azimuth  angles  -  0.5  milrad,  range  rate  -  2m/sec  and 
angle  rates  -  1%  of  the  actual  rate.  The  observations  z_  which  include  position 
and  velocity  in  the  X,  Y  and  Z  direction  of  a  particular  Kalman  filter  are  simply 
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transformations  of  the  noisy  (R,  E,  A,  R,  E,  A) f 


Since  the  work  in  this  study  was  performed  in  the  XYZ  coordinate  system  while 
the  noise  occurs  in  REA,  the  2x2  matrices  that  enter  into  the  three  XYZ  # 
Kalman  filters  must  be  derived  from  the  ocvariance  matrix  S_n  for  (Rr  E,  A,  R, 

E,  A) .  Since  the  noise  added  to  the  different  quantities  is  assumed  to  be 
uncorrelated,  this  matrix  will  be  zeros  everywhere  except  on  the  diagonal  which 
will  be  made  up  of  the  appropriate  variances.  Notice  that  since  the  angle  rate 
variances  are  a  function  of  the  data,  2.  (and  hence  the  R  's)  will  change 
with  time.  Now,  the  full  6x6  covariance  matrix  in  XYZ  is  given  by 

R  =  A  E  A  (12) 

— n  “*n  ”n 

where  A  is  the  Jacobian  matrix  of  the  transformation  at  time  tn.  In  order  to 
simplif?  the  calculations,  a  decoupling  of  the  XYZ  direction  filters  is  forced 
on  R  .  For  example,  the  R  matrix  for  X  has  the  variance  of  the  X-position 
noise  and  the  variance  of~?he  X-velocity  noise  on  the  diagonal  and  the  covariance 
between  the  two  on  the  off-diagonals.  The  covariances  between  directions  are 
lost.  Preliminary  work  indicates  that  filtering  in  REA  which  would  eliminate  the 
need  for  this  decoupling  will  provide  better  estimates.  However,  it  still 
seems  that  prediction  is  best  done  in  XYZ. 


This  study  will  present  results  obtained  using  two  different  filters.  The 
first  is  based  on  the  type  of  model  used  in  the  benchmark  predictor.  In  fact, 
a  first-order  AR  model  of  acceleration  with  the  coefficient  o< =  0.995  will  be 
used.  Tests  show  that  the  performance  of  this  filter  is  rather  insensitive  to  the 
exact  value  of  this  coefficient.  The  vectors  and  matrices  in  the  Kalman  filter 
formulation  defined  above  are  the  following  for  the  X-direction  case: 
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where  E  [u  21  =  q.  Only  position,  velocity  and  acceleration  are  being  estimated 

so  only  predictors  that  require  these  variables  can  be  used  with  this  filter. 

The  results  in  Table  III  were  produced  by  the  benchmark  predictor,  which  is  based 
on  a  first-order  Markov  model  of  acceleration,  operating  on  the  estimates 
produced  by  three  third-order  Kalman  filters.  The  results  in  these  histograms 
form  the  benchmark  in  this  scenario  of  noisy  measurements. 

The  second  filter  is  based  on  the  fifth-order  AR  model  of  acceleration-dot.  This 
means  that  the  state  transition  matrix  F  (now  an  8x8  matrix)  is  the  same  as  the  pre¬ 
dictor  matrix  of  Equation  6.  The  results  presented  here  were  generated  by 
three  Kalman  filters  —  a  different  set  of  coefficients  for  each  direction 
(produced  from  pass  #2) .  However,  the  same  filters  were  used  on  all  three 
flightpaths.  Equation  6  also  defines  the  state  vector  £.  The  8x8  residual  co- 
variance  matrix  £  wil1  be  a11  zeros  except  for  the  element  q(4,4)  which  is  equal  to  the 
variance  of  the  residuals  of  the  X,  Y  or  Z  AR  model  in  pass  #2.  Again  the 
observations  are  of  position  and  velocity  so  the  H_ matrix  is  now  2x8.  Since  this 
filter  produces  estimates  of  position  through  acceleration-dot,  the  benchmark 
and  higher  order  predictors  can  both  be  used. 

Figure  1  contains  the  spectral  density  estimates  of  the  velocity  innovations 
processes  for  the  X-direction  of  pass  #2  for  this  and  the  third-order  Kalman  filters. 

As  was  previously  stated,  if  a  filter  is  properly  tuned,  these  processes  should 
be  white.  It  is  clear  from  the  figure  that  the  higher  order  filter  has  a  flatter 
spectrum  which  indicates  that  it  is  more  closely  tuned  to  the  data  and  hence 
that  it  will  provide  better  estimates. 

Table  IV  contains  the  results  obtained  using  this  eigth-order  Kalman  filter  with 
the  benchmark  predictor.  A  large  increase  in  accuracy  for  short  prediction  times 
(less  than  one  second)  is  evident  in  all  three  flight  paths.  Also,  except  for 
pass  #2,  there  is  a  significant  over-all  improvement  in  effectiveness  for 
prediction  times  between  one  and  three  seconds.  It  seems  clear  that  given 
reasonable  noise  levels  little  can  be  expected  for  prediction  times  greater 
than  three  seconds. 

A  single  third-order  AR  model  of  acceleration-dot  with  coefficients 
P1  =  1.451057 
£  =-0.725526 

=  0.125C00 

provides  the  basis  for  a  position  predictor  that  was  used  against  all  directions 
and  flight  paths.  This  predictor  together  with  the  three-eigth  order  Kalman 
filters  used  above  gives  the  results  in  Table  V.  Note  there  is  a  further  increase 
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in  the  accuracy  for  short  prediction  times.  But  here,  for  the  longer 
prediction  times,  only  pass  #3  shows  any  significant  improvement  over  the 
results  in  Table  IV.  The  results  for  the  other  flight  paths  are  only 
minimally  better . 

PREDICTION  AGAINST  A  "WORSE  CASE  FLIGHT  PATH" 

One  can  associate  with  a  pt^  order  autoregression  process  a  characteristic 
equation  whose  roots  lie  inside  the  unit  circle  of  the  complex  plane  if  the 
process  is  stable.  Thus,  for  (4),  the  corresponding  characteristic  equation 
is  * 


JT 


A-0 


(14) 


When  the  roots  associated  with  (14)  for  each  of  the  flight  passes  are  plotted 
in  the  complex  plane,  it  is  seen  that  they  congregate  in  clusters.  One  can 
identify  five  clusters  per  coordinate.  Furthermore,  the  groupings  of  the 
clusters  is  remarkably  similar  in  each  coordinate.  This  is  suggestive  of 
a  robust  property  of  the  models  which  resulted  in  the  single  set  of  coefficients 
for  each  coordinate  exhibited  above. 

One  can  partition  each  cluster  and  design  a  fifth  order  A-R  process  which  has 
"worse  case"  properties  from  the  point  of  view  of  the  A-A  artillery  weapon. 

The  resulting  class  of  "worse  case  flight. paths"  formed  by  determining  the 
coefficients  from  the  root  locations  can  then  be  used  to  test  the 
effectiveness  of  the  fifth  order  predictors  based  on  real  flight  data. 

One  can  also  use  the  WCFP  predictor  to  see  how  well  it  performs  against  the 
real  aircraft  data. 

The  WCFP  was  designed  in  the  following  way:  One  root  per  cluster  was  chosen 
to  produce  a  set  of  yS  coefficients  such  that  the  steady  state  variance 
of  the  velocity  estimate  is  maximized.  The  condition  of  maximum  variance 
is  achieved  by  roots  located  on  the  boundary  of  each  cluster.  A  magnified 
view  of  two  clusters  and  the  WCFP  root  locations  is  shown  in  Figure  2.  A 
computer  search  over  each  boundary  determined  the  desired  roots.  The  ground 
track  of  a  WCFP  is  exhibited  in  Figure  3. 

Table  VI  shows  the  performance  of  three  predictors  against  the  WCFP.  Observe 
that  the  a-dot  predictor  with  the  WCFP  coefficients  achieves  the  highest 
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average  hit  probability ,  with  the  fifth  order  in  a-dot  with  x  coefficients 
from  flight  pass  10  a  close  second.  The  benchmark  predictbr  is  a  distant 
third. 

It  was  found  that  the  WCFP  predictor  also  performs  well  against  the  real  attack 
flight  paths. 

CONCLUSIONS :  , 

Higher  order  AR  Models  of  acceleration-dot  provide  robust  models  of  high 
performance  aircraft  maneuvers  which  are  much  improved  over  the  models 
based  on  a  first-order  model  of  acceleration.  The  filtering  capabilities 
of  a  fire  control  system  could  be  greatly  enhanced  by  the  incorporation  of 
the  higher  order  models.  Further  enhancement  can  be  achieved  by  the  use 
of  the  higher  order  model  in  the  predictor  as  well. 

These  improvements  are  noise  level  dependent/  of  course.  Higher  noise  levels 
mean  a  bigger  payoff  by  changing  to  the  higher  order  filter  with  less  of 
a  return  for  the  higher  order  predictors.  Lower  noise  levels  mean  the 
reverse  with  the  payoff  for  the  higher  order  predictors  approaching  the 
noiseless  case.  However,  in  the  case  of  lower  noise  levels,  the  higher 
order  filter  is  still  expected  to  yield  significantly  better  estimates  than 
the  lower  order  filter.  In  addition,  it  is  needed  to  provide  the  higher 
order  predictor  with  the  estimated  quantities  that  it  requires. 
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FICDEE  1.  SPECTRAL  DENSITY  ESTIMATES  OF  X-DIRECTION 
VELOCITY  INNOVATIONS  PROCESSES  FOR  PASS  ?2 
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SPECTRAL  DENSITY  ON  LOGARITHMIC  SCALE  VS  FREQUENCY  (HZ) 


TABLE  I.  BENCHMARK  PRF.DICTOR 
(NOISELESS) 


TABLE  II.  AR  A-DOT  PREDICTOR 
(NOISELESS) 
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TABLE  III.  THIRD-ORDER  FILTER  WITH 
BENCHMARK  PREDICTOR 


FLIGHT  PATH  1 


Pred.. 

Number  of  rounds 

tint 

w/  miss  between 

Mean 

between 

0-5  -10-15-20  T 

iiiss 

"  3. "37 

0  -  JT 

34  0  0  -  46 

1  -  2 

4  9  7  13  -  42 

14.43 

2-3 

1  4  3  9  -  36 

36.78 

3  -  4 

0  1  2  0  -  39 

60.63 

4-5 

0  0  0  0  -  29 

112.99 

Total-T97 

'40:75 

FLIGHT  PATH  2 

Pred. 

Humber  of  rounds 

time 

w /  miss  between 

Mean 

between 

0-5  -10-15-20  T 

Miss 

0-1 

urn  3  o  =“5? 

1-2 

9  15  10  2  -  41 

10.10 

2  -  3 

0  4  2  2  -  37 

60.34 

3-4 

0  0  0  0  -  31 

105.92 

4-5 

0  0  0  0  -  35 

146.63 

Total-!?! 

6d73£ 

FLIGHT  PATH  3 

Pred. 

Number  of  rounds 

time 

w/  miss  between 

Mean 

between 

0-5  -10-15-20  T 

Miss 

0~T“ 

27  17  0  0  -  44 

1  -  2 

3  10  28  4  -  45 

11.07 

2-3 

0  0  0  1  -  36 

57.72 

3-4 

0  0  0  0  -  39 

130.05 

4-5 

0  0  0  0  -  37 

212.25 

Total-201 

78.08 

TABLE  IV.  EIGHTH -dhDER  FILTER  WITH 
BENCHMARK  PREDICTOR 


Pred. 

tine 

between 

FLIGHT  PATH  1 
Number  of  rounds 
w /  miss  between 

0-5  -10-15-20  T 

Mean 

Miss 

T-  1  ' 

40 

b 

0  -  46 

7T4T 

1-2 

11 

6 

10 

15  -  42 

11.21 

2  -  3 

2 

5 

3 

7  -  36 

31.63 

3  -  4 

0 

1 

3 

1-39 

53.45 

4-5 

0 

0 

0^ 

0  -  33 

107.54 

T 

otal-ITF 

77757 

FLIGHT 

PATH  2 

Pred. 

Number 

of 

rounds 

time 

w/  miss 

between 

Mean 

between 

0-5 

-10-15-20  T 

Miss 

6-1 

35  14 

0-49 

Til 

1-2 

9  16 

12 

0-41 

9.45 

2-3 

0 

2 

5 

3-37 

53.16 

3-4 

0 

0 

0 

0-32 

36.33 

4-5 

0 

0 

0 

0-34 

150.54 

Total-177 

54.92 

FLIGHT 

PATH  3 

Pred. 

Humber 

of 

rounds 

tine 

w/ 

miss 

between 

Mean 

between 

0-5 

-10-15-20  T 

Miss 

0-1 

ITT" 

0-44 

2.72 

1-2 

7  : 

Z1 

17 

0-45 

8.57 

2-3 

0 

0 

3 

2-37 

50.16 

3-4 

0 

0 

0 

0-33 

122.41 

4-5 

0 

0 

0 

0-36 

209.39 

Total-ZOT 

77777 

TABLE  V.  EIGHTH-ORDER  FILTER  WITH 
THIRD-ORDER  A- DOT  PREDICTOR 


FLIGHT  PATH  1 


Pred. 

Number  of  rounds 

time 

w/  miss  between 

Mean 

between 

0-5  -10-15-20  T 

Miss 

T37 

0-1" 

43  1  0  0  -'46 

1-2 

11  8  14  9  -  42 

10.20 

2-3 

1  4  3  8  -  36 

29.10 

3-4 

0  0  1  1  -  40 

62.15 

4-5 

0  0  11  -  33 

87.04 

Total-197 

33774 

FLIGHT  PATH  2 

Pred. 

Number  of  rounds 

time 

w/  miss  between 

Mean 

between 

0-5  -10-15-20  T 

Miss 

“7779 

0  -  1 

50  9  0  0  =-57 

1-2 

10  16  11  2  -  41 

8.09 

2-3 

2  2  4  4  -  37 

56.46 

3-4 

0  0  1  0  -  32 

76.33 

4-5 

0  0  0  0  -  36 

153.82 

Total-175- 

54. 32 

FLICHT  PATH  3- 

Pred. 

Huaber  of  rounds 

time 

w/  miss  between 

Mean 

between 

0-5  -10-15-20  T 

Miss 

“7771 

0  -  I 

3?  5  0  0  -  44 

1-2 

14  13  13  0  -  45 

7.02 

2-3 

1  I  3  3  -  38 

43.84 

3-4 

0  0  0  0  -  38 

116.09 

4-5 

0  0  0  0  -  37 

232.39 

Total-207 

•75750 
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TABLE  iv 


WORST  CASE  FLIGHT  PATH  (NOISELESS) 


FIFTH-ORDER  AR  IN  A-DOT  PREDICTOR  (WCFP  COEFFICIENTS) 


Prediction  Number  of  Rounds  with  Closest  Miss 

Time  Approach  Miss  (in  meters)  between:  RMS 


between 

1 

in 

l 

o 

5-10 

10-15 

15-20 

20-25 

>25 

Total 

Error 

31 

0 

0 

0 

0 

0 

s: 

31 

0.01 

43 

0 

0 

0 

0 

0 

= 

43 

0.48 

2.0  -  3.0 

30 

8 

1 

0 

0 

0 

ss 

39 

4.88 

3.0  -  4.0 

2 

12 

8 

4 

5 

4 

= 

35 

16.40 

4.0  -  5.0 

1 

0 

0 

1 

3 

24 

s 

29 

64.27 

(seconds)  -  - - - -  - -  - 

Total  Rounds  =  177 

Average  Hit  Probability  *  0.61578?  Total  RMS  Error  =  27.116 


FIRST-ORDER  IN  A  PREDICTOR  (BENCHMARK) 


Prediction 

Time 

between 

Number 
Approach 
0-5  5-10 

of  Rounds  with  Closest 

Miss  (in  meters)  between: 
10-15  15-20  20-25  >25  Total 

Miss 

RMS 

Error 

0.0  -  1.0 

0 

0 

0 

0 

0 

= 

31 

2.69 

1.0  -  2.0 

16 

21 

6 

0 

0 

0 

SS 

43 

6.74 

2.0  -  3.0 

1 

2 

8 

4 

6 

19 

ss 

40 

24.68 

3.0  -  4.0 

4 

2 

1 

1 

1 

20 

=r 

29 

93.31 

4.0  -  5.0 

0 

0 

0 

0 

0 

35 

- 

35 

111.84 

(seconds) 

— 

— 

— — — - 

—  — 

—  — — 

— - 

— 

Total  Rounds  *  178 

Average  Hit  Probability  =  0.30626?  Total  RMS  Error  =  63.458 


FIFTH-ORDER  AR  IN  A-DOT  PREDICTOR  (PASS#10  X  COEFFICIENTS) 


Prediction 

Time 

between 

Number  of  Rounds  with  Closest 
Approach  Miss  (in  meters)  between: 

0-5  5-10  10-15  15-20  20-25  >25  Total 

Miss 

RMS 

•  Error 

0.0  -  1.0 

31 

0 

0 

0 

0 

0 

= 

31 

0.04 

1.0  -  2.0 

44 

0 

0 

0 

0 

0 

= 

44 

0.65 

2.0  -  3.0 

19 

8 

11 

0 

0 

0 

ss 

38 

7.80 

3.0  -  4.0 

0 

15 

6 

1 

1 

11 

s 

34 

30.49 

4.0  -  5.0 

0 

0 

0 

2 

3 

27 

= 

32 

71.85 

(seconds) 

•  * 

Total  Rounds 

ss 

179 

'  Average  Hit 

Probability 

=  0. 

55486?  Total 

RMS 

Error 

=  33.352 

Next  page  is  blank 
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A  DESIGN  METHODOLOGY  FOR  ESTIMATORS  AND  PREDICTORS 
IN  FIRE  CONTROL  SYSTEMS 
BY 

JAMES  F.  LEATHRUM 
CONSULTANT  TO  AMSAA 

ABSTRACT 


The  state-of-the-art  in  the  design  of  Kalman  filters  for  fire  control 
systems  leaves  the  designer  with  several  parameters  to  be  used  to  over¬ 
come  the  effects  of  modelling  errors.  The  fixing  of  these  parameters 
usually  requires  extensive  simulation  and  trial-and-error  searching  for 

satisfactory  operating  conditions.  In  the  process,  the  effects  of  the 
various  modelling  errors  are  easily  confounded  and  the  intuitive  under¬ 
standing  of  target  behavior  often  lost. 

The  purpose  of  the  analysis  reported  here  is  to  establish  a  design 
methodology  which  begins  with  the  allowable  variances  in  miss  distance 
and  leads  directly  to  filter  parameters  for  an  optimal  filter.  Structural 
mismatching  between  the  filter  and  the  actual  target  are  left  to  analysis 
by  simulation.  The  advantage  of  this  methodology  is  that  the  performance 
of  the  fire  control  system  in  terms  of  miss  distance  enters  the  design 
process  at  the  outset  rather  than  as  a  "take  what  we  get"  outcome.  Math¬ 
ematically,  the  process  involves  explicit  solution  of  the  steady  state 
matrix  Riccati  equation. 
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Conventional  Deslqn  Methodolo 


The  conventional  approach  to  the  design  of  estimators  and  predictors 
for  fire  control  systems  is  best  illustrated  by  the  following  development 
of  models  and  parameters.  One  would  start  by  formulating  target  and 
observer  models  of  the  form. 

(a)  Target  Model 

Xk+1  3  $kXk  +  BkUk 
(2)  Observer  Model 

vk  •  Vk +  \ 

These  models  immediately  involve  a  linearization  approximation.  The  target 
model  captures  the  well  defined  motion  in  the  state  transition  matrix,  <5^, 
and  leaves  the  less  defined  part  of  the  motion  to  a  noise  term,  B^U^.  The 
observer  is  usually  a  statement  that  not  all  the  state  components  are  visible, 
and  that  the  observations  are  corrupted  by  an  error,  V. .  (The  index,  k,  is  a 
discrete  time  index). 

If  one  can  further  approximate  and  by  white  gaussian,  zero  mean 
processes,  an  estimator  of  can  be  formulated  as: 


Xk+1  3  \  Xk 


Xk  3  xk  +  Kk  (Yk  ‘  Hk  V 


Which  is  the  Kalman  Filter  wherein 


Kk  *  pkHVRk +  V  H'k> 


p  3  Wk +  Wk 
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In  the  most  sophisticated  fire  control  systems,  the  target  noise  is 
represented  in  a  target  oriented  coordinate  system.  Thus,  the  given  Qk 
will  rotate  as  the  target  moves  which  in  turn  leads  to  a  nonsteady  K^. 

The  Kalman  gains  tend  to  change  throughout  the  estimation  process.  In 
addition,  may  be  range  dependent  which  leads  to  further  variability 
1"  Kk. 

In  designing  such  a  filter,  the  implementor  is  left  with  choices  of 
the  magnitude  of  Qk  and  (i.e.,HQkll  and  HRjH).  A  conventional  design 
process  would  require  assessing  j(  Rk/|  from  the  accuracy  of  the  instrumentation 
used  by  the  observer.  Since  llQ^/l  represents  unmodelled  behavior,  it  is 
usually  adjusted  to  achieve  some  other  objective,  such  as  white  innovation,  or 
minimum  ensemble  miss  distances.  Whatever  the  objective,  the  last  phase  is 
unguided  by  the  theory  and  thus  usually  requires  extensive  simulation.  The 
design  process  for  a  filter-predictor  in  tandom  is  illustrated  by 


A  Direct  Design  Methodology 

A  methodology  which  could  utilize  maximum  allowable  miss  distances 
to  assess  the  design  parameters  directly  would  have  some  obvious  advantages 
over  the  conventional  process.  Such  a  methodology  is  proposed  here  with  the 
following  features. 


Allowable 
Miss  Distance 

i 

(i.e.  lc  miss) 


Power  Spectrum  of  the  target  motion 
_ _ 1 _ 


Filter 

||  Qjl  ,  II  Rk/I 

Design 

► 

Predictor  Algorithm 

The  principles  of  the  design  will  be  illustrated  by  restricting  the 
discussion  to  a  single  dimension  and  further  restricing  models  to 
$  :  Upper  Triangular 

Q^:  Scalar  constant,  q 
Rk:  Scalar  constant,  r 

These  restrictions  do  not  limit  variability  of  the  gains  in  the  final 
design,  but  only  allow  one  to  focus  attention  on  the  magnitudes  of  the 
parameters  in  each  direction.  The  design  process  requires  solution  of 
the  steady  state  filter  operations  which  become 


P  =  $P4  +  BqB'  -$RT(R  +  HPH'fV* 

'Ir 

The  solution  for  P  In  terms  of  r  and  q  requires  iteration.  However 
a  closed  form  solution  for  P/q  and  r/q  in  terms  of  band  width  of  the  target 
motion  is  possible.  It  requires  the  observation  from  the  analogous 
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continuous  models  that 


% 


P33/q  =  (^/l)  It 


Where  wv  is  the  bandwidth  for  velocity. 

(a)  is  the  bandwidth  for  acceleration 

a 

At  s  t^  -  tk;  the  time  increment  between  observations 
The  required  bandwidths  can  be  assessed  from  the  power  spectrum  observed 
in  field  tests  of  generic  targets.  The  ratios  (P/q  and  r/q)  completely 
specify  the  one  dimensional,  steady  state  design,  but  they  do  not  product 
the  parameters  needed  for  |{  Q^Hand  ((  R^lj  in  a  multidimensional  design.  The 

required  magnitudes  are  obtained  from  the  variance  of  miss  distance 

r 


amiss 


f  rq 


-f>n/q  +\fP22/q,tf4'fP 


33/ q 


tJ/2 


Where  in  an  optimal  design,  the  P  is  interpreted  as  the  variance  of  the 
estimator  error.  The  above  equation  represents  the  variance  propagation 
through  a  second  order  predictor  (t^  is  the  time  of  flight.)  Other 
predictor  algorithms  could  be  used  at  this  point.  Since  the  ratio  P/  , 
is  computed,  by  asserting  a  la  miss,  one  can  directly  determine  q.  Further, 
knowledge  of  r/q  leads  directly  to  r. 

The  logical  outcome  of  this  process  is  the  question  of  whether  an 
observer  with  an  accuracy  on  the  order  of  fr  is  achievable.  The  power 
spectrum  of  the  target  motion,  and  the  limits  of  miss  distance  in 
conjunction  may  force  the  enhancement  of  instrumentation  technology. 

A  Typical  Design 

The  purpose  of  this  section  is  to  illustrate  the  design  process 
by  starting  with 

1)  A  Second  order  predictor. 
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2)  Bandwidths  of  the  target  motion 

wv  =  0.128hz 
=  0.160hz 

3)  1  amiss  distance  =  1  meter. 

4)  Time  offlight  *  2.0  sec. 

In  the  order  which  they  would  be  computed,  the  variance  ratios  are 


i »  J 

PU/q 

3.3 

.0737 

2,3 

.02716 

1,2 

.00367 

1,3 

.00534 

2,2 

.0146 

1,1 

.00126 

and  then  r/q  =  1.59X10"3 

In  this  computation,  the  model  coeffi cents  are 


1. 


$  “ 


H  = 


1. 


10"3/6 

.005 

0.1 


where  the  $  and  B  parts  of  the  model  are  determined  by  the  data  rate  of 
At  *  0.1  sec.  From  the  la  miss  distance,  the  q  is  fixed  at 


q  »  1.88 

thence  =  2.99  X  lo~3 

Given  the  parameters  of  this  example,  the  technological  conclusion  is 
that  an  observer  with  a  la  accuracy  of  r  =  0.0547  meters  is  needed  to 
achieve  a  1  a  miss  distance  of  1  meter. 

The  observer  accuracy  vs.  miss  distance  is  summarized  in  the  following 
table. 
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la  Miss  Distance,  meters 

q 

fr,  meters 

1.0 

1.88 

0.0547 

1.5 

4.23 

0.0820 

2.0 

7.52 

0.1094 

2.5 

11.75 

0.1367 

3.0 

16.92 

0.1641 

Conclusions 

The  design  process  is  completely  characterized  by  assessing  the 
proportionality  constant  in 

r  =  C  .  amiss 

where  C  is  determined  by  the  bandwidths  of  the  target  motion  and  by  the 
predictor  algorithm.  The  computation  of  C  may  be  tedious,  but  it  is 
straight  forward  and  free  of  iteration. 

Although  the  impact  of  this  methodology  is  clear  from  the  models 
used  here,  its  full  generalization  is  yet  to  be  worked  out.  More  general 
forms  of  $,  B,  and  H  need  to  be  considered  in  the  interest  of  establishing 
the  theoretical  limits  of  the  design  process. 
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AN  ADAPTIVE  LEAD  PREDICTION  ALGORITHM  FOR  MAN I^C:M^AY31 

Pak  T.  Yip  &  Norman  P.  Coleman 
USA  ARRADCOM 
Dover,  NJ  07801 

ABSTRACT .  An  algorithm  concept  which  processes  target  bearing 
and  range  input  data  and  provides  "optimal”  estimates  of  target  position, 
velocity  and  acceleration  a  time-of-f light  in  the  future  is  discussed. 

Since  the  algorithm  concept  involves  certain  important  statistical  assump¬ 
tions  about  target  acceleration  dynamic  models,  these  assumptions  will  be 
discussed  in  detail  along  with  several  important  methods  used  in  the  model 
identification  process.  Secondly,  the  filter  algorithm  itself  will  be  dis¬ 
cussed.  This  algorithm  involves  the  parallel  processing  of  target  range  and 
bearing  data  by  several  extended  Kalman  Filters  corresponding  to  distinct 
maneuver  characteristics  of  anticipated  target  vehicles.  At  time  of.  fire 
the  filter  with  the  largest  computed  likelihood  function  is  selected  for 
lead  prediction.  Finally,  results  of  simulation  studies  in  which  actual 
target  path  data  is  used  to  generate  filter  input  data  for  hit  probability 
evaluation  is  discussed.  Comparisons  are  made  between  the  adaptive  al¬ 
gorithm  and  non-adaptive  first  order  algorithms. 

I.  INTRODUCTION.  This  paper  describes  a  multiple  model  adaptive 
Kalman  Filter  approach  to  the  problem  of  estimating  and  predicting  the 
position,  velocity  and  acceleration  states  of  tank  targets  of  varying 
maneuverability.  The  estimation  and  prediction  problem  presupposes  that  the 
range  and  angle  DATA  (measurements  corrupted  by  Gaussian  white  noise)  is 
available.  The  target  dynamics  is  described  by  a  system  equation.  Our 
solution  to  this  problem  is  an  adaptive  algorithm  implementable  in  real  time 
with  a  microprocessor  to  compute  target  position  a  projectile  time  of  flight 
in  the  future.  This  study  begins  with  the  selection  of  the  Antitank  Miss¬ 
ile  Test  (ATMT)  Phase  II  data}  to  identify  the  filter  acceleration  models. 
It  consists  of  three  dimensional  (x,y,z)  position  data  recorded  at  approx¬ 
imately  10  samples  per  second.  Maximum  likelihood  identification  method  is 
applied  to  this  data  to  identify  a  finite  set  of  Markov  Acceleration  Models 
which  are  representative  of  a  broad  spectrum  of  vehicle  maneuvers  consider¬ 
ed  likely  to  occur  in  actual  engagements.  These  models  provide  the  requir¬ 
ed  state  variable  description  of  the  target  dynamics  used  in  the  formulation 
of  the  multiple  model  extended  Kalman  Filter  Algorithm  for  lead  prediction. 
The  extended  Kalman  Filter  is  required  in  this  application  as  a  result 
of  nonlinearities  induced  by  target  coordinate  transformations  and  non¬ 
linear  measurement  equation. 

The  adaptive  lead  prediction  concept  is  based  on  the  simultaneous 
(parallel)  processing  of  the  discrete  extended  Kalman  Filters  corresponding 
to  the  distinct  target  models  identified  from  the  ATMT  data.  The  likeli¬ 
hood  function  associated  with  each  filter  is  computed  up  to  the  time  of 
fire  of  the  weapon,  and  the  filter  having  the  greatest  likelihood  is 
automatically  selected  for  lead  prediction. 
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In  the  present  study,  only  the  azimuth  and  range  information  of  the  tar¬ 
get  is  processed  in  the  filter  with  the  target  elevation  considered  con¬ 
stant.  The  performance  of  this  design  is  examined  with  a  Monte  Carlo 
simulation  and  the  sensitivity  of  the  lead  estimates  to  measurement  noise, 
level  of  target  maneuver,  range  sampling  rate,  and  time  of  flight  of  pro¬ 
jectile  are  analyzed  to  determine  the  feasibility  of  using  this  algorithm 
for  fire  control  lead  prediction  against  various  maneuvering  targets. 


II.  DATA  ANALYSIS.  The  ATMT  data  consists  of  six  tracks  pro¬ 
duced  by  a  M60AI  tank,  a  Scout  Vehicle  and  a  Twister  Vehicle  undergoing 
evasive  maneuvers.  The  M60A1  tank  is  capable  of  speeds  of  10  to  16  miles 
per  hour  and  with  a  maximum'  acceleration  of  approximately  .3g.  The  Scout 
is  an  armored  reconnassiance  vehicle  capable  of  moving  at  a  speed  of 
15  to  25  miles  per  hour  and  a  maximum  acceleration  of  approximately  . 5g. 
Since  our  only  interest  is  in  modeling  the  acceleration,  the  position  data 
is  sampled  at  a  frequency  of  2  cps  and  twice  differentiated  to  obtain  the 
acceleration  estimates  which  are  then  resolved  into  along-track  and  cross¬ 
track  components.  The  power  spectral  density  of  this  data  is  computed  by 
the  maximum  entropy  method3  which  assumes  the  data  is  generated  by  an  auto¬ 
regressive  process.  The  power  spectral  desity  S (f )  is  given  by 


S(f) 


_ 2  2* _ 

i  -  f 

C  =  1 


where  C\  is  the  standard  deviation  of  a  Gaussian  noise  process ;  o<£ is  the  i-th 
coefficient  of  the  autoregressive  process;  M  is  the  number  of  coefficients, 
and  the  coef f icients  ©<£ are  estimated  recursively3. 

The  number  of  the  autoregressive  coefficients  is  usually  larger  than  3 
which  is  not  desirable  for  Kalman  Filtering.  However,  the  power  density 
spectrum  affords  enough  information  for  estimating  essential  poles  and 
zeros  of  a  simplier  model  structure.  Later  the  maximum  likelihood  identif¬ 
ication  program  is  used  to  fine  tune  the  pole  and  zero  estimates. 


The  simplified  model  determined  from  the  spectral  analysis  has  the  follow¬ 
ing  form: 

40- _ -  -o 


where  q(s)  is  the  Gaussian  noise  process;  A(c)  is  the  system  acceleration; 

/Qt  and  are  parameters  to  be  identified  for  the  chosen  tracks  and 
each  of  the  along-track  and  cross-track  formulations. 

III.  DISCRETE  EXTENDED  KALMAN  FILTER.  The  system  and  the  measurement 
equations  are  readily  defined  as  follows: 

xk  =  i(ik-i  >  Ji)  +  i* 

-  b  ( xK)  +*k 
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where  Xk  system  state  vector  at  the  discrete  time  kdt  in  the  cartesian 

coordinate  system,  $  the  system  function  containing  all  information  about 
the  system  dynamics,  tj^the  plant ‘noise  vector,  15* the  measurement  vector, 
h(Xk)  a  vector  containing  the  true  range  and  azimuth  angle  of  the  target 
position  at  time  kdt,  r^  the  measurement  noise  vector,  and  dt  the  time 
between  two  samples.  ..  ... 

The  necessary  statistics  and  conditions  are: 

cov(ii  >5 j)  =  Qi  fij 

cov(r;  >r j)  =  Ri  sii 

cov(iitrj)  =o.-  . ;  Vi, j 

E(x.)  =  2. 

cov(£„)  =  p„ 

where  £ij  is  the  Kronecker  Delta. 


Given  the  above,  the  discrete  Extended  Kalman  Filter  equations  can  be 
written  as  follows:  The  predicted  state  estimate  vector  is  given  by 

><k*r|k  =  £(&<Ji) 

and  the  state  error  a  priori  covariance  matrix  by 

Pk+l|k  **  $  Pk  $T  + 


where 


a ±(2,  Jt) 

®  "  - ST” 


x  =xk 


$  =  \  +  *kdt  +  5?  (d-t)/2 


/«>  a 

— rrXk 


Xk  “  Xk 


The  updated  state  estimate  vector  can  be  written  as  follows: 

^k+l  88  ^k+ljk  K  2k  +  | 

~  2k  +  i  ~  h  (  Xk+ljk) 


where 


k+l 


K  -  Pk+I|k  HT(  H  Pk+,|k  Hr  +  Rk+i) 

m  _  Shte) 

"  t  ”^k  +  lfk 


-I 


h(w)  ={C(x,)1  +0<ifj 


aV/4 


tan 


-i 


(X,/Xa)} 
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x^,  represents  x,y  position  state  estimates  respectively  in  fixed  cartesian 
coordinates.  The  state  error  a  posteriori  covariance  matrix  is  given  by 


K H  ~~  fk  +  lfk  k  H  + 

and 

Qk  =  f  "  i(tk-t)Qs  iT(tk-*)  di 

V, 

where  the  continuous  case  plant  noise  covariance  matrix,  Qs,  is  known. 

The  continuous  time  system  dynamic  equations  used  in  deriving  the  dis¬ 
crete  time  equations  are  given  by 

x,  =x3  .  x2=x4 

X3  =  (X3  A&  +  X4  Ac)/  V 

X4  =  ( X4  Aa  ~  ^zAc)/ V 
Xs  =  "T^ai  X5  — fizLiXi  » 

X7  —  “7®ci  X7  ~7®c2  Xg  > 

S2  +/S^i  S  +  /Saa 

•s  +  yc 
Sz  +  ^6C|  S  + 

(xi  +  xS)'/4 


a4  = 

Ac  *“ 
V  = 


X6  *  X* 
Xg  —  X7 


where  X3  and  are  the  corresponding  X  and  Y  components  of  the  velocity 
vector;  Aa  is  the  target  acceleration  along  the  velocity  vector;  Ac  is  the 
target  acceleration  perpendicular  to  the  velocity  vector. 

With  this  filter,  target  range  and  angle  measurements  may  be  processed  to 
generate  target  state  estimate  recursively.  Before  defining  an  adaptive 
filter  procedure,  the  parameters  of  the  Markov  model  need  to  be  identified. 


IV.  LIKELIHOOD  FUNCTION  &  MAXIMUM  LIKELIHOOD  IDENTIFICATION  OF 

PARAMETERS 9 .  Given  a  parameter  vector  ox  ,  the  probability  of 
occurence  of  the  measurement  vector  sequence  Z^can  be  represented  by  a  multi¬ 
variate  gaussian  distribution. 


f  (zk ;  *) 

t(s  k|?K"'* 


? dais'  ;<) 

a _  exp^-TjU^K) 


—  H  PkJk-iH  +  Rk 


85 


where 


P(2l^;£)  -..the  likelihood  function 

n  =  number  of  elements  in  the  measurement  vector  2^. 

k 

In  order  to  identify  the  best  parameter  vector  c*  to  give  a  maximum  p(Z  ;o<), 
we  can  equivalently  minimize  the  negative  log  likelihood  function: 

M  «)  =  h  {i  S7  sr'  %  +  i  Ik  (  j  e 4  Si)} 

Since  the  term  (2.7T  )n^  in  the  likelihood  function  does  not  contribute  any 
interesting  information  it  has  been  eliminated  in  forming  M(.Z*;£j!)  .  The 
Gauss -Newton  method  is  used  in  the  minimization  procedure. 


sSj+i  -Kj  -  f>  D 


-I  b  M  (zk;^i) 

^  —.1* 


where 


1-for  this  method, 


D  =  .E 


( 


and  D,  the  expected  Hessian 

} 


The  test  for  convergence  is  given  by 

(sSj+i  j)r  D  (—j+i  - 


V.  PARALLEL  FILTERS  &  ADAPTIVE  ESTIMATION.  ’  Target  state  prediction 
for  maneuvering  ground  targets  have  never  been  a  simple  task  to  undertake. 

The  major  uncertainty  comes  from  the  target  driver’s  (stochastic)  decision 
to  maneuver.  However,  it  appears  there  exists  a  maximum  level  of  maneuver 
that  the  ground  vehicles  studied  can  attain.  This  maximum  level  provides  a 
non-trivial  range  of  dynamic  motion  that  can  be  quantized  to  a  finite  number 
of  maneuver  levels.  In  this  study,  five  filters  are  incorporated  into  the 
multiple  model  filter.  Model  Ml  (Filter  1)  is  a  simple  4  states  constant 
velocity  filter.  The  remaining  4  filters  are  identified  with  various  maneuver 
levels. 


The  adaptive  estimation  is  a  straight  forward  decision  making  process. 
Measurement  in’ range  and  azimuth  angle  are  processed  through  the  parallel 
filters.  The  filter  having  the  largest  likelihood  function  is  automatically 
chosen  to  provide  the  best  estimate  for  lead  prediction  and  gun  orders.  Two 
concepts  of  adaptive  prediction  are  examined.  In  concept  A  the  likelihood 
functions  account  for  the  entire  measurement  history  up  to  the  time  of  fire. 
Thus  this  adaptive  prediction  concept  is  good  against  targets  with  constant 
maneuver  level.  In  concept  B,  only  .the  last  ten  samples  prior  to  the  firing 
time  are  used  to  compute  the  likelihood  functions.  This  adaptive  filter 
concept  tends  to  be  more  sensitive  to  changes  in  target  maneuver  levels. 

VI.  SIMULATION.  A  Monte  Carlo  simulation  of  100  runs  was  set  up 
to  process  a  number  of  10  second  segments  from  the  ATMT  data  representing 
various  maneuver  levels  for  the  M60,  Twister  and  Scout  Vehicles.  These 
segments  of  data  are  different  from  those  used  for  the  parameter  identification 
tasks  discussed  earlier. 
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For  evaluating  the  system  performance,  the  perpendicular  miss  distance  of 
the  predicted  line  of  sight  from  the  real  target  position  is  defined  as  the 
prediction  error  Ep  in  meters.  The  firing  time  points  are  fixed  for  each 
segment  under  process.  The  performance  indicator  ph  at  each  firing  time 
point  is  defined  as  the  ratio  of  the  number  of  times  that  the  prediction 
error  Ep  is  less  than  1.15  meters' to  the  total  number  of  runs.  Actually,  they 
are  hit  probabilities  considering  the  prediction  errors  alone. 

Assuming  engagement  range  of  approximately  2000  meters,  45°  cross  range 
(across  the  range  vector),  1  cr  range  measurement  error  of  2  meters,  1  cr 
azimuth  tracking  error  of  0.3  mils,  a  projectile  speed  of  1500  meters  per 
second  and  using  the  adaptive  prediction  concept  A,  the  hit  probability  re¬ 
sults  are  illustrated  in  Figure!  and  summarized  in  the  following  table: 


BH I 

Number  of  Cases, 
7  Firing  Points 
per  Case _ 

Nean 

Const.  Uelocity 
Prediction _ 

ph 

Adaptive 

Prediction 

M60A1 

13 

.41 

.49 

Scout 

10 

.27 

.38 

Twister 

8 

.20 

.26 

For  an  engagement  range  of  approximately  1158,  60°  cross  range,  1  cr 
range  measurement  error  of  3  meters,  1  cr  azimuth  tracking  error  of  0.3  mils, 
a  projectile  speed  of  1158  meters  per  second  and  using  the  adaptive  prediction 
concept  B,  the  hit  probability  results  are  summarized  in  the  following  table: 


Target 

Number  of  Cases, 

Nean 

ph 

Type 

7  Firing  Points 

Const.  Uelocity 

Adaptive 

oer  Case 

_ Prediction. 

Prediction 

M60A1 

6 

.51 

.56 

Twister 

6 

.31 

.37 

With  the  latter  conditions,  the  sensitivities  of  the  system  are  observed 
for  a  particular  maneuvering  segment  as  shown  in  Figure  2.  Figure  3 
illustrates  the  system  range  (hence  the  time  of  flight  of  projectile)  sen¬ 
sitivity.  Figure  4  illustrates  the  system  sensitivity  to  angular  measurement 
noise.  Figure  5  illustrates  the  system  sensitivity  to  range  measurement 
noise.  Figure  6  illustrates  the  system  sensitivity  to  range  sampling  rate. 

VII.  DISCUSSION  &  FUTURE  PLAN.  This  study  has  demonstrated  that 
maneuvering  target  acceleration  may  be  adequately  modeled  as  a  discrete 
set  of  stationary  Markov  processes  whose  parameters  can  be  identified  off 
line.  Parallel  discrete  extended  Kalman  filters  have  been  used  to  success¬ 
fully  process  range  and  angle  measurements.  The  adaptive  selection  of  the 
most  appropriate  filter  at  each  time  step,  based  on  its  largest  likelihood 
function,  has  been  accomplished  on  line.  Representatives  maneuver  patterns 
and  levels  used  in  this  study  were  taken  from  the  ATMT  data  base.  The 
results  from  the  Monte  Carlo  simulations  indicate  that  the  performance  of 
the  multiple  model  adaptive  filter  design  is  generally  comparable  to  a 
filter  which  is  tuned  to  the  target  dynamics  of  that  particular  tracking 
interval.  In  particular,  the  results  show  that  the  adaptive  prediction 
consistently  performed  better  than  the  constant  velocity  prediction  with  an 
improvement  in  prediction  ranging  from  10  to  40  percent. 

Since  the  range  data  is  currently  not  a  uniformly  accessable  measurement, 
the  range  sampling  rate  has  been  examined  as  an  area  of  uncertainty 
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together  with  range,  angular  measurement  noise,  and  range  measurement  noise. 
The  results  indicate  that  the  system  performance  for. the  azimuth  channel  is 
heavily  dependent  of  angular  measurement  noise  and  projectile  time  of 
flight  in  terms  of  range,  and  is  not  very  sensitive  to  range  measurement 
noise  and  range  sampling  rate.  The  results  also  indicate  that  higher 
probability  of  hit  can  be  obtained  in  the. cross  range  geometry  than  in 
the  down  range  (coming  down  along  the  range  vector)  geometry. 

Implementation  of  this  filter ; algorithm  in  real  time  with  a  state  of 
the  art  microprocessor  is  in  the  planning  stage.  We  have  noticed  that 
Bierman1 s  UD  factorization^ for  the  state  error  covariance  propagation  is  a 
desirable  feature  considering  computation  accuracy  and  stability.  Several 
variations  of  the  existing  filter  algorithm  are  also  under  consideration. 
Finally,  a  complete  real  time  simulation  of  the  fire  control  system  with 
the  auto-tracker  or  human  operator  in  the  loop  and  filter  modifications 
to  improve  maneuver  detection  will  be  subjects  of  our  future  work. 
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APPLICATIONS  OF  DELAY  FEEDBACK  IN  CONTROL  SYSTEMS  DESIGN 

N.  P.  Coleman,  E. . Carroll,- D.  Lee  and  K.  Lee 
US  ARRADCOM 
Dover,  NJ  07801 


ABSTRACT:  Necessary  and  sufficient  conditions  for  exact  state  recon¬ 

struction  using  delays  are  discussed  together  with  an  example  in  which  the 
technique  is  implemented  in  real  time  using  an  8080/8085  microprocessor. 

Also,  a  frequency  domain  technique  for  synthesizing  certain  feedback  control 
laws  with  delays  is  developed  and  several  examples  discussed. 

I.  INTRODUCTION:  In  designing  a  control  system  using  optimal  control 
theory  or  classical  frequency  domain  techniques,  one  often  encounters  sit¬ 
uations  in  which  certain  required  signals  or  states  of  the  system  are 
unavailable  by  direct  measurement.  In  modern  control  design  this  problem  is 
usually  handled  by  implementing  some  form  of  reduced  order  or  full  order 
observer  which  provides  an  asymtotic  estimate  of  the  unmeasured  state.  In 

this  paper  a  technique  is  developed  for  exact  state  reconstruction  of  unmeasured 
system  states  using  values  of  the  measured  variables,  their  delayed  values 
and  the  control  variables  on  the  maximum  delay  interval .  Several  examples 
are  discussed  which  demonstrate  the  application  of  this  technique  on  a 
laboratory  servo  system  using  an  8080  microprocessor. 

A  second  application  of  delay  feedback  for  frequency  domain  compensation 
is  also  discussed.  A  frequency  domain  technique  is  developed  for  selecting 
appropriate  gain  and  delay  parameters  for  synthesizing  a  feedback  controller 
using  delays  in  the  output  and  several  applications  as  discussed. 

II.  REAL  TIME  STATE  RECONSTRUCTION  USING  DELAYS:  In  this  section  a 
technique  is  presented  for  exact  state  reconstruction  using  delay  feedback 
of  measured  states  of  a  control  system  and  the  values  of  the  control  input 
over  the  delay  interval.  A  real  time  application  of  this  technique  in  a 
servo  control  system  using  an  8080  microprocessor  is  also  discussed.  For 
simplicity,  consider  the  linear  time  invariant  system: 

x ( t )  =  Ax(t)  +  Bu(t)  (1) 

where  x  is  an  nxl  state  vector,  u  is  an  rxl  control  vector,  A  is  an  nxn 
constant  matrix,  and  B  is  an  nxr  constant  matrix.  Let  the  observation  vector 
y (t)  be  given  by: 

y (t)  =  Hx(t) 

where  y  is  a  mxl  vector,  and  H  is  an  mxn  constant  matrix.  Let  o  <  ^2 
be  time  delays. 

The  problem  is  to  reconstruct  the  state  x(t)  from  the  measurements  y(t), 
y(t-h),  •  *  • ,  y(t-h^)  and  the  measureable  control  vector  u(s),  t-h^s.s<Lt. 

The  following  argument  due  to  D.  H.  Chyung,  Reference  (  I  )  provides  the 
basis  for  a  real  time  state  reconstruction  algorithm  discussed  in  the  examples. 
This  argument  makes  use  of  the  well  known  variation  of  parameter  expression 
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for  the  time  response  x(t)  of  the  system  (1)  given  by 


x(t)  =  eA^  hi^x(t»hi)  +  jeA<^t  s^Bu(s) 


=  eAhix(t-hi)  +Je‘AsBu(t+s)ds 

-hi  i  =  1,2, -‘A 

Multiplying  both  sides  of  equation  (2)  by  He-A^i  results  in  the  equation 

;  ,° 

He~Ahlx(t)  =  Hx(t-h±)  +  He“Ahi j e“AsBu(t+s)ds 

-hi 


=  y(t-hi)  +  He  ^i  j  e~AsBu(t+s)ds 


i  =  1,2, ---A 


in  which  the  right  hand  side  is  completely  known.  Letting  C  denote  the  matrix 
given  by: 

Hje-AhTI 


C  =  He 


fe-AhX  • 


we  can  now  write  equation  (3)  in  the  form 
Cx(t)  =  Z(t) 


where; 


y(t~hi)  +  He^l  J^e~AsBu(t+s)ds 


z(t)  y(t-h2)  +  He  ^2  J^e  AsBu(t-fs)ds 

J  ~h2 

y(t-hj()  +  H e“AsBu(t+s)ds 


is  a  known  mXxl  vector  and  C  is  an  mj(xn  constant  matrix  depending  on  the 
parameters  hp  h2,  hj|.  If  the  matrix  C  has  rank  n,  then  equation  (4) 
can  be  written  as: 


x(t)  =  CTcj  -1CTz (t) 


where  C1  denotes  matrix  transpose. 
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Equation  (5)  has  several  important  implications;  First,  if  the  matrix  JC  CJ 
exists,  then  the  state  x(t)  can  be  exactly  reconstructed  from  the  measurement 
y(t),  its  delayed  values  and  the  input  signal  u(t),  o<  tsh^  ;  secondly,  the 
C  matrix  depends  only  on  the  delays  h^,  h^  so  the  right  hand  matrix  cal¬ 

culation  can  be  performed  completely  off  line.  This  leaves  only  the  relatively 
straight  forward  calculation  of  x(t)  and  a  matrix  multiplication  for  on-line 
microprocessor  computation.  This  latter  comment  is  of  particular  importance 
in  real  time  control  applications  in  which  relatively  low  speed  microprocessors 
are  utilized  for  control  law  implementation.  The  following  result  establishes 
the  condition  under  which  the  matrix  C  has  rank  n. 


Result:  There  exists  a  set  of  n  delays  o^h\-<h2  ^  a,  for  any  a>o 

such  that  the  matrix  C  has  rank  n,  if  and  only  if  rank  (Q)  =  n,  where 


Q 


HA 

0 


HAn_1 


Proof i  Let  a>o  and  assume  rank  (Q)  =  n.  Then  the  row  vectors  of  the  matrices 
He-Ahhg- [£ , a|  contains  n  independent  vectors  since,  if  not,  there  exists  b  £•  R 
such  that  He-Ahb=0  for  all  he  fo.al  .  Repeated  dif ferentation  with  respect  to 
h  gives  He_Ahb=HA-Ahb=HAn-1e  “®b=0.  This  implies  that  the  non  zero  vector 
e-Ahb  is  in  the  null  space  of  the  matrix  Q  and  hence  rank  (Q)<n. 


Conversely,  assume  rank  C  =  n,  then  rank  (Q)  -  n  since,  if  not,  exists 

b  +  o£Rn  such  that  Hb  =  HAb  =  •••  =  HAn_1b  =  o.  This  implies  He  ^b  =  o  for 
all  h  and  hence  rank  C<n,  a  contradiction. 


Example:  Evaluation  of  the  state  reconstruction  technique  given  by  equation  (5) 

was  carried  out  on  an  8080  microprocessor  development  system  which  was  in  turn 
interfaced  with  a  laboratory  servo  system  as  shown  in  Figure  1.  In  this  example 


the  system  state  vector  is  given  by 


r*f 

=  J2. 


where-  xj  is  the  motor  shaft  output 


The  measured  signal  is  xi  and  X2  is 
software  was  developed  and  debugged 


position  and  x2  is  the  motor1  shaft  velocity 
reconstructed  using  equation  (5) .  Once  the 

the  program  was  down-loaded  to  a  single  board  8085  microprocessor  shown  in 
Figure  2,  for  faster  execution.*  The  block  diagram  of  the  servo  system  without 
tack  feedback  is  shown  in  Figure  3. 


2406 

S+4 


-x2- 


Figure  3 


*  The  8085  configuration  shown  in  Figure  3  is  currently  being  used  to  evaluate 
digital  control  concepts  for  the  XM97  turret  system  shown  in  Figure  21. 
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The  state  space  equation  is  given  by: 


xi  0 


x2  -2406  .-4 


y  =  x 


The  state  transition  matrix  for  this  system  is  readily  computed  to  b 


-  e-2t 


cos49t  4^|-eln49t  jg-sin49t 


49 . Isin49t 


cos49t  -  ^~sin49t 


the  associated  C  matrix  of  equation  (4)  being  given  by 


with  hi  =  0  and  h2  =  h 


e^1(cos49h  -  sin49h)  49  s^n4^ 


For  values  of  h  ^  2nTT,  the  matrix  C  is  non  singular  and  we  may  compute  (C^C) 
-1  ^ 

C  directly  as  .  - 1 


49cot (49h)  -  2 


-49e  21lcsc49h 


Using  equation  (5)  one  obtains  the  required  state  reconstruction  equation  for 
x2(t)  in  terms  of  the  measurements  xjCt),  x^Ct-h)  and  u(s),  k-h<s^t. 

x2(t)  =  |49cot(49h)  -  I]  Xl(t)  -  (49e-2hcsc49h)Xl(t-h) 

— ip  (10) 

+  2406  e2hcos49h  -  2e2hsin49h  (-e2s (sin49s)u(t+s)ds 


+  2406  -e2hsin49h  Ue2scos49s  +  2e2ssin49s)u(t+s) ds 


The  implementation  of  this  state  reconstruction  algorithm  was  carried  out  on  an 
8080  microprocessor  with  a  delay  value  h  —  .01  sec.  The  position  output  state 
was  sampled  at  2.2  milisecond  intervals  and  the  accuracy  of  the  A/D  and  D/A 
converters  was  12  binary  bits. 
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+ sr 
^  a 


Recjonstikiction  with '  80*80/8085 
j  ||l  Micfco  IPrbcess’oi: 


Figure  4a  compares  the  actual  tach  output  signal  representing  the  X2(t)  state 
with  the  microprocessor  output  signal  which  attempts  to  reconstruct  X2(t)  via 
equation  (10)  using  only  the  first  two  terms  of  this  expression.  Note:  In 
this  case  equal  weighings  must  be  used  for  x^(t)  and  x^(t-h).  The  effects 
of  measurement  noise  are  readily  apparent  in  this  figure.  Figure  4b  again 
compares  measured  tach  output  with  the  microprocessor  output  signal,  however, 
in  this  case  the  full  state  reconstruction  equation  (10)  is  implemented.  This 
implementation  is  seen  to  give  a  very  accurate  state  reconstruction  which  is 
less  sensitive  to  measurement  noise. 

III.  FREQUENCY  DOMAIN  CONTROL  SYNTHESIS  USING  DELAY  FEEDBACK:  Several 
papers,  (see  Reference  (5  )  and  (6  ))  have  appeared  in  the  recent  literature 
which  address  the  problem  of  feedback  control  using  delays.  Reference  (6  ) 
develops  several  feedback  control  laws  using  delays  in  the  state  and  derivative 
of  the  state  which  are  shown  to  drive  the  full  state  of  the  system  to  zero 
and  keep  it  there.  The  constructions,  however,  have  limited  utility  in  servo 
control  applications  since  they  assume  first  that  the  control  space  has  the 
same  dimenstion  as  the  state  and  all  states  of  the  system  are  accessible  for 
on-line  measurement. 

In  this  section  we  consider  a  restricted  class  of  delay  feedback  controllers 
shown  in  Figure  5.  This  configuration  has  proved  quite  useful  in  turret  and 
servo  control  applications  in  which  G(s)  represents  the  open  loop  transfer 
function  between  the  command  input  and  the  position  output.  The  two  design 
parameters  introduced  by  delayed  feedback  are  seen  to  be  K,  the  feedback  gain, 

and  T,  the  feedback  time  delay.  The  reason  for  chosing  the  two  feedback  gains 

in  the  form  K  and  K-l  differing  by  unity  in  the  general  case,  will  be  made  clear 
below.  The  equivalent  feedback  transfer  function,  H(s),  for  the  system  in 
Figure  is: 

H(s)  =  K  -  (K-l)e“Ts  (11) 

We  may  represent  the  e“Ts  term  by  its  equivalent  Taylor  series  form  as: 

e-Ts  -  1  -  Ts  +  T2s2  -  T3s3  +  ... 

2  6 

The  frequency  band  of  primary  interest  from  a  stability  and  transient  response 
point  of  view  in  |5:  |G(s)H(s)|  F]  or  [s:  s— w3  wliere  wc  denotes  the 

gain  crossover  frequency  of  the  compensated  open  loop  system.  Setting  s  =  jw 

and  assumming  jwTl^^l,  we  may  approximate  e~^s  by  the  first  two  terms  of  its 

Taylor  series  expansion  or; 

e~Ts  =  1  -  Ts  =  1  -  jwT  (12) 

Substituting  (12)  into  equation  (11)  yeilds; 

H(s)  =  K  -  (K-l) ( 1-jwT)  =  1  +  j (K-l )Tw  (13) 

Since  K>1  will  be  required,  this  corresponds  to  a  phase  lead  network  (on  a 
first  order  approximation  basis)  in  the  controller  feedback  path.  If  this  phase 
lead  term  is  properly  positioned  in  frequency,  it  will  produce  a  stabilizing 
effect  upon  the  control  systems  unit  step  or  impulse  response  characteristics. 

As  will  be  seen  in  the  examples,  the  time  delay  or  feedback  gain  can  be  adjusted 
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to  provide  any  desired  damping  in  the  system  response.  The  procedure  for  intro¬ 
ducing  a  lead  network  effect  around  w  =  wc  using  delayed  feedback  can  now  be 
developed. 

First,  choose  wc  such  that  j G ( j wc ) |  =  1 

XT 

Second,  select  the  feedback  time  delay,  T,  such  that  Twc<^l.  The  choice  of 
Twc  =  .1  is  reasonable  and  is  used  in  the  examples.  For  this  choice,  the  first 
term  disregarded  in  the  series  expansion  has  magnitude  .005  at  wc  and  rapidly 
becomes  smaller  for  higher  frequencies.  Third,  select  the  feedback  gain  para¬ 
meter,  K,  such  that  the  lead  time  constant  becomes  effective  at  or  near  w  =  wc 
i.e.  (K-1)T'=  1  .  Note  under  this  condition  using  step  1  and  equation  (8),  that; 

wc 

G(wc)H(wc)|  =  | G(wc) |  j  H(wc) 

^ _ XT 

=  1  and  K  =  11 

The  delayed  feedback  design  procedure  thus  is  seen  to  be  straight  forward  in 
concept.  The  effect  of  the  particular  delayed  feedback  configuration  discussed 
here  is  to  replace  the  more  standard  tach  feedback  stabilization  loop.  When 
the  delay  time  and  gain  parameters  are  properly  chosen,  system  response  charact¬ 
eristics  may  be  improved  substantially. 

Example: 

We  consider  first  a  simple  laboratory  servo  system. whose  open  loop  transfer 
function,  G(s),  is  given  by; 

G(s)  =  600 

s(l+s)  (14) 

4 

The  -3db  crossover  frequency,  wc,  of  the  open  loop  transfer  function  G(s)  is 
56  rad/sec  and  the  delay  time,  T,  is  computed  from  step  3  and  satisfies  10T  =  1 
or  T  =  .0017  sec.  The  gain  K  is  fixed  and  satisfies  the  relation;  56 

K  -  1  =  1  =  1  =  10 

Twc  .1 

Due  to  limitations  of  the  8080  microprocessor,  the  above  design  using  a  delay  of 
1.7  ms  could  not  be  implemented.  The  smallest  delay  which  could  be  implemented 
with  the  8080  was  2.2  ms.  The  performance  of  this  design  for  a  step  input  command 
is  shown  in  Figure  6b.  Figure  6  illustrates  that  the  effective  damping  intro¬ 
duced  by  the  feedback  delay  can  be  further  increased  by  increasing  the  delay 
parameter  T.  The  desired  damping  can  also  be  achieved  by  adjusting  the  gain 
parameter  K. 

To  evaluate  the  effects  of  delay  parameters  which  were  too  small  for  imple¬ 
mentation  on  an  8080  microprocessor,  simulations  were  run  for  values  of  T  =  .8 
msec,  1.7  msec,  2.2  msec,  4.4  msec,  8.8  msec  and  17.6  msec,  using  the  servo 
transfer  function  (14).  These  results  are  shown  in  Figures  7  -  12.  Deficiencies 
in  the  linear  model  of  the  servo  system  are  readily  apparent  since  the  simulations 
indicate  more  damping  than  is  evident  in  the  test  results  of  Figure  6  and  Figure  12 
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indicates  an  instability  with  the  17.6  msec  delay  in  constraht  to  the  over 
damped  response  in  the  hardware  test  shown  in  Figure  6e. 

Example: 

In  this  example  we  illustrate  the  application  of  the  delay  feedback  control 
synthesis  technique  to  the  design  of  a  controller  for  an  XM97  helicopter  turret 
control  system  shown  in  Figure  13.  The  transfer  function  block  diagram  of  this 
system  is  shown  in  Figure  14.  The  -3db  crossover  frequency  for  the  open  loop 
system  (tach  loop  open)  was  computed  to  be  20  rad/sec  resulting  in  a  feedback 
time  delay  of  .003  sec.  The  step  response  of  the  original  XM97  design  is 
shown  in  Figure  15  and  that  of  the  delay  feedback  design  in  Figuia  16.  The 
latter  design  exhibits  a  dramatic  improvement  with  respect  to  overshoot  and 
settling  time.  This  improvement  can  be  explained  partially  by  the  fact  that  the 
original  system  uses  motor  tachometer  feedback  for  stabilization  while  the  delay 
feedback  design  effectively  uses  actual  turret  position  and  rate  for  feedback 
stabilization.  Figures  17  -  22  also  show  the  effects  of  increasing  and  decreasing 
the  delay  feedback  parameter.  Saturation,  columb  friction  and  deadband  non- 
linearities  are  included  in  the  simulation. 

_ CONCLUSION :  Applications  of  delay  feedback  for  state  construction 

and  feedback  control  are  presented  together  with  simulation  results  and  examples 
of  actual  implementations  using  Intel  8080  and  8085  microprocessors.  These 
examples  demonstrate  the  practicality  of  the  ideas  and  suggest  that  these  tech¬ 
niques  may  provide  a  useful  adjunct  to  the  more  standard  frequency  domain  and 
state  variable  techniques  for  estimation  and  control  applications. 
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